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1.  Introduction 


1.1  Background 

This  work  involves  the  development  of  a  delay-throughput  performance 
evaluator  for  distributed  systems.  Currently  planned  and  future  Navy  distributed 
integrated  computer  and  communications  systems  involve  the  extensive  use  of 
medium  access  control  procedures  for  sharing  distributed  communications, 
processing  and  computing  resources  among  distributed  stations.  This  work 
contributes  to  the  development  of  methods  and  tools  for  carrying  out  modeling, 
performance  evaluation,  analysis  and  design  of  such  systems. 

The  principal  area  of  work  for  this  effort  involves  research  and  development  of 
delay-throughput  performance  evaluation  tools  for  distributed  communications  and 
computer  multiple-access  schemes.  This  development  involves  the  modeling  and 
analysis  of  generic  TDMA  and  token-ring  type  medium  access  control  architectures 
that  are  of  key  importance  to  present  and  future  Navy  distributed  systems.  t 

i  v  ’>s 

It  has  been  well  demonstrated  that  the  performance  features  of  distributed 
computer  systems  and  their  associated  distributed  computer  communications 
networks  are  not  separable.  In  considering  distributed  systems  and  networks  that 
are  essential  to  the  execution  of  mission  critical  tasks,  it  is  essential  to  be  able  to 
evaluate  the  message/task  delay  and  system  throughput  performance  under  a 
military  priority  based  traffic  and  messaging  environment.  The  performance 
evaluation  tools  to  be  developed  under  this  effort  are  to  be  integrated  and  used  as 
key  elements  of  the  POD  (performance  oriented  design)  system.  They  would  provide 
efficient  tools,  based  upon  academically  developed  analytical  and  combined 
analytical/simulation  approaches,  for  the  evaluation  of  key  generic  medium  access 
control  procedures  which  are  critical  ingredients  of  Navy  distributed  computer  and 
communications  systems  and  networks.  They  will  allow  the  users  and  designers  of 
such  distributed  systems  to  assess  the  delay-throughput  performance  behavior  of 
their  systems  under  various  traffic  loading  and  system  and  channel  parameter 
conditions. 

Our  works  as  outlined  in  this  report  incorporate  the  following  elements. 

a.  Development  of  performance  evaluation  models  and  delay-throughput 
analysis  techniques  for  distributed  systems  involving  TDMA  and  token-ring  multiple- 
access  procedures;  with  applications  to  Navy  systems,  including:  MILSTD1553B, 
SAFENET II  and  JTIDS. 
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b.  Development  of  analytical  and  joint  analytical/simulation  programs  for  the 
multiple-access  systems  investigated  in  (a)  for  the  calculation  of  their  delay- 
throughput  performance. 

c.  The  capability  to  carry  out  performance  evaluation  tradeoffs  and  analyses 
using  the  programs  developed  in  (b). 

1.2  The  Delav-Throuohout  Performance  Evaluation  Programs 

The  delay  performance  evaluation  programs  developed  are  classified  as 
follows. 

A.  Time-Division  Multiple-Access  (TDMA)  models,  which  are  further 
categorized  into  three  separate  models  and  programs: 

A.1  PS-TDMA  (Packet-Switched  TDMA)  performance  evaluation 
program,  for  a  TDMA  station  which  shares  a  communications  (or  processing, 
computing,  buffering,  etc.)  channel  on  a  PS-TDMA  basis,  so  that  it  transmits  its 
packets  in  a  packet-switched  fashion  during  its  dedicated  time  slots.  We  obtain  the 
system  throughput  and  the  average,  variance  and  distributions  of  the  message  delay 
and  the  station  queue-size,  using  both  simulation  and  analytical  techniques.  This 
program  can  also  been  used  to  evaluate  the  performance  of  a  demand-assigned 
PS-TDMA  system. 

A.2  CS-TDMA  (Circuit-Switched  TPM  A)  performance  evaluation 
program,  for  a  TDMA  station  which  shares  a  communications  (or  processing, 
computing,  buffering,  etc.)  channel  on  a  CS-TDMA  basis,  so  that  it  transmits  its 
established  sessions  (connections)  in  a  circuit-switched  fashion  during  its  edicated 
time  slots.  We  obtain  the  system  throughput  and  the  average,  variance  and 
distributions  of  the  message  delay  and  the  station  queue-size,  using  both  simulation 
and  analytical  techniques.  This  program  can  also  been  used  to  evaluate  the 
performance  of  a  demand-assigned  CS-TDMA  system. 

A.3  l-TDMA  (integrated  CS/PS  TDMA)  performance  evaluation  program, 
for  a  TDMA  station  which  shares  a  communications  (or  processing,  computing, 
buffering,  etc.)  channel  on  a  multiplexed  integrated  PS-TDMA  and  CS-TDMA  basis, 
so  that  it  transmits  its  established  sessions  (connections)  in  a  circuit-switched  fashion 
during  its  dedicated  CS  time  slots  and  transmits  its  packets  on  a  PS-TDMA  during  its 
dedicated  PS  slots,  as  well  as  during  the  CS  frame  slots  which  remain  unused.  We 
obtain  the  system  throughput  for  the  PS  and  CS  components,  the  session  blocking 
probability  for  the  CS  connections,  and  the  average,  variance  and  distributions  of  the 
packet  delay  and  the  station  queue-size  for  the  PS  subsystem,  using  both  simulation 
and  analytical  techniques.  This  model  also  represents  the  delay-throughput 
performance  of  a  demand-assigned  integrated  CS/PS  services  TDMA  system. 
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B.  Timed  Token  Rotation  Protocol  (FDDI-I  Type)  model,  as  applied  to 
the  operation  of  a  token  ring  network.  This  protocol  is  of  the  same  type  employed 
by  FDDI  and  the  Navy  SAFENET  II  networks.  It  also  allows  performance  analysis  for 
numerous  distributed  computer  and  communications  systems  that  employ  polling 
access  protocols  of  this  type  to  share  access  to  system  resources  (such  as  system 
buses,  communications  channels,  processors,  memory,  etc.) 

An  extensive  simulation  model  has  been  set-up  which  allows  the  investigation 
of  the  performance  of  such  systems  under  various  network  traffic  loading  conditions 
and  system  configurations.  Performance  results  are  obtained  for  the  system 
throughput  and  for  message  delays  and  buffer  queue-sizes. 

The  presented  programs  constitute  the  Version  1  of  the  corresponding 
developments.  Further  extensions  and  expansions  are  being  planned. 
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2.  Delay-Throughput  Performance  Evaluator  for  Packet-Switched  Time- 
Division  Multiple-Access  (PS-TDMA)  Distributed  Systems 

2.1  Introduction 

In  this  Chapter,  we  describe  the  structure  and  provide  the  operation  directions 
for  the  Packet-Switched  Time  Division  Multiple-Access  (PS-TDMA)  program.  In 
Section  2.2,  the  structures  of  the  service,  channel,  message  and  traffic  models  are 
described,  and  the  performance  measures  are  defined.  In  Section  2.3,  we  provide 
detailed  instructions  for  running  the  PS-TDMA  program.  Examples  and  the  source 
codes  are  given  in  the  appendices.  The  output  is  discussed  in  Section  2.4. 

2.2  The  Structure  of  the  PS-TDMA  Model 

The  Channel  Structure 

A  Time  Division  Multiple  Access  (TDMA)  channel  structure  is  assumed.  The 
(service,  processing  or  communications)  channel  is  shared  among  multiple  user 
stations,  on  a  synchronous  basis.  The  channel  time  is  divided  into  time  slots  (or, 
slots)  so  that 

x  =  duration  of  a  slot  [sec.]  « transmission  time  of  a  single  packet. 

Time  frames  (or  simply  frames,  or  cycles)  are  consecutively  identified  across 
the  channel,  so  that 

A  time  frame  consists  m  slots. 

TF  =  duration  of  a  time  frame  =  mt. 

Each  station  is  dedicated  a  number  of  slots  within  each  time  frame,  during 
which  it  can  transmit  its  packets.  Since  each  station  is  then  operating  independently 
of  each  other  station,  and  its  performance  is  independent  of  that  of  any  other  station, 
we  need  to  examine  the  performance  of  only  an  arbitrary  single  station,  which  is  the 
station  subsequently  identified  here  (as  “the  station”). 

We  set 

The  number  of  slots  per  frame  dedicated  to  the  station  *  n; 
where:  nsm. 

The  channel  structure  is  illustrated  in  Fig.  1. 
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time  frame 
—  m=5  - 


time  frame 
—  m=5  — 


n=2 


n«2  — ► 


Fig.  1  Illustration  of  a  TDMA  channel  structure  whereby  a  station 
is  allocated  n«2  slots  during  each  time  frame  of  duration 
m=5  slots. 


The  Message  Structure 

Messages  arrive  at  random  at  the  station’s  buffer  (whose  capacity  is  not 
limited)  and  are  served  on  a  packet-switched  store-and-forward  FCFS  (First  come 
First  Served)  basis. 

A  message  contains  a  random  number  (say  P)  of  packets  with  an  average 
given  by: 

Average  number  of  packets  per  message  =  1/q, 
where  0<q£1 . 

Each  packet  is  assumed  to  be  a  segment  of  fixed  length,  where  the  packet 
transmission  time  is  equal  to  a  single  slot. 

Since  the  transmission  time  of  a  single  packet  is  equal  to  a  single  slot,  if  the 
message  would  have  been  transmitted  continuously  across  the  channel,  its 
transmission  rate  would  be 

message  transmission  rate  =  q  [messages/slot]. 

In  the  PSTDMA  program  it  is  further  assumed  that  the  number  of  packets  per 
message  (P)  is  governed  by  a  Geometric  distribution,  so  that  q  also  designates  the 
probability  that  a  packet  transmission  is  the  last  transmission  of  the  underlying 
packet  (i.e.,  with  probability  q  this  transmitted  packet  is  the  last  packet  and  with 
probability  1-q  the  message  contains  more  packets).  The  probability  that  a  message 
contains  k  packets  is  given  by  the  Geometric  distribution: 
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Prob{P=k}=(1  -q)k_1q,  k=1,2,3,... 


Example:  A  message  contains  an  average  of  5  packets.  Then,  q=1/5=0.2,  and 

the  transmission  rate  is  equal  to  5  [mess./slot].  If  the  transmission  rate  across  the 
channel  is  equal  tc  10  Kbps  and  each  packet  (of  fixed  size)  contains  1000  bits,  then 

the  packet  transmission  time=s!ot  duration=T=1 000/1 0K=0.1  sec. 

The  message  average  transmission  time  of  5  slots  would  take  0.5  sec  if  5 
consecutive  slots  are  allocated  to  the  transmission  of  this  message.  Otherwise,  a 
longer  transmission  time  would  be  required. 

The  Message  Arrival  Process 

Messages  are  assumed  to  arrive  into  the  station  in  accordance  with  a 
Geometric-Batch  arrival  point  process  described  as  follows.  Message  arrivals  are 
recorded  at  the  end  of  the  arrival  slot. 

A  batch  of  messages  will  arrive  in  a  slot  independently  of  arrival  into  other 
slots.  We  set: 

P{a  batch  of  messages  arrives  in  a  slot}=p; 

P{no  messages  arrive  in  a  slot}=1-p; 
where  0<p<1 . 

The  number  of  messages  in  an  arrival  batch  (say  B)  is  a  random  variable,  with 
an  average  that  is  set  to  be  given  by 

The  average  size  of  the  arrival  messages  batch  =  b  [mess./batch]. 

The  PSTDMA  program  developed  by  IRI  allows  the  user  to  select  one  of  three 
distributions  for  the  number  of  messages  in  the  batch: 

1.  A  Deterministic  distribution;  under  which  the  number  of  messages  in  the 
arrival  batch  is  always  fixed,  equal  to  b. 

2.  A  Geometric  distribution  with  mean  set  equal  to  b. 

3.  A  uniform  distribution,  whereby  the  number  of  messages  in  a  batch  is 
uniformly  distributed  between  a  specified  lower  value  Ui  and  upper  level  U2;  note 
the  mean  value  is  now  equal  to  (U]  +U2  )/2. 

Performance  Measures 

Ihroutthoul 

The  station’s  throughput  (THS)  is  equal  to  the  number  of  messages  per  unit 


8 


Delay-throughput  Evaluator,  IRI  Corp.,  Dr.  Izhak  Rubin 


time  (or  per  slot)  transmitted  by  the  station  across  the  channel.  Under  the  PS-TDMA 
model  the  station  throughput  rate  is  equal  to  the  message  arrival  rate  (no  messages 
are  blocked)  so  that  we  have 

THS  =  bp  [mess./slot]. 

The  packet  throughput  is  given  by 

Packet  throughput  =  bp/q  [packets/slot]. 

Note  that  since  each  slot  is  equal  to  x  [sec.],  where 

T[sec]  =(average  #bits/packet)/(channel  transmission  rate  in  bps), 
we  nave,  per  station, 

Station  Message  Throughput  =  bp/T  [mess/sec] 

Station  Packet  Throughput  =  bp/qf  [packets/sec]. 

The  station’s  normalized  throughput,  or  traffic  Intensity,  is  equal  to  the 
ratio  of  the  traffic  rate  of  arriving  packets  at  the  station  and  the  channel’s  service  rate 
dedicated  to  the  station,  and  is  thus  given  by 

Normalized  Throughput=p=(bp)/(nq/m). 

For  system  stability  (i.e,  to  ensure  finite  steady  state  system  delays  and  buffer  queue 
sizes)  it  is  necessary  to  ensure  that 

p<1. 


Message  Delay  and  Queue-Size 

The  message  waiting-time  is  defined  as: 

W=Message  Waiting  Time=Total  time  elapsed  since  the  message  arrival  slot 
to  the  time  that  the  first  packet  of  the  message  starts  transmission 

The  message  effective  transmission  time  is 

T-message  effective  transmission  time=Time  elapsed  since  the  start  of 
transmission  of  the  session’s  first  packet  to  the  time  that  the  last 
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packet  ends  transmission 
The  messge  delay  time  (D)  is  defined  as: 

D=Message  DeIay=Total  time  elapsed  since  the  message  arrival  slot  to  the 
time  that  the  last  packet  of  the  message  is  transmitted 


Thus,  we  have 


D=W+T. 


The  Message  queue-size  is: 

X=Message  queue-size=number  of  messages  resident  in  the  station’s  buffer 
(including  the  message  in  the  process  of  transmission,  if  any) 

The  PSTDMA  program  provides  results  for  the  mean  and  standard-deviation 
of  X  and  D  as  we!!  as  for  the  distributions 
P(j)=P{X=j):  D(j)=P{D=j). 

These  results  are  obtained  by  us  through  the  use  of  two  methods: 

a.  We  develop  a  simulation  program  that  employs  analytic  recurrence 
relationships  expressing  the  evolution  of  the  system  states.  Random  generators  are 
used  to  generate  the  traffic  loading.  The  sample  mean,  variance  and  distributions  of 
the  queue  size  and  delay  variables  are  then  obtained  through  statistics  collections 
and  computations.  The  user  is  able  to  select  the  number  of  simulation  runs  (slots)  for 
the  stop  time  as  well  as  the  start  time  for  the  collection  of  statistics  to  incorporate  in 
the  calculation  of  the  X  and  D  performance  measures. 

b.  We  use  methods  developed  by  us  and  presented  in  the  References  which 
provide  for  analytical  derivations  of  exact  and/or  lower  and  upper  bound  formulas  for 
the  steady  state  mean  queue  sizes  and  message  delays. 

Note  that  the  analytical  results  are  steady  state  results  and  thus  describe 
performance  results  after  the  system  has  been  running  for  a  long  time.  Hence,  the 
analytical  based  output  results  can  differ  from  simulation  results  which  are  based  on 
shorter  run  times. 
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2.3  Instructions  for  Running  the  PSTDMA  Program 

In  the  following  we  provide  instructions  for  running  the  PSTDMA  program 
(Version  1 ).  It  is  noted  that  the  specific  name  of  the  program  used  can  be  PSTDMAx 
where  x  is  a  number  designating  a  version  of  the  program  compiled  to  incorporate 
certain  limits  on  the  program  size.  For  example,  PSTDMA2  involves  the  limit  m£500, 
and  printed  values  for  the  queue-size  (x)  and  for  the  delay  (d)  which  are  in  the  range 
x^lOO,  d<100. 

Step  bv  Step  Instructions  (See  Appendix  A  for  examples) 

1 .  Enter  the  name  of  the  program  (such  as  PSTDMA2) 

Subsequent  inputs  are  entered  in  response  to  program  prompts. 

2.  Enter  the  name  of  an  output  file;  no  more  than  12  characters;  in  PC,  quote 
the  ‘file  name’ 

3.  Enter  a  number  to  designate  the  run  method.  Select 

1  -  for  simulation  only 

2  -  for  analysis  only 

3  -  for  simulation  and  analysis 

4.  Enter  the  start  time  for  simulation  collection  start 

5.  Enter  the  stop  time  for  the  simulation  length  duration 

6.  Enter  the  frame  duration  (m,  an  integer,  where  l£m<specified  upper  oound 
such  as  500);  then  also  type  (following  a  space  key  entry)  the  number  of  slots 
allocated  to  the  station  (n,  an  integer,  where  Knsm) 

7.  Enter  the  batch  arrival  rate  (p=p1 ,  so  that  0<p<1 ,  where  p  designates  the 
probability  that  a  batch  arrives  in  a  slot) 

8.  Enter  the  message  transmission  rate  (q=q1,  0<q£l,  where  q  is  the 
probability  that  the  transmitted  packet  is  the  last  one  belonging  to  the  underlying 
message,  and  1/q  is  the  average  number  of  packets  belonging  to  a  message,  or  the 
average  number  of  slots  required  to  transmit  a  message) 

9.  Enter  the  batch  size  distribution  index;  select 

1  -  for  deterministic  (fixed)  batch  size 

2  -  for  Geometric  f  *.ch  size 

3  -  for  uniform  be»ch  size 

10.  If  selected  deterministic  or  Geometric  batch  distributions,  enter  next  the 
mean  batch  size  b. 

If  selected  a  uniform  distribution,  select  the  lower  and  upper  levels  for 
the  batch  sizes  U]  and  U2;  the  program  then  ca.  jlates  the  mean  batch  size  as 
b=(U|+u2)/2. 

To  ensure  system  stability,  so  that  the  message  arrival  rate  is  lower  than 
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the  channel  message  service  rate,  it  is  necessary  to  ensure  that  the  normalized 
throughput  (rho=p)  is  less  than  one,  or  that 

bp<(nq/m). 

If  the  parameters  selected  violate  this  condition,  the  program  indicates  so  and 
requests  the  user  to  re-select  the  system  parameters. 

11.  When  simulation  is  used,  enter  the  levels  x  and  d  for  the  program  to 
explicitly  provide  at  the  output  the  probabilities 

P{X>x},  P{D>d}. 

12.  The  program  then  responds  with  ‘please  wait’  and  proceeds  with  its  run. 
When  done,  the  output  is  shown  on  the  screen  as  well  as  provided  in  the  specified 
output  file. 

Extension:  An  extended  version  PSTDMA5  includes  in  addition  to  the  above 
described  batch  message  arrival  model  also  a  Poisson  arrival  model.  The  structure 
of  the  program  is  similar.  The  user  is  asked  to  select  which  arrival  model  is  desired. 
If  the  batch  message  arrival  model  (selection  1)  is  chosen,  the  remainder  of  the  input 
is  as  described  above.  If  a  Poisson  arrival  process  is  chosen  (selection  2),  the  only 
arrival  parameter  that  must  be  selected  is  the  message  arrival  rate  A  expressed  in 
[mess/slot]  units.  Note  that  A=bp  represent  the  average  message  arrival  rate  per 
slot.  The  message  length  is  Geometrically  distributed  with  mean  1/q.  as  for  the  batch 
arrival  model.  The  program  proceeds  otherwise  as  described  above  and  in  the 
following. 
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2.4  Output  of  the  PSTDMA  Program 

Appendix  A  provides  two  examples  of  the  inputs  and  outputs  for  the  PSTDMA 
program. 

The  output  of  the  program  contains  the  following  information: 

1.  Statement  of  the  input  parameters. 

2.  The  normalized  throughput  (rho=p)  for  the  station  under  consideration,  in 
relation  to  the  transmission  time  provided  for  this  station 

3.  When  simulation  used,  a  table  presenting  the  queue-size  and  message 
delay  distributions,  P(j)  and  D(j),  respectively. 

4.  For  simulation:  the  means  and  standard  deviations  for  the  station  queue 
size  and  for  the  message  delay;  the  probabilities  P(X>x),  P(D>d),  for  the  selected  x 
and  d  values 

5.  For  analysis:  the  means  or  the  station  queue  size  and  for  the  message 
delay,  exact  and  upper  and  lower  bound  results. 

As  indicated  in  tne  References,  the  analytical  method  developed  for  the 
calculation  of  the  upper  and  lower  bounds  is  much  more  computationally  efficient 
than  the  method  developed  for  the  analytical  calculation  of  the  exact  result. 

Further  note  that  the  analytical  steac  ;tate  mean  delay  E(D)  and  mean 
queue  size  E(X)  are  related  through  Little's  Theorem  so  that 

E(X)=bpE(D). 
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3.  Delav-ThrouahDut  Performance  Evaluator  for  Circuit-Switched  Time- 
Division  Multiple-Access  fCS-TDMAl  Distributed  Systems 

3.1  Introduction 

In  this  Chapter,  we  describe  the  structure  and  provide  the  operation  directions 
for  the  Cacket-Switched  Time  Division  Multiple-Access  (CS-TDMA)  program.  In 
Section  3.2,  the  structures  of  the  service,  channel,  message  and  traffic  models  are 
described,  and  the  performance  measures  are  defined.  In  Section  3.2,  we  provide 
detailed  instructions  for  running  the  PS-TDMA  program.  Examples  and  the  source 
codes  are  given  in  the  appendices.  The  format  of  the  output  is  discussed  in  Section 
3.4. 

3.2  The  Structure  of  the  CS-TDMA  Model 

The  Channel  Structure 

A  Time  Division  Multiple  Access  (TDMA)  channel  structure  is  assumed.  The 
(service,  processing  or  communications)  channel  is  shared  among  multiple  user 
stations,  on  a  synchronous  basis.  The  channel  time  is  divided  into  time  slots  (or, 
slots)  so  that 

t  *  duration  of  a  slot  [sec.]  =  transmission  time  of  a  single  segment 

Time  frames  (or  simply  frames,  or  cycles)  are  consecutively  identified  across 
the  channel,  so  that 

A  time  frame  consists  m  slots. 

Tf  =  duration  of  a  time  frame  =  mx. 

Each  station  is  dedicated  a  number  of  slots  within  each  time  frame,  during 
which  it  can  transmit  its  established  session’s  segments.  Since  each  station  is  then 
operating  independently  of  each  other  station,  and  its  performance  is  independent  of 
that  of  any  other  station,  we  need  to  examine  the  performance  of  only  an  arbitrary 
single  station,  which  is  the  station  subsequently  identified  here  (as  “the  station”). 

We  set 

The  number  of  slots  per  frame  dedicated  to  the  station  =  n; 
where:  nsm. 

The  channel  structure  is  illustrated  in  Fig.  2. 
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time  frame 
—  m=5 


time  frame 
—  m=5  — 


n=3 - ► 


n=3 - ► 


Fig.  2  Illustration  of  a  TDMA  channel  structure  whereby  a  station 
is  allocated  n=3  slots  during  each  time  frame  of  duration 
m=5  slots. 


TheSession  Structure 

Sessions  (or  connection  requests)  arrive  at  random  at  the  station's  buffer 
(whose  capacity  is  not  limited)  and  are  served  on  a  circuit-switched  store-and- 
forward  FCFS  (First  come  First  Served)  basis. 

Each  session  (connection)  requires  for  its  support  across  the  channel  the 
allocation  of  a  single  slot  during  each  frame,  for  as  long  as  the  session  lasts  (i.e.,  for 
the  session’s  holding  time).  The  station  transmits  a  single  segment  of  information 
from  the  session's  output  during  each  session’s  designated  slot. 

Since  the  station  is  allocated  only  n  slots  per  frame,  it  can  simultaneously 
support  at  most  n  sessions.  If  a  larger  number  of  session  requests  for  support  arrive, 
they  will  be  queued  on  a  FCFS  basis,  and  served  when  other  sessions  are 
terminated.  Under  the  underlying  CSTDMA  model,  a  delayed  rather  then  blocking 
service  is  provided,  so  that  sessions  arriving  at  a  station  when  all  of  its  dedicated 
slots  are  occupied  would  not  be  blocked  (as  in  the  ITDMA  CS  component  model), 
but  rather  queued  until  served.  The  CSTDMA  program  provides  performance 
analysis  for  determining  this  session  delay. 

A  session  lasts  for  a  random  number  (say  C)  of  frames  with  an  average  given 
by: 


Average  session  holding  time 

■  average  number  of  frames  per  session  connection  =  1/q, 
where  0<q£l . 

The  session  transmission  rate  is  thus  also  denoted  as 
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session  transmission  rate  =  q  [sessions/frame]. 

In  the  CSTDMA  program  it  is  further  assumed  that  the  session  holding  time 
(number  of  frames  per  session)  is  governed  by  a  Geometric  distribution,  so  that  q 
also  designates  the  probability  that  a  segment  frame  transmission  is  the  last 
transmission  of  the  underlying  session(i.e.,  with  probability  q  this  transmitted 
segment  is  the  last  segment  and  with  probability  1-q  the  session  contains  more 
segments).  The  probability  that  a  session  lats  for  k  frames  is  given  by  the  Geometric 
distribution: 


Prob{C=k}=(1  -q)k_1q,  k=1,2,3,... 

Example:  A  session  average  holding  time  is  equal  to  500  frames.  Then, 

q=1/500=0.002,  and  the  transmission  rate  is  equal  to  q=500  [frames/session].  If  the 
transmission  rate  across  the  channel  is  equal  to  R=640  Kbps  and  each  segment  (of 
fixed  size)  contains  1 00  bits,  then 

the  segment  transmission  time=slot  duration=T=100/640K=0. 15625  msec. 
Assume  that  each  session  requires  a  channel  circuit  operating  at  an  average  rate  of 
16  Kbps.  Then,  the  number  of  frame  slots  must  be  set  equal  to 
m=channel  rate/circuit  rate=640K/1 6K=40. 

The  frame  thus  contains  m=40  slots,  and  the  frame  duration  is  then  equal  to 

TF=rm:=6.25  msec. 

The  session’s  average  holding  time  is  equal  to  50  frames,  so  that  the  average 
session  holding  time  is  equal  to  500x  6.25  msec=3.125  sec. 

IheSessipn  Arrival  Process 

Sessions  (i.e.,  requests  for  allocation  of  a  circuit  to  establish  a  session)  are 
assumed  to  arrive  into  the  station  in  accordance  with  a  Geometric-Batch  arrival  point 
process  described  as  follows.  Session  arrivals  are  recorded  at  the  end  of  the  arrival 
slot. 

A  batch  of  sessions  will  arrive  in  a  slot  independently  of  arrivals  into  other 
slots.  We  set: 

P{a  batch  of  sessions  arrives  in  a  slot}=p; 

P{no  sessions  arrive  in  a  slot}=1-p; 
where  0<p<1 . 

The  number  of  sessions  in  an  arrival  batch  (say  B)  is  a  random  variable,  with 
an  average  that  is  set  to  be  given  by 

The  average  size  of  the  arrival  sessions  batch  =  b  [sessions./batch]. 
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The  CSTDMA  program  developed  by  IRI  allows  the  user  to  select  one  of  three 
distributions  for  the  number  of  sessions  in  the  batch: 

1.  A  Deterministic  distribution;  under  which  the  number  of  sessions  in  the 
arrival  batch  is  always  fixed,  equal  to  b. 

2.  A  Geometric  distribution  with  mean  set  equal  to  b. 

3.  A  uniform  distribution,  whereby  the  number  of  sessions  in  a  batch  is 
uniformly  distributed  between  a  specified  lower  value  Uj  and  upper  level  U2;  note 
the  mean  value  is  now  equal  to  (Uj  +U2  )/2. 

Performance  Measures 

Throughput 

The  station’s  throughput  (THg)  is  equal  to  the  number  of  sessions  per  unit 
time  (or  per  slot)  transmitted  by  the  station  across  the  channel.  Under  the  CS-TDMA 
model  the  station  throughput  rate  is  equal  to  the  session  arrival  rate  (no  sessions  are 
blocked)  so  that  we  have 

THS  =  bp  [sess./slot]. 

Note  that  since  each  slot  is  equal  to  t  [sec.],  where 

'ctsec]  =(average  #bits/segment)/(channel  transmission  rate  in  bps), 
we  have,  per  station, 

Station  Message  Throughput  *  bp/T  [sess/sec]. 

The  station’s  normalized  throughput,  or  traffic  Intensity,  is  equal  to  the 
ratio  of  the  traffic  rate  of  arriving  sessions  at  the  station  and  the  channel’s  session 
service  rate  dedicated  to  the  station,  and  is  thus  given  by 

Normalized  Throughput=p«(bp)/(nq/m). 

For  system  stability  (i.e,  to  ensure  finite  steady  state  system  delays  and  buffer  queue 
sizes)  it  is  necessary  to  ensure  that 

p<1. 
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Message  Delay  and  Queue-Size 

The  session  waiting-time  is  defined  as: 

W=Session  Waiting  Time»Total  time  elapsed  since  the  session’s  arrival  slot  to 
the  time  that  the  first  segment  of  the  session  starts  transmission 

The  session  transmission  time  is 

T=Session  Transmission  Time=Time  elapsed  since  the  start  of  transmission  of 
the  session’s  first  segment  to  the  time  that  the  last  segment  ends 
transmission 

The  session  delay  time  (D)  is  defined  as: 

D=Session  Delay=Total  time  elapsed  since  the  session’s  arrival  slot  to  the 
time  that  the  last  segment  of  the  session  ends 
transmission 


Thus,  we  have 


D=W+T. 


The  station’s  queue-size  is: 

X=station  queue-size=number  of  sessions  (or  session  requests)  resident  in 
the  station’s  buffer  (including  sessions  to  which  a  circuit  is  currently 
allocated,  if  any) 

The  CSTDMA  program  provides  results  for  the  mean  and  standard-deviation 
of  X  and  D  as  well  as  for  the  distributions 
P(j)=P{X»j);  D(j)=P{D=j). 


These  results  are  obtained  by  us  through  the  use  of  two  methods: 
a.  We  develop  a  simulation  program  that  employs  analytic  recurrence 
relationships  expressing  the  evolution  of  the  system  states.  Random  generators  are 
used  to  generate  the  traffic  loading.  The  sample  mean,  variance  and  distributions  of 
the  queue  size  and  delay  variables  are  then  obtained  through  statistics  collections 
and  computations.  The  user  is  able  to  select  the  number  of  simulation  runs  (slots)  for 
the  stop  time  as  well  as  the  start  time  for  the  collection  of  statistics  to  incorporate  in 


18 


Delay-throughput  Evaluator,  IRI  Corp.,  Dr.  Izhak  Rubin 


the  calculation  of  the  X  and  D  performance  measures. 

b.  We  use  methods  developed  by  us  and  presented  in  the  References  which 
provide  for  analytical  derivations  of  exact  and/or  lower  and  upper  bound  formulas  for 
the  steady  state  mean  queue  sizes  and  message  delays. 

Note  that  the  analytical  results  are  steady  state  results  and  thus  describe 
performance  results  after  the  system  has  been  running  for  a  long  time.  Hence,  the 
analytical  based  output  results  can  differ  from  simulation  results  which  are  based  on 
shorter  run  times. 
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3.3  Instructions  for  Running  the  CSTDMA  Program 

In  the  following  we  provide  instructions  for  running  the  CSTDMA  program 
(Version  1).  It  is  noted  that  the  specific  name  of  the  program  used  can  be  CSTDMAx 
where  x  is  a  number  designating  a  version  of  the  program  compiled  to  incorporate 
certain  limits  on  the  program  size.  For  example,  CSTDMA2  involves  the  limit  m£500, 
and  printed  values  for  the  queue-size  (x)  and  for  the  delay  (d)  which  are  in  the  range 
x^lOO,  d£100. 

Step  bv  Step  Instructions  (See  Appendix  C  for  examples) 

1 .  Enter  the  name  of  the  program  (such  as  CSTDMA2) 

Subsequent  inputs  are  entered  in  response  to  program  prompts. 

2.  Enter  the  name  of  an  output  file;  no  more  than  12  characters;  in  PC,  quote 
the  ‘file  name’ 

3.  Enter  a  number  to  designate  the  run  method.  Select 

1  -  for  simulation  only 

2  -  for  analysis  only 

3  -  for  simulation  and  analysis 

4.  Enter  the  start  time  for  simulation  collection  start 

5.  Enter  the  stop  time  for  the  simulation  length  duration 

6.  Enter  the  frame  duration  (m,  an  integer,  where  ism^specified  upper  bound 
such  as  500);  then  also  type  (following  a  space  key  entry)  the  number  of  slots 
allocated  to  the  station  (n,  an  integer,  where  l£n£m) 

7.  Enter  the  batch  arrival  rate  (p=p1 ,  so  that  0<p<1 ,  where  p  designates  the 
probability  that  a  batch  arrives  in  a  slot) 

8.  Enter  the  message  transmission  rate  (q=q1 ,  0<q£l ,  where  q  is  the 
probability  that  the  transmitted  segment  is  the  last  one  belonging  to  the  underlying 
message,  and  1/q  is  the  average  session  holding  time,  measured  in  frames) 

9.  Enter  the  batch  size  distribution  (for  the  number  of  sessions  arriving  in  a 
batch)  index;  select 

1  -  for  deterministic  (fixed)  batch  size 

2  -  for  Geometric  batch  size 

3  -  for  uniform  batch  size 

10.  If  selected  deterministic  or  Geometric  batch  distributions,  enter  next  the 
mean  batch  size  b. 

If  selected  a  uniform  distribution,  select  the  lower  and  upper  levels  for 
the  batch  sizes  Uj  and  u2;  the  program  then  calculates  the  mean  batch  size  as 
b=(Ui+u2)/2. 

To  ensure  system  stability,  so  that  the  session  arrival  rate  is  lower  than 
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the  channel's  session  service  rate  for  the  underlying  station,  it  is  necessary  to  ensure 
that  the  normalized  throughput  (rho=p)  is  less  than  one,  or  that 

bp<(nq/m). 

If  the  parameters  selected  violate  this  condition,  the  program  indicates  so  and 
requests  the  user  to  re-select  the  system  parameters. 

1 1 .  When  simulation  is  used,  enter  the  levels  x  and  d  for  the  program  to 
explicitly  provide  at  the  output  the  probabilities 

P{X>x},  P{D>d}. 

12.  The  program  then  responds  with  ‘please  wait’  and  proceeds  with  its  run. 
When  done,  the  output  is  shown  on  the  screen  as  well  as  provided  in  the  specified 
output  file. 

Extension:  An  extended  version  CSTDMA5  includes  in  addition  to  the  above 
described  batch  session  arrival  model  also  a  Poisson  arrival  model.  The  structure  of 
the  program  is  similar.  The  user  is  asked  to  select  which  arrival  model  is  desired.  If 
the  batch  session  arrival  model  (selection  1)  is  chosen,  the  remainder  of  the  input  is 
as  described  above.  If  a  Poisson  arrival  process  is  chosen  (selection  2),  the  only 
arrival  parameter  that  must  be  selected  is  the  session  arrival  rate  A  expressed  in 
[mess/slot]  units.  Note  that  A=bp  represent  the  average  session  arrival  rate  per  slot. 
The  session  length  is  Geometrically  distributed  with  mean  1/q.  as  for  the  batch  arrival 
model.  The  program  proceeds  otherwise  as  described  above  and  in  the  following. 
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3.4  Output  of  the  CSTDMA  Program 

Appendix  C  provides  two  examples  of  the  inputs  and  outputs  for  the  CSTDMA 
program. 

The  output  of  the  program  contains  the  following  information: 

1 .  Statement  of  the  input  parameters. 

2.  The  normalized  throughput  (rho=p)  for  the  station  under  consideration,  in 
relation  to  the  transmission  time  provided  for  this  station 

3.  When  simulation  used,  a  table  presenting  the  queue-size  and  session 
delay  distributions,  P(j)  and  D(j),  respectively. 

4.  For  simulation:  the  means  and  standard  deviations  for  the  station  queue 
size  and  for  the  session  delay;  the  probabilities  P(X>x),  P(D>d),  for  the  selected  x 
and  d  values 

5.  For  analysis:  upper  and  lower  bounds  for  the  means  of  the  station  queue 
size  and  for  the  session  delay 

As  indicated  in  the  References,  the  analytical  method  developed  for  the 
calculation  of  the  upper  and  lower  bounds  is  much  more  computationally  efficient 
than  the  method  developed  for  the  analytical  calculation  of  the  exact  result. 

Further  note  that  the  analytical  steady-state  mean  delay  E(D)  and  mean 
queue  size  E(X)  are  related  through  Little's  Theorem  so  that 

E(X)=bpE(D). 
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4.  Delav-Throuahput  Performance  Evaluator  for  Integrated  Circuit- 
Switched  and  Packet-Switched  Time-Division  Multiple-Access  flTDMAl 
Distributed  Systems 

4.1  Introduction 

In  this  Chapter,  we  describe  the  structure  and  provide  the  operation  directions 
for  the  Integrated  Circuit-Switched  and  Packet-Switched  Time  Division  Multiple- 
Access  (ITDMA)  program.  In  Section  4.2,  the  structures  of  the  service,  channel, 
message  and  traffic  models  are  described,  and  the  performance  measures  are 
defined.  In  Section  4.3,  we  provide  detailed  instructions  for  running  the  ITDMA 
program.  Examples  and  the  source  codes  are  given  in  the  appendices.  The  output 
is  discussed  in  Section  4.4. 

4.2  The  Structure  of  the  ITDMA  Model 

The  Channel  Structure 

A  Time  Division  multiplexing  channel  structure  is  assumed.  The  (service, 
processing  or  communications)  channel  is  shared  among  the  multiple  user  Cs 
sessions  and  PS  packets,  on  a  synchronous  basis. 

A  single  station  is  considered.  This  station  is  allocated  the  global  channel 
bandwidth  and  time  for  the  transmission  of  its  information.  This  station  can  serve  as 
the  network  control  station  of  an  integrated  CS/PS  demand-assignment/TDMA 
(DA/TDMA),  which  is  the  basis  for  many  military  and  commercial  system  and  network 
operations. 

The  station  accommodates  two  traffic  types: 

1.  Sessions  which  are  accommodated  on  a  circuit-switched  (CS)  basis;  a 
supported  CS  session  (connection)  requires  the  allocation  of  a  single  slot  per  frame 
for  the  duration  of  the  session;  during  a  slot  the  session  transmits  a  segment 

2.  Messages  which  are  served  on  a  packet-switched  (PS)  FCFS  basis.  Each 
message  is  assumed  here  to  consist  of  a  single  packet.  The  packet  transmission 
time  is  equal  to  the  slot  duration. 

The  channel  time  is  divided  into  time  slots  (or,  slots)  so  that 

t  e  duration  of  a  slot  (sec.]  =  transmission  time  of  a  single  PS-packet. 

Time  frames  (or  simply  frames,  or  cycles)  are  consecutively  identified  across 
the  channel,  so  that 

A  time  frame  consists  m  slots. 
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Tf  =  duration  of  a  time  frame  =  mt. 

The  station  provides  constrained  higher  priority  to  the  support  of  CS  sessions, 
by  allocating  within  each  frame  a  maximum  of  n  slots  which  are  used  for  the  support 
of  CS  sessions,  so  that 

the  maximum  number  of  slots  per  frame  dedicated 

to  the  support  of  CS  sessions  =  n; 
where:  n^m. 

The  remainder  of  the  slots  in  the  frame,  m-n  slots,  are  always  dedicated  to  the 
transmission  of  PS-packets.  However,  in  addition,  any  slots  within  the  frame  which 
are  not  utilized  for  the  support  of  CS-sessions  are  also  made  available  for  the 
transmission  of  PS-packets. 

The  channel  structure  is  illustrated  in  Fig.  3. 


time  frame  time  frame 

^  . -  m=5  ►  ^  m*=5  —  ^ 

CS  PS  CS  PS 


n=3 


n-3 


Fig.  3  Illustration  of  a  Integrated  CS-PS  TDMA  channel  structure 
where  a  maximum  of  n=3  slots  during  each  frame  of  duration 
m«5  slots  are  allocated  to  CS  traffic;  during  the  first  frame 
all  3  CS  slots  are  occupied  by  3  CS  sessions;  during  the  second 
frame,  only  1  CS  slot  is  occupied,  and  the  remaining  4 
frame  slots  are  used  for  PS-packet  transmissions 


TheSession  Structure 

Sessions  (or  connection  requests)  arrive  at  random  at  the  station’s  buffer 
(whose  capacity  is  not  limited)  and  are  served  on  a  circuit-switched  store-and- 
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forward  FCFS  (First  como  First  Served)  basis  using  the  n  slots  per  frame  allocated  to 
the  support  of  CS  session  traffic. 

Each  session  (connection)  requires  for  its  support  across  the  channel  the 
allocation  of  a  single  slot  during  each  frame,  for  as  long  as  the  session  lasts  (i.e.,  for 
the  session’s  holding  time).  The  station  transmits  a  single  segment  of  information 
from  the  session’s  output  during  each  session’s  designated  slot. 

Since  the  station  is  allocated  only  n  slots  per  frame,  it  can  simultaneously 
support  at  most  n  sessions.  A  blocking  CS  service  is  assumed,  so  that  a  session 
arriving  whem  all  the  CS  slots  are  allocated  will  be  blocked  and  assumed  lost.  The 
ITDMA  program  uses  a  Poisson/Geom/n/n  queueing  system  model  to  carry  out  the 
analysis  of  the  CS  session  support  subsystem.  In  particular,  the  following  functions 
and  probabilities  are  computed: 

v,=P{i  CS  slots  per  frame  are  occupied  by  sessions},  0£i£n; 

V=average  number  of  slots  per  frame  occupied  by  CS  session 
connectios; 

PB=blocking  pt  obability=probabili**'  *’  ~t  a  session  is  blocked. 

A  session  lasts  for  a  rand  -n  number  (say  C)  of  frames  with  an  average  given 
by: 


Average  session  holding  time 

=  average  number  of  frames  per  session  connection  =  1/q, 
where  0<q<1 . 

The  session  transmission  rate  is  thus  also  denoted  as 
session  transmission  rate  =  q  [sessions/frame]. 

In  the  ITDMA  program  it  is  further  assumed  that  the  session  holding  time 
(number  of  frames  per  session)  is  governed  by  a  Geometric  distribution,  so  that  q 
also  designates  the  probability  that  a  segment  frame  transmission  is  the  last 
transmission  of  the  underlying  session(i.e.,  with  probability  q  this  transmitted 
segment  is  the  last  segment  and  with  probability  1-q  the  session  contains  more 
segments).  The  probability  that  a  session  lats  for  k  frames  is  given  by  the  Geometric 
distribution: 


Prob{C=k}=(1  -q)k_,q,  k=1,2,3,... 
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TheSession  Arrival  Process 

Sessions  (i.e.,  requests  for  allocation  of  a  circuit  to  establish  a  session)  are 
assumed  to  arrive  into  the  station  in  accordance  with  a  Poisson  arrival  point  process 
characterized  by  the  arrival  rate  A 

A  =  Session  arrival  rate  =  Avearge  number  of  session  arrivals  per  slot. 

The  PS  Message  and  PacketStructure 

Messages  arrive  at  random  at  the  station’s  buffer  (whose  capacity  is  not 
limited)  and  are  served  on  a  packet-switched  store-and-forward  FCFS  (First  come 
First  Served)  basis. 

A  message  is  assumed  here  to  consist  of  a  single  packe. 

Each  packet  is  assumed  to  be  a  segment  of  fixed  length,  where  the  packet 
transmission  time  is  equal  to  a  single  slot. 

The  PSPacket  Arrival  Process 

Packets  (or  messages)  are  assumed  to  arrive  into  the  station  in  accordance 
with  a  Geometric-Batch  arrival  point  process  described  as  follows.  Packet  arrivals 
are  recorded  at  the  end  of  the  arrival  slot. 

A  batch  of  packets  will  arrive  in  a  slot  independently  of  arrival  into  other  slots. 
We  set: 

P{a  batch  of  messages  arrives  in  a  s!ot}=p; 

P{no  messages  arrive  in  a  slot}=1-p; 
where  0<p<1 . 

The  number  of  packets  in  an  arrival  batch  (say  B)  is  a  random  variable,  with 
an  average  that  is  set  to  be  given  by 

The  average  size  of  the  arrival  messages  batch  =  b  [mess./batch]. 

The  ITDMA  program  developed  by  IRI  allows  the  user  to  select  one  of  three 
distributions  for  the  number  of  packets  in  the  batch: 

1.  A  Deterministic  distribution;  under  which  the  number  of  packets  in  the 
arrival  batch  is  always  fixed,  equal  to  b. 

2.  A  Geometric  distribution  with  mean  set  equal  to  b. 

3.  A  uniform  distribution,  whereby  the  number  of  packets  in  a  batch  is 
uniformly  distributed  between  a  specified  lower  value  U,  and  upper  level  U2;  note 
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the  mean  value  is  now  equal  to  (Uj  +U2  )/2. 

Performance  Measures 
CS  Session  Throughput 

The  station’s  CS  session  throughput  (THCS)  is  equal  to  the  number  of 
sessions  per  unit  time  (or  per  slot)  transmitted  by  the  station  across  the  channel. 
Under  the  ITDMA  model  the  station  throughput  rate  is  equal  to  the  fraction  of  session 
arrival  rate  corresponding  to  unblocked  sessions,  so  that  we  have 

THCS  =X(1-PB)  [sess./slot]. 

Note  that  since  each  slot  is  equal  to  t  [sec].,  we  have 

Station  Message  Throughput  =  7v(1-Pb)/t  [sess/sec]. 

The  station  s  normalized  CS  session  throughput,  or  CS  session 
traffic  intensity,  is  equal  to  the  ratio  of  the  traffic  rate  of  arriving  sessions  at  the 
station  and  the  channel’s  session  service  rate  dedicated  to  the  station,  and  is  thus 
given  by 


Normalized  Throughput=pcs=[A(1  -PB)]/(nq/m). 

PS  Packet  Throughput 

The  PS  packet  throughput  (THPS)  is  equal  to  the  number  ofpackets  per 
unit  time  (or  per  slot)  transmitted  by  the  station  across  the  channel.  Under  the  PS- 
TDMA  model  the  station  throughput  rate  is  equal  to  the  message  arrival  rate  (no 
messages  are  blocked)  so  that  we  have 

THPS  =  bp  [mess./slot]. 

Note  that  since  each  slot  is  equal  to  t  [sec.],  where 

T[sec]  =(average  #bits/packet)/(channel  transmission  rate  in  bps), 
we  have,  per  station, 

Station  Packet  Throughput  =  bp/T  [packets/sec] 
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The  station’s  normalized  PS  packet  throughput,  or  packet  traffic 
Intensity,  is  equal  to  the  ratio  of  the  traffic  rate  of  arriving  packets  at  the  station  and 
the  channel’s  PS  packet  service  rate  allocated  to  the  station.  Since  an  average  of  V 
slots  per  framo  are  used  for  CS  support,  the  number  of  slots  per  frame  available  for 
PS  packet  transmissions  is  equal  to  m-n.  Hence  we  have 

Normalized  PS  packetThroughput=pPS=(mbp)/[m-V]. 

For  system  stability  (i.e,  to  ensure  finite  steady  state  system  delays  and  buffer  queue 
sizes  for  PS  traffic)  it  is  necessary  to  ensure  that 


Pps<^ • 


Packet  Delay  and  Queue-Size 

The  performance  of  the  CS  subsystem  is  expressed  by  the  CS  throughput 
presented  above  and  by  the  session  blocking  probability  obtained  by  the  ITDMA 
program.  The  performance  of  the  PS  subsystem  is  expressed  interms  of  the  station 
packet  queue-size  (also  termed  as  system-size)  and  the  packet  delay. 

The  message  waiting-time  is  defined  as: 

W=Message  Waiting  Time=Total  time  elapsed  since  the  message  arrival  slot 
to  the  time  that  the  first  packet  of  the  message  starts  transmission 

The  message  effective  transmission  time  is 

T=message  effective  transmission  time=Time  elapsed  since  the  start  of 
transmission  of  the  session’s  first  packet  to  the  time  that  the  last 
packet  ends  transmission 

The  messge  delay  time  (D)  is  defined  as: 

D=Message  Delay=Total  time  elapsed  since  the  message  arrival  slot  to  the 
time  that  the  last  packet  of  the  message  is  transmitted 


Thus,  we  have 
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D=W+T. 


The  Message  queue-size  is: 

X=Message  queue-size=number  of  messages  resident  in  the  station’s  buffer 
(including  the  message  in  the  process  of  transmission,  if  any) 

Note  that  for  this  ITDMA  program  a  message  consists  of  a  single  packet. 

The  ITDMA  program  provides  results  for  the  mean  and  standard-deviation  of 
the  PS  X  and  D  variables  as  well  as  for  the  distributions 
P(j)=P{X=j}:  D(j)=P{D=j}. 

These  results  are  obtained  by  us  through  the  use  of  two  methods: 

a.  We  develop  a  simulation  program  that  employs  analytic  recurrence 
relationships  expressing  the  evolution  of  the  system  states  for  the  PS  subsystem. 
Random  generators  are  used  to  generate  the  traffic  loading.  The  sample  mean, 
variance  and  distributions  of  the  queue  size  and  delay  variables  are  then  obtained 
through  statistics  collections  and  computations.  The  user  is  able  to  select  the 
number  of  simulation  runs  (slots)  for  the  stop  time  as  well  as  the  start  time  for  the 
collection  of  statistics  to  incorporate  in  the  calculation  of  the  X  and  D  performance 
measures. 

b.  We  use  methods  developed  by  us  and  presented  in  the  References  which 
provide  for  analytical  derivations  of  lower  and  upper  bound  formulas  for  the  steady 
state  mean  queue  sizes  and  message  delays  for  the  PS  subsystem.  The  Cs 
subsystem  is  analyzed  by  using  a  M/Geom/n/n  queueing  analysis. 

Note  that  the  analytical  results  are  steady  state  results  and  thus  describe 
performance  results  after  the  system  has  been  running  for  a  long  time.  Hence,  the 
analytical  based  output  results  can  differ  from  simulation  results  which  are  based  on 
shorter  run  times. 
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4.3  Instructions  for  Running  the  ITDMA  Program 

in  the  following  we  provide  instructions  for  running  the  ITDMA  program 
(Version  1 ).  It  is  noted  that  the  specific  name  of  the  program  used  can  be  PSTDMAx 
where  x  is  a  number  designating  a  version  of  the  program  compiled  to  incorporate 
certain  limits  on  the  program  size.  For  example,  ITDMA2  involves  the  limit  m£50, 
and  printed  values  for  the  queue-size  (x)  and  for  the  delay  (d)  which  are  in  the  range 
x£20,  d£20. 

Step  bv  Step  Instructions  (See  Appendix  E  for  examples) 

1 .  Enter  the  name  of  the  program  (such  as  ITDMA2) 

Subsequent  inputs  are  entered  in  response  to  program  prompts. 

2.  Enter  the  name  of  an  output  file;  no  more  than  12  characters;  in  PC,  quote 
the  ‘file  name’ 

3.  Enter  a  number  to  designate  the  run  method.  Select 

1  -  for  simulation  only 

2  -  for  analysis  only 

3  -  for  simulation  and  analysis 

4.  Enter  the  start  time  for  simulation  collection  start 

5.  Enter  the  stop  time  for  the  simulation  length  duration 

6.  Enter  the  frame  duration  (m,  an  integer,  where  1<m£specified  upper  bound 
such  as  50);  then  also  type  (following  a  space  key  entry)  the  number  of  maximum 
slots  allocatedfor  CS  support  (n,  an  integer,  where  nsm) 

7.  Enter  the  CS  session  arrival  rate  (A>  0  [sess./slot]) 

8.  Enter  the  CS  session  transmission  rate  (q=q1  [frames/session),  where  the 
average  session  holding  time  is  equal  to  1/q  [frames/session]) 

9.  Enter  the  batch  arrival  rate  of  PS  packets  (p=p1 ,  so  that  0<p<1 ,  where  p 
designates  the  probability  that  a  batch  of  packets  arrives  in  a  slot);  each  message  is 
assumed  to  contain  a  single  packet 

10.  The  program  computes  and  exhibits  the  value  of  V,  which  expresees  the 
average  number  of  slots  per  frame  used  by  the  CS  session  traffic;  this  value  must  be 
used  to  ensure  that  the  PS  subsystem  is  stable,  as  noted  in  the  following 

1 1 .  Enter  the  PS  packet  arrival  batch  size  distribution  index;  select 

1  -  for  deterministic  (fixed)  batch  size 

2  -  for  Geometric  batch  size 

3  -  for  uniform  batch  size 

12.  If  selected  deterministic  or  Geometric  batch  distributions,  enter  next  the 
mean  batch  size  b. 

If  selected  a  uniform  distribution,  select  the  lower  and  upper  levels  for 
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the  batch  sizes  Uj  and  u2;  the  program  then  calculates  the  mean  batch  size  as 
b=(u,+u2)/2. 

To  ensure  system  stability,  so  that  thePS  packet  arrival  rate  is  lower 
than  the  channel  average  packet  service  rate,  it  is  necessary  to  ensure  that  the 
normalized  throughput  (rho=pPS)  is  less  than  one,  or  that 

mbp<m-V. 

If  the  parameters  selected  violate  this  condition,  the  program  indicates  so  and 
requests  the  user  to  re-select  the  PS  subsystem  parameters. 

13.  When  simulation  is  used,  enter  the  levels  x  and  d  for  the  program  to 
explicitly  provide  at  the  output  the  probabilities 

P{X>x},  P{D>d) 

for  the  PS  subsystem  queue  size  (X)  and  aelay  (V)  measures. 

14.  The  program  then  responds  with  ‘please  wait’  and  proceeds  with  its  run. 
When  done,  the  output  is  shown  on  the  screen  as  well  as  provided  in  the  specified 
output  file. 


31 


Delay-throughput  Evaluator,  IRI  Corp.,  Dr.  Izhak  Rubin 


4.4  Output  of  the  ITDMA  Program 

Appendix  E  providesan  example  of  the  input  and  output  for  the  ITDMA 
program. 

The  output  of  the  program  contains  the  following  information: 

1 .  Statement  of  the  input  parameters. 

2.  The  normalized  throughput  (rho=pPS)  for  the  PS  subsystem 

3.  When  simulation  used,  a  table  presenting  the  queue-size  and  message 
delay  distributions,  P(j)  and  DQ),  respectively. 

4.  For  simulation:  the  means  and  standard  deviations  for  the  station  PS 
queue  size  and  for  the  PS  packet  delay;  the  probabilities  P(X>x),  P(D>d),  for  the 
selected  x  and  d  values 

5.  The  blocking  probability  for  CS  sessions 

6.  For  analysis:  upper  and  lower  bounds  for  the  means  of  the  station  queue 
size  and  for  the  message  delay 

As  indicated  in  the  References,  the  analytical  method  developed  for  the 
calculation  of  the  upper  and  lower  bounds  is  much  more  computationally  efficient 
than  the  method  developed  for  the  analytical  calculation  of  the  exact  result. 

Further  note  that  the  analytical  steady-state  mean  delay  E(D)  and  mean 
queue  size  E(X)  are  related  through  Little’s  Theorem  so  that 

E(X)=bpE(D). 
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5.  Analytical  Performance  Equations  for  PSTDMA  and  References  for 

the  TPMA  Analyses. 


5.1  References  for  the  TDMA  Analyses 

The  TDMA  analyses  used  in  the  PSTDMA,  CSTDMA  and  ITDMA  programs 
are  based  upon  the  analytical  models,  derivations  and  performance  evaluations 
presented  in  the  following  references  based  upon  research  investigations  carried 
out  by  Professor  Izhak  Rubin  and  his  research  group  at  the  Electrical  Engineering 
department  of  UCLA. 
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5.2  Analytical  Performance  Eve.uatlon  for  PSTDMA:  The  System  and 

Traffic  Model 

Under  the  PSTDMA  program,  the  following  model  has  been  assumed. 

1.  A  station  is  considered.  THe  station  is  allocated  n  (or  N)  Slots  per  frame. 
The  frame  (cycle)  consists  of  m  (or  M)  slots. 

2.  Messages  arrive  at  the  station  at  random  according  to  an  independent 
(Geometric  Batch)  arrival  process,  so  that  if 
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An=number  of  messages  arriving  during  the  n-th  slot 
then  {An.n^l}  is  a  sequence  of  i.i.d.  random  variables.  The  moments  of  A  are  defined 
as: 


a(i)=P{A=i},  teO;  (1) 

with  moments: 

a=E(A)=mean  number  of  message  arrivals  per  slot=A  (2) 

a2=E(A2)  (3) 

a3=E(A3).  (4) 


Equivalently,  the  arrival  process  can  be  represented  as  a  Geometric  batch 
arrival  process,  characterozed  as  follows. 

2.1  A  batch  of  messages  arrives  in  a  slot  with  probability  p,  0<p<1 ;  no 
messages  arrive  in  a  slot  with  probability  1-p. 

2.2  The  message  batches  are  independent  identically  distributed  (i.i.d.) 
random  variables  {Gn,  n£1},  with  distribution 


P(Gn=i)=g,,  tel  (5) 

and  moments 

b=E(G)=Average  batch  size;  k>2=E(G2);  t>3=E(G3)  (6) 

In  relating  the  two  above-mentioned  presentations  of  the  independent  arrival 

process,  the  following  relationships  hold: 

g,«a(i)/[1-a(0)].k1;  (7) 

p=1-a(0);  (8) 

a=A=pb  (9) 

a2=pb2  (10) 

a3=pb3.  (11) 


In  the  following,  we  present  examples  of  batch  distributions. 
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Geometrically  Distributed  Batches 

For  a  Geometric  batch  size  distribution  with  mean  batch  size  set  equal  to  b,  we 
have: 


g1,=(1-b-'r1(b-,),ial; 

(12) 

E(G)=b 

(13) 

b2=E(G2)=(2-b_,)b2 

(14) 

b3=E(G3)=b3(6-6b-,+b-2) 

(15) 

Deterministic  (Fixed!  Batch  Sizes 

For  batch  sizes  .  .  _,i  are  assumed  to  be  of  fixed  (deterministic)  size  equal  to  b, 
we  have 


b2=t>3=b 


(16) 


Uniformly  Distributed  Batch  Sizes 

For  a  batch  size  which  is  assumed  continuously  uniformly  distributed  in  (u1,u2), 
we  have: 


b=(u1+u2)/2 

(17) 

b2=(u23-u,3)/l3(u2-u,)] 

(18) 

b3=(u24-ui4)/l4(u2-u1)]. 

(19) 

Poisson  Arrival  Process 

Under  a  Poisson  arrival  process  with  intensity  A  {mess/slot],  we  have 


a(i)=P(A=i)=exp(-/\)A7i!  ,  iaO; 

(20) 

a=E(A)=A 

(21) 
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a2=A+A2  (22) 

a3=A+3A2+A3.  (23) 

3.  A  message  contains  a  random  number  of  packets.  The  packet  transmission 
time  is  equal  to  a  single  slot  duration. 

Let  Bn  denote  the  number  of  packets  contained  in  the  n-th  messsage.  The 
packet  length  sequence  {Bn,n>1}  consists  of  i.i.d.  random  variables,  characterized  by 
the  distribution  and  moments: 

p(i)=P(B=i),  tel ;  (24) 

(3-fJ|=E(B);  |32=E(B2);  {33=E(B3).  (25) 


The  following  special  message  length  distributions  are  noted. 

Geometrically  Distributed  Message  Length 

For  a  Geometrically  distributed  message  length  with  parametr  q  and  mean 
message  length  1/q,  we  have: 


p(i)«P(B=i)=(1-qr'q,  tel;  (26) 

p-|3pE(B)=1/q;  (27) 

02*E(B2)=(2-q)/q2;  (28) 

P3»E(B3)=(6-6q+q2)/q3.  (29) 

Fixed  (Deterministic)  Message  Length 

If  the  message  contains  a  fixed  number  of  packets,  set  equal  to  (3,  we  have: 

B=p;  (3,=(3;  (32=P2;  P3=|33-  (30) 


Uniformly  Distributed  Message  Length 

For  amessage  length  which  is  assumed  continuously  uniformly  distributed  in 
(U),u2),  we  have: 


p=(U]+u2)/2 

(31) 

p2=(u23-u,3)/l3(u2-u,)] 

(32) 

p3=(u24-u14)/I4(u2-u1)]. 

(33) 
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The  Frame  Arrival  Process 

The  number  of  arrivals  during  the  n-th  frame  is  set  equal  to  Nn,n>1.  The  frame 
arrival  sequence  {Nn,n>1}  consists  of  i.i.d.  random  variables,  with  distribution  and 
moments  given  by: 


n(i)=P(N=i),  teO;  (34) 

n,=E(N);  n2=E(N2);  n3=E(N3).  (35) 

When  the  slot  arrival  process  {An,n>1}  is  assumed  to  be  an  independent 
process  (to  consist  of  i.i.d.  random  variables),  as  assumed  here,  we  have: 

N„=T/V  (36) 

where  M  denotes  the  number  of  slots  per  frame.  Hence, 

n,=E(N)=Ma=MA;  (37) 

n2=E(N2)=(M2-M)a2+Ma2;  (38) 

n3=E(N3)=M(M-1  )(M-2)a3+3M(M-1  )aa2+Ma3.  (39) 


5.3  Analytical  Performance  Evaluation  for  PSTPMA:  ContiauousSlot 
Asglflnmg.ni,Pgr  .station 

The  PSTDMA  program  assumes  contiguous  assignment  of  station  slots  in  each 
frame.  The  station  is  thus  allocated  N  (or  n)  contigouguous  slots  during  each  frame. 
THe  frame  contains  M  (or  m)  slots. 

Under  such  a  slot  allocation,  our  analysis  yields  the  following  performance 
results.  Assumed  here  is  a  Geometric  batch  arrival  process  with  batch  arrival 
intensity  p,  0<p<1,  batch-size  mean  b  and  moments  and  and  message  length 
which  is  Geometrically  distributed  with  parameter  q  and  mean  1/q;  see  Section  5.2 
for  details  of  the  traffic  load  definitions  and  relationships. 

The  steady-state  mean  queue-size  expressing  the  average  number  of 
messages  resident  in  the  station  buffer,  E(X),  is  obtained  to  be  given  by  the  following 
equations. 
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E(X)=E(X1)+E(X2)+E(X3),  (1) 

where 

E(X,)-[2pb(1  -pb)+pb22]/[2(q-pb)];  (2) 

E(X2)={q2[b22-pb2+pb(2-q)])/I2(q-pb)]H(M-N)/(Nq-Mpb));  (3) 

N-t 

E(X3)=((M-N)q2/(2M)KX;  lAtWjJQfWjJ+wJIAlWjJQlWjJ-Wjl  (4) 

J-l 

and 

A(z)=1-p+pG(z)  (5) 

Q(z)=q+(1-q)z  (6) 

where  B(z)  is  the  generating  function  of  the  batch  size  distribution 

00 

G(z)-y  iP(G-i)  (7) 


and 

b=E(G),b22=E{G(G-1)}=b2-b.  (8) 

The  numbers  {Wj,j=1 ,2,...N}  are  the  N-1  roots  of 

zn*A(z)mQ(z)n  inside  |z|<1.  (9) 

The  mean  message  delay  is  then  obtained  by  Little's  formula  to  be: 

E(D)=E(X)/pb.  (10) 

The  station’s  normalized  traffic  intensity  is 

p=  (Mpb)/(Nq).  (11) 
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The  conditiori  for  system  stability,  so  that  a  finite  message  delay  is  attained,  is 

p<1,  or  Mpb<Nq.  (12) 

Upper  and  Lower  Bounds  for  the  Mean  Queue=Size  and  Message  Delay 

We  have  developed  equations  for  upper  and  lower  bounds  for  the  mean 
queue-size  and  mean  message  delays,  which  are  easy  to  compute.  These 
approximations  are  given  as  follows.  They  are  used  in  the  PSTDMA  program.  We 
have: 


Q<E(X)£Q+[(N-1  )/2][(M-N)q/M], 


(13) 


where 

Q=[2a(1-a)+a22]/I2(q-a)]+{Lq[a22-a2+a(2-q)M2(q-a)(Nq-Ma)].  (14) 

and 


a=E(A)=pb;  a22=E(A2)-E(A)=a2-a=p(b2-b).  (1 5) 

We  note  that  the  upper  and  lower  bounds  for  the  mean  queue-size  differ  by  at  most  a 
constant  value  of  (N-1  )/2.  For  N=1 ,  a  station  is  allocated  a  single  slot  per  cycle  and 
the  upper  and  lower  bounds  coincide,  yielding  E(X)=G. 

The  mean  message  delay  upper  and  lower  bounds  are  then  computed  using 
the  relationship: 


E(D)=E(X)/pb=E(X)/a.  (16) 


5.4  Analytical  Performance  Evaluation  for  PSTDMA;  Single  and 
Uniform  Slot  Per  Frame  Assignments 

In  this  section  we  make  the  assumption  that  the  station  is  assigned  a  single  slot 
(N=1)  during  a  frame  containing  M  slots. 

If  the  station  is  assigned  slots  within  a  frame  in  a  uniformly  distributed  manner, 
so  that  the  station  is  allocated  N  slots  in  a  cycle  of  duration  L  slots  and 


L/N=M 


(1) 
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where  M  is  an  integer,  the  same  channel  configuration  conditions  occur  so  that  the 
station  is  in  fact  assigned  1  slot  every  M  slots,  and  the  results  of  this  section  can  be 
applied. 

Additionally,  the  results  of  this  section  can  also  be  used  as  am  approximation  if 
the  ration  L/N  is  not  an  integer,  whereby  we  use  M  to  denote 

M=average  number  of  frame  slots  per  each  allocated  station  slot.  (2) 

The  key  advantage  of  assuming  a  uniform  distribution  of  station  slots  per  frame 
(or,  equivalently  a  single  slot  per  cycle  allocation)  is  that  it  induces  a  much  simpler 
analytical  model  that  we  employ  to  compute  exact  results  for  the  mean  and  the 
variance  of  the  message  delay,  E(D)  and  Var(D)=cr2(D),  resp.  The  ability  to 
analytically  compute  the  standard-deviation  <x(D)  is  of  prime  importance  in  many 
applications  in  which  the  tail  probability  of  the  delay  must  be  computed  and 
constrained.  Thus,  for  example,  we  can  use  the  analytically  computed  value 

D.99  ~  E(D)+3cr  (D),  (3) 

to  provide  an  estimate  of  the  99-th  percentile  delay,  for  which  we  guarantee 

P(D^D  93)2:0.99  (4) 

so  that  the  message  delay  D  can  be  guaranteed  to  be  lower  than  D  99  for  99%  of  the 
served  messages.  (Note  that  for  a  Gaussian  distribution  estimate  3  yields  a  99.9% 
probability,  while  for  an  exponential  distribution,  for  which  cr(D)=E(D),  this  estimate 
represents  a  98%  probability.) 

In  deriving  the  results  for  the  PSTDMA  system,  we  assume: 

1.  An  independence  arrival  provess  with  parameters  as  described  in  Section 
5.2.  Thus,  the  arrival  process  parameters  involve  the  moments  a=A,  a2  and  a3. 
Equivalently,  the  Geometric  batch  arrival  model  description  induces  the  batch  arrival 
probability  p  and  the  batch  moment  parameters  b,  b2  and  b3.  We  also  use  the  frame 
arrival  moments  n,  r^  and  n3.  (See  Section  5.2.) 

2.  The  message  can  contain  a  random  number  of  packets.  The  message 
length  distribution  has  the  moments  P,  P3-  (See  Section  5.2.) 

3.  The  station  is  assigned  a  single  slot  every  M  slots,  as  described  above.  The 
packet  transmission  time  is  equal  to  the  duration  of  a  single  slot. 


40 


Delay-throughput  Evaluator,  IRI  Corp.,  Dr.  Izhak  Rubin 


The  mean  message  delay  E(D)  is  cjven  by  the  following  formulas. 

E(D)=ME(Wd)+[(M-1  )/2]+(M(3+1  -M),  (5) 

where 

E(WD)-(p/[2(1-p)]){(f>2/f>)+[p(n2-n,)/ni2l-1)+(p/2)t(n2/n,)-1],  (6) 

and  the  normalized  throughput  (traffic  intensity)  is 

p=n1(3=MAp<1.  (7) 

For  example,  for  the  special  case  of  a  Poisson  arrival  process  with  intensity 
A[mess/slot],  we  obtain 


E(D)={Mp/[2(1  -p)]}(f32/p)+M|3+1  -M/2.  (8) 

In  deriving  Eq.  (8),  we  have  added  1/2  slot  to  the  delay  to  account  for  the  continuous¬ 
time  random  nature  of  the  Poisson  arrival  process,  noting  that  a  message  may  arrive 
anywhere  within  the  slot,  but  according  to  the  the  discrete  model  its  arrival  is 
recorded  at  the  end  of  the  slot. 

The  mean  delay  expression  can  also  been  written  as: 

E(D)={Mp/[2<1  -p)M(M3)-1  )+{p2/I2A(1  -P)])(MA-A-1  +a  2I\) 

+(Mf3/2)[(M- 1  )A+a2/A]+(M|3 +1  -M)/2.  (9) 


The  variance  of  the  message  delay,  Var(D),  is  given  by  the  following  formulas. 


Var(D)=Var(W,)+Var(W2)+Var(F)+Var(T). 

(10) 

where 

VarfFMMp-IJ/IZ 

(11) 

Var(T)=M2(p2-p2). 

(12) 
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We  have: 


where 


Var(W,)=M2{Wt2+Wn  -  Wn2}, 

(13) 

Var(W2)=M2{W22+W21  -  W2]2}, 

(14) 

W„  =  {n2p  12-n1(Pi2+p1-p2)}/[2(1-p)]; 

(15) 

W21  =  (n^n^p^n,. 

(16) 

Note  that  E(WD)=Wn+W2i  .  For  the  computation  of  W,2  and  W22,  we  present  the 
following  formulas. 

W]2  ={M6(1'p)2]},{3{n2-n])p]2+3n1((32-p1)+2(1-p)[(n3-3n2+2n1)p13 

+3(n2-n,)pt(p2-p|)+n((|33-3p2+2^()j};  (1 7) 

W22  ={1/I6pPi2]}*{3Pi(p2-Pi)[(n2*n|)p,2+n|(p2-P])]+2pPj(P3-3p2+2p,) 
-2p)2[(n3-3n2+2ni)p]3+3(n2-ni)p1(p2-p1)+r,)(p3-3p2+2pi)] 

'  3p(02-f3i)2}.  (18) 
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6.  Delav-Throuahout  Performance  Evaluator  for  Timed  Token  Protocol 
(FPPI-I  Type)  Token-Ring  Systems 

6.1  Introduction 

In  this  Chapter,  we  describe  the  structure  and  provide  the  operation  directions 
for  the  FDDI  program,  which  provides  a  simulation  based  delay-throughput  analysis 
of  token  ring  networks  employing  the  timed  token  protocol  (FDDI-I  type)  medium 
access  control  (MAC)  scheme,  as  used  by  SAFENET  II..  In  Section  6.2,  the 
structures  of  the  service,  channel,  message  and  traffic  models  are  described,  and  the 
performance  measures  are  defined.  In  Section  6.3,  we  provide  detailed  instructions 
for  running  the  FDDI  program.  Examples  and  the  source  codes  are  given  in  the 
appendices.  The  output  is  discussed  in  Section  6.4. 

6.2  Hie  Structure  of  the  FDDI  Medium  Access  Control  Model 

The  Network  Topology 

The  network  is  configured  as  a  token  ring,  where  each  station  is  actively 
inserted  into  the  ring  through  its  ring  interface  unit.  Each  pair  of  neighboring  stations 
is  connected  by  a  unidirectional  fiberoptic  link.  (  A  contra  revolving  ring  is  also  used 
to  increase  the  network  reliability;  however,  its  existence  does  not  affect  the  delay- 
throughput  performance  of  the  network,  which  is  based  on  the  regular  use  of  a  single 
link.) 

The  network  topology  is  illustrated  in  Fig.  4. 

We  set: 

N=Number  of  Stations. 

The  access  protocol  is  based  on  a  token-based  distributed  polling  procedure. 
A  single  token  (packet  of  a  determined  pattern)  exists  across  the  ring.  A  station  is 
allowed  to  transmit  some  of  its  packets  (in  accordance  with  the  timed  token  protocol 
described  below)  only  when  it  acquires  the  token.  Upon  the  termination  of 
transmission  of  its  messages,  the  station  immediately  (early  token  release) 
generates  a  new  token  and  transmits  it  across  the  unidirectional  link  to  its 
neighboring  station. 

A  critical  parameter  determining  the  delay-throughput  efficiency  of  the  token 
ring  system  is  the  underlying  walk  time.  The  ring  walk  time  represents  the  time  it 
takes  for  the  token  to  walk  around  the  ring  when  no  station  wishes  to  transmit.  It 
involves,  as  components,  the  time  it  takes  for  the  station  to  generate  the  token, 
transmit  it  to  the  neighboring  station  and  the  time  it  takes  for  the  latter  station  to  fully 
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receive  the  token,  recognize  it  and  act  properly  (i.e.,  immediately  release  it  if  the 
station  does  not  wish  to  transmit  across  the  channel,  or  remove  the  token  from  the 
ring  if  the  station  wishes  to  transmit  its  own  messages  across  the  ring). 

We  set: 

r(i)  =  walk  time  for  the  token  from  station  i  to  station  i+1 . 

Key  components  included  in  the  walk  time  r(i)  are: 

-  the  token  transmission  time  (which  is  equal  to  the  token  length  divided  by 
the  transmission  rate  across  the  channel) 

-  the  token  propagation  delay  from  station  i  to  station  i+1 ,  which  depends 
upon  the  length  of  the  fiber  between  these  two  stations,  and  is  computed  by  using 
the  medium  propagation  rate  of 

Propagation  rate  =  5p.sec/km. 

Example:  For  example,  for  a  system  which  uses  a  token  which  is  100  bits  long,  a 
channel  transmission  rate  of  10  Mbps  and  an  average  interstation  fiber  length  of  5 
km,  we  have: 

Token  transmission  time  =  100/100M  =  1p.sec; 

Interstation  propagation  time  *  5)jsec/km  x  5km  *  25  jasec; 
so  that  (when  other  delay  are  neglected)  we  have: 

r(i)  =  26  ^isec. 


S4 
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The  Message  Structure 

Messages  arrive  at  random  at  the  station’s  buffer  (whose  capacity  is  not 
limited)  and  are  served  on  a  packet-switched  store-and-forward  priority  basis,  in 
accordance  with  the  timed  token  protocol. 

Messages  arriving  at  a  station  are  categorized  into  the  following  priority 
classes: 

-  Synchronous  messages,  which  are  granted  the  highest  priority 

-  Asynchronous  messages  which  are  further  classified  into  multiple  priority 
classes,  so  that  priority-j  asynchronous  messages  are  granted  higher  priority  than 
priority-(j+1 )  asynchronous  messages.  Thus,  priority  -1  messages  are  provided 
higher  priority  for  accessing  the  ring  than  priority-2  messages. 

We  set 


Np(i)  =  number  of  message  priority  classes  for  station  i. 

In  our  current  (Version  1)  FDDI  program,  each  message  is  assumed  to  consist 
of  a  single  packet.  The  packet  is  assumed  to  require  a  transmission  time  across  the 
link  which  is: 

Packet  transmission  time  =  Uniformly  distributed,  with  a  specified  mean. 

In  the  FDDI  program,  we  set: 

plens  (i)  ^average  transmission  time  for  synchronous  packet  at  station  i  [msec] 

plena(i.j)  =  average  transmission  time  for  an  asynchronous  packet  at  station  i, 
considering  a  packet  of  priority  class  j. 


The  Message  Arrival  Process 

Messages  are  assumed  to  arrive  into  the  station  in  accordance  with  a  Poisson 
process,  with  the  arrival  rates  given  as  follows: 

as(i)  ■  arrival  rate  of  synchronous  messages  (packets)  at  station  i 
[packets/msec] 

aa(i,j)  *  arrival  rate  of  priority-j  asynchronous  packets  at  station  i 
[packets/msec] 
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The  Timed  Token  Rotation  access  protocol  and  its  parameters 


The  Timed  Token  Rotation  access  protocol  operates  as  described  in  the 
following.  We  also  present  here  the  associated  system  and  program  parameters  that 
must  be  selected. 

During  initialization,  a  target  time  is  selected  for  token  rotation,  which  is  set  in 
our  FDDI  program  to  be: 

Target  Token  Rotation  Time  =  TTRT  [msec]. 


Upon  the  receipt  of  a  token,  the  station  can  transmit  its  synchronous  packets, 
provided  the  total  time  it  spends  for  this  transmission  of  synchronous  packets  during 
this  visit  of  the  token  is  not  longer  than  a  specified  time  limit  called  the  bandwidth 
lima.  We  thus  set  in  the  FDDI  program: 


Bandwidth  time  =  BWT  [msec].  % 


b 

O 


0 

a 


To  ensure  that  if  all  stations  transmit  synchronous  packets  for  the  maximum 
duration  BWT  the  token  rotation  time  will  not  exceed  TTRT,  we  must  select  BWT  so 
that 


NxBWT+  ]jT  r(i)  +  Max[asynchronous  packet  length]  £TTRT  7 

/ 

where  £  r(i)  represents  the  total  walk  time  around  the  ring  (which  is  equal  to  Nr, 
when  r(i)=r,  each  i,  for  symmetrically  located  stations  around  the  rinc^. 

Each  station  keeps  a  record  of  the  token  rotation  time  (TRT)  representing  the 
total  time  elapsed  from  the  instant  that  it  was  visited  by  the  token  for  the  last  time  to 
the  instant  of  current  visit  of  the  token. 

Following  the  receipt  of  a  token,  assume  the  station  to  have  transmitted 
synchronous  packets  for  a  total  transmission  time  of  TT(0).  Assume  the  station  holds 
in  its  buffer  priority-1  asynchronous  packets  which  it  wishes  to  transmit,  to  determine 
if  the  latter  can  be  transmitted  the  following  time  threshold  is  selected: 

T_pri(l )  ■  time  threshold  for  priority-1  packet  transmissions  [msec]. 
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The  station  then  compares  the  time  duration  TRT+TT(0)  with  this  threshold.  If 
TRT+TT(0)  <  T_pri(1 ), 

the  station  can  transmit  a  priority-1  packet.  Another  priority  1  packet  can  be 
transmitted  provided 

TRT+TT{0)+TT(1 )  <  T_pri(1 ), 

where  TT(1)  represents  the  total  transmission  time  (up  to  this  time)  of  priority-1 
packets  during  this  token  visit. 

Upon  the  termination  of  transmission  of  priority-1  packets,  the  station  is 
allowed  to  transmit  its  priority-2  packets  (if  it  has  any)  provided 

TRT+TT(0)+TT(1  )+TT(2)  <  T_pri(2), 

where  TT(2)  represents  the  total  transmission  time  of  priority-2  packets  heretofore 
during  this  token  visit,  and 

T_pri(2)  =  time  threshold  for  priority-2  packet  transmissions  [msec]. 

Similarly,  for  the  priority  based  access  control  for  asynchronous  packets  of 
lower  priority  levels. 

Our  FDDI  program  provides  for  a  controlled  exhaustive  service  at  a  station,  so 
that  synchronous  packet  transmissions  can  take  place  after  the  termination  of  an 
asynchronous  packet  transmission,  provided  their  underlying  transmission  threshold 
has  not  been  crossed  yet.  Similarly,  for  asynchronous  packet  transmissions. 


Performance  Measures 
IhfflUQhPUl 

The  station’s  throughout  (THS)  is  equal  to  the  number  of  messages  per  unit 
time  transmitted  by  the  station  across  the  channel.  The  station  i  priority-j  throughput 
rate  is  then  equal  to 


THs(i.j)  =aa(i,j)  [packets/msec]. 
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while  the  normalized  throughput  for  station  i  and  priority-j  packets  is  given  by 

Normalized  Throughput  (i,j)=p(ilj)=aa(i,j)plena(i,j)  [Erlangs]. 

For  synchronous  traffic  (which  is  also  denoted  here  as  priority-0  traffic  class) 
we  have: 

THs(i.O)  =as(i)  [packets/msec]. 

Normalized  Throughput  (i,0)=p(i,0)=as(i)plens(i)  [Erlangs]. 

The  network  normalized  throughput  for  priority-j  traffic  is  equal  to 


pG)-2>('.i).  £0. 

I 

For  system  stability  (i.e,  to  ensure  finite  steady  state  system  delays  and  buffer 
queue  sizes)  it  is  necessary  to  ensure  that 

P=I>(i)<1.  £0. 

j 

Under  the  current  version  of  the  FDDI  program,  the  station  buffer  is  limited  (to 
200  packets),  leading  to  a  blocking  probability  PB(i,j)  for  station-i  and  priority-j 
packets.  Unless  the  traffic  intensity  levels  are  very  high,  the  blocking  probabilities  are 
kept  very  low.  In  turn,  when  very  high  traffic  intensities  (such  as  p>0.95  )  are 
assumed,  buffer  overflows  can  occur  and  distinctive  blocking  probabilities  result.  In 
this  case,  the  above  formulas  represent  the  offered  traffic  loads.  THe  carried  traffic 
loads  are  computed  by  multiplying  the  offered  traffic  loads  by  (1-Pb).  for  the  proper 
blocking  probability  factor  (PB). 

Message  Delay  and  Queue-Size 

The  message  waiting-time  is  defined  as: 

W=Packet  Waiting  Time=Total  time  elapsed  since  the  packet  arrival 

Instant  to  the  instant  that  the  packet  starts  transmission 

The  packet  transmission  time  is: 
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T=packet  transmission  time 

The  packet  delay  time  (D)  is  defined  as: 

D=Packet  Delay=Total  time  elapsed  since  the  message  arrival  time  to  the 
time  that  the  packet  is  transmitted 


Thus,  we  have 


D=W+T. 


ThePacket  queue-size  is: 

X=Packet  queue-size=number  of  packets  resident  in  the  station’s  buffer  at  the 
time  of  a  token  arrival  to  the  station 

The  FDDI  program  measures  also  other  queue-sizes: 

Xa=Packet  queue-size  at  packet  arrival  times  to  the  station 

X^Packet  queue-size  at  packet  departure  times  from  the  station 

Xt=Packet  queue-size  at  a  station  at  an  arbitrary  time 

The  FDDI  program  provides  results  for  the  mean  and  standard-deviation  of  X, 
W  and  D  as  well  as  for  the  distributions 

P(j)=P{X=j};  W(w)=P{W£w},  D(d)=P{D£d}. 

In  addition,  the  FDDI  program  provides  statistical  information  about  other 
system  variables: 

Vj=dwell  time  for  priority-j  messages  at  a  station 
■  total  transmission  time  of  priority-j  packets  during  a  station’s  dwell  time 
(token  visit  time) 

V«dwell  time  of  a  station  =  total  time  used  by  a  station  to  transmit  packets 
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across  the  ring  during  a  single  token  visit 

C  =  cycle  duration  (relative  to  a  station)  =  time  elapsed  between  successive 
visits  of  the  token  to  a  station 

In  addition  to  obtaining  performance  results  for  each  station  and  for  each 
message  priority  class,  we  compute  global  performance  measures  that  represent  the 
means  and  standard  deviations  of  X,  W  and  D  over  all  network  stations,  for  each 
priority  class. 

To  obtain  these  results,  we  have  developed  a  simulation  program  that 
employs  analytic  recurrence  relationships  expressing  the  evolution  of  the  system 
states.  Random  generators  are  used  to  generate  the  traffic  loading.  The  sample 
mean,  variance  and  distributions  of  the  queue  size,  packet  delay  and  other  system 
variables  are  then  obtained  through  statistics  collections  and  computations.  The 
user  is  able  to  select  the  simulation  stop  time  as  well  as  the  start  time  for  the 
collection  of  statistics  to  incorporate  in  the  calculation  of  the  performance  measures. 
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6.3  Instructions  for  Running  the  FDDI  Program 

In  the  following  we  provide  instructions  for  running  the  FDDI  program  (Version 
1).  It  is  noted  that  the  specific  name  of  the  program  used  can  be  FDDIx  where  x  is  a 
number  designating  a  version  of  the  program  compiled  to  incorporate  certain  limits 
on  the  program  size. 

Step  bv  Step  Instructions  (See  Appendices  G-l  for  examples) 

1 .  Enter  the  name  of  the  program  (such  as  FDDI) 

Subsequent  inputs  can  be  entered  by  a  data  file  or  in  response  to  program 
prompts.  Appendices  G-H  describe  both. 

2.  Enter  the  name  of  an  output  file 

3.  Enter  a  number  to  designate  the  network  feature.  Select 

1  -  symmetric  system 

2  -a  system  with  2  classes  of  stations 

3  -a  system  for  which  the  traffic  loading  can  change  from  station  to 

station 

Under  a  symmetric  system  feature,  the  stations  are  assumed  to  be  statistically 
identical,  so  that  the  same  traffic  loads  are  offered  to  each  station.  All  interstation 
walk  times  are  assumed  identical.  Appendix  G  illustrates  the  input  and  output 
formats  of  this  case. 

Under  system  feature  2,  the  network  consists  of  two  classes  of  stations. 
Stations  belonging  to  the  same  class  are  assumed  statistically  identical.  This  model 
allows  the  evaluation  of  networks  where  class  1  stations  are  characterized  with 
certain  service  types  and  traffic  loads  while  class  2  stations  experience  different 
service  mixtures  and/or  higher  traffic  loads,  or  serve  as  gateways  to  other  networks 
and  thus  involve  different  traffic  characteristics.  Appendix  H  illustrates  the  input  and 
output  formats  of  this  case. 

Under  system  feature  3,  an  asymmetric  network  configuration  can  be  defined, 
where  each  station  can  have  its  own  uniquely  defined  traffic  loading  features. 
Appendix  I  illustrates  the  input  and  output  formats  of  this  case. 

In  the  following,  we  illustrate  the  inputs  for  a  symmetric  system;  they  are 
similar  for  the  other  cases. 

4.  Enter  the  start  time  for  simulation  collection  start  [msec] 

5.  Enter  the  stop  time  for  the  simulation  length  duration  [msec] 

6.  Enter  w,d,x  for  computing  P(W>w),  P(D>d),  P(X>x). 

7.  Enter  the  number  of  stations  (N) 

8.  Enter  the  walk  time  (r  [msec]) 
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9.  Enter  the  value  selected  for  TTRT  [msec] 

10.  For  synchrnous  traffic,  one  then  enters: 

10.1  packet  arrival  rate  (as,  packets/msec  per  station) 

10.2  mean  packet  transmission  time  (plens,  m^c) 

10.3  Bandwidth  time  (BWT,  msec) 

1 1 .  For  asynchrnous  traffic,  for  each  priority  class  j  (and  for  each  station  type 
for  feature  cases  other  than  symmetric)  one  then  enters: 

11.1  packet  arrival  rate  (aa(i.j),  packets/msec  per  station) 

1 1 .2  mean  packet  transmission  time  (plena(i.j),  msec) 

11.3  Threshold  time  (T_pri(j),  msec) 
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6.4  Output  of  the  FDDI  Program 

Appendices  G-l  provide  three  examples  of  the  inputs  and  outputs  for  the  FDDI 
program. 

The  output  of  the  program  contains  the  following  information: 

1.  Statement  of  the  input  parameters. 

2.  The  normalized  throughput  load  offered  and  realized.  Also  computed  is 
the  theoretical  maximum  normalized  throughput. 

3.  Total  number  of  cycles  realized  during  the  simulation  run  is  presented. 

4.  We  then  present  the  means  and  standard  deviations  for  the  station  queue- 
size  (X),  for  message  wait-time  (W)  and  for  message  delay  (D).  These  values  are 
first  averages  over  all  the  stations,  for  each  message  priority  class  (and  for  each 
station  class  for  configuration  feature  2).  These  values  are  then  presented  for  each 
station. 

5.  For  each  station  and  each  message  priority  class,  performance  results  are 
then  exhibited  for: 

5.1  mean  and  standard  deviation  of  vj,  the  station  dwell  time  for  the 
transmission  of  priority  j  messages 

5.2  the  probabilities  P(X>x),  P(W>w),  P(D>d),  for  the  selected  values  of 

x,  w  and  d 

6.  For  each  station,  presented  are  results  for 

6.1  mean  and  standard  deviation  of  V,  the  station’s  dwell  time 

6.2  mean  and  standard  deviation  of  C,  the  station’s  cycle  time 

6.3  mean  and  standard  deviation  of  C-V,  the  station’s  vacation  time 

7.  For  each  station  and  for  each  message  priority  class  served  by  each 
station,  performance  results  are  presented 

7.1  for  the  queue  size  distributions:  P(X=j),  P(Xa=j).  P(Xd=j),  P(Xt=j) 

7.2  for  the  message  wait-time  distribution  P(W^w) 

7.3  for  the  message  delay-time  distribution  P(D^d) 

8.  Also  noted  is  the  realized  normalized  throughput  per  station. 
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Appendix  A. 


The  PS-TDMA  Program:  Input  and  Output 
Examples 
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pstdma  2 

Enter  the  f ile_name  for  data: 

(no  more  than  12  characters) 

(in  PC,  QUOTE  the  'file_name') 

' pout  1  ' 

What  would  you  Tike  to  do  ? 

1:  Simulation  only 
2 :  Analysis  only 
3:  Simulation  and  analysis 
Please  enter  a  number: 

3 

Enter  the  start  time  (when  statistic  starts) 

10 

Enter  the  stop  time  (total  simulation  time) 

1000 

Enter  the  frame  duration  (m)  (<=500)  and  number 
of  slots  (n)  allocated  to  a  station  (n<=m) 
m  and  n  should  be  integer 
1 0  5 

Enter  the  batch  arrival  rate  [batch. /slot]  pi 
0  .  1 

Enter  the  transmission  rate  [mess. /slot]  ql 
0 . 5 

Enter  the  batch  size  distribution  index 
1:  deterministic 
2:  geometric 
3:  uniform 


2 

Enter  the  mean  batch  size,  b  ( b*p1 <n*q 1/m) : 

2 

To  find  the  probabilities  P(X>x)  and  P(D>d),  please 
enter  x  and  d  (x<=100,  d<=100) 
x  and  d  should  be  integer 
10  10 
please  wait 
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**  Simulation  and  Analysis  ** 

PS-TDMA  System  Parameters: 
m ,  n  =  10  5 

the  mean  batch  arrival  rate  pi  1 . OOOOOOOOOOOOOOOE-OO 1 

the  mean  service  rate  ql  5 . OOOOOOOOOOOOOOOE-OO 1 

system  throughput  rho  8 . OOOOOOOOOOOOOOOE-OO 1 

the  batch  size  is  ;  geometric 

the  mean  batch  size  is  :  2.000000000000000 

The  queue-size  probabilities  are: 


j 

P(j)  D(j) 

0 

.200101010 

.000000000 

1 

.091616162 

.329718739 

2 

.072525253 

.  0  1  5249068 

3 

. 083030303 

.  0  1  1  860386 

4 

. 062828283 

. 0 17282277 

5 

.071919192 

.016265673 

6 

.068989899 

.019654354 

7 

. 052020202 

.020332091 

8 

. 037676768 

.018976618 

9 

. 042929293 

.013893595 

1  0 

. 027272727 

.019315486 

1  1 

.025252525 

.014910200 

12 

. 020404040 

.013554727 

13 

.017777778 

.014571332 

14 

.015858586 

.010504914 

1  5 

.015757576 

.01  1521518 

1  6 

.012828283 

.014232464 

1  7 

. 013737374 

.013893595 

1  8 

.  007272727 

.014571332 

1  9 

.  006363636 

.015926804 

20 

.002323232 

.010504914 

21 

. 003232323 

.012876991 

22 

.  003232323 

.012538123 

23 

. 002424242 

.012538123 

24 

.001313131 

.  0  1  321  5859 

25 

.  002929293 

.010843782 

26 

. 002525253 

.012876991 

27 

.  002222222 

.010504914 

28 

.003737374 

.012876991 

29 

.001010101 

.010504914 

30 

.  002424242 

.008471705 

31 

.001010101 

.012199254 

32 

.001212121 

.012538123 

33 

.002525253 

.  005083023 

34 

.004040404 

.010166045 

35 

.001515152 

.  006777364 

36 

.002121212 

.006438495 

.002121212 

.  005421  891 

38 

.001818182 

.007116232 

57 


39 

. 003434343 

.005083023 

40 

.001212121 

. 005421891 

4  1 

.002424242 

. 006099627 

42 

.001919192 

.  0074551  00 

43 

.000808081 

.005421891 

44 

.  000202020 

.007455100 

45 

.000101010 

.006099627 

46 

. 000000000 

. 004405286 

47 

. 000000000 

. 003727550 

48 

.000000000 

.0047441 55 

49 

.  000000000 

.005421891 

50 

. 000000000 

. 006438495 

51 

. 000000000 

.006099627 

52 

. 000000000 

. 006438495 

53 

. 000000000 

.003727550 

54 

. 000000000 

.  003049814 

55 

. 000000000 

.002710945 

56 

. 000000000 

. 0040664 1 8 

57 

. 000000000 

.004066418 

58 

.  000000000 

. 0071 16232 

59 

. 000000000 

. 0040664 1 8 

60 

.  000000000 

. 002372077 

61 

.  000000000 

. 002372077 

62 

.000000000 

.  003388682 

63 

. 000000000 

.  003049814 

64 

.  000000000 

.002033209 

65 

.  000000000 

. 002372077 

66 

.000000000 

. 002033209 

67 

. 000000000 

.002033209 

68 

.  000000000 

.001355473 

69 

.000000000 

. 003388682 

70 

.  000000000 

.  002033209 

71 

.000000000 

.003388682 

72 

.  000000000 

.000677736 

73 

.  000000000 

. 000677736 

74 

.  000000000 

.  002372077 

75 

.000000000 

.001016605 

76 

.  000000000 

.  00  1  355473 

77 

.  000000000 

.001694341 

78 

.  000000000 

. 000000000 

79 

.  000000000 

.  000677736 

80 

. 000000000 

.001694341 

81 

.  000000000 

.000338868 

82 

.  000000000 

.001355473 

83 

.  000000000 

. 000677736 

84 

.000000000 

.  000338868 

85 

.  000000000 

.  000677736 

86 

.  000000000 

.  000677736 

87 

.000000000 

.001016605 

88 

. 000000000 

.  000338868 

89 

.  000000000 

.001016605 

90 

. 000000000 

. 000000000 

58 


91 

. 000000000 

.000338868 

92 

.000000000 

.000338868 

93 

.000000000 

.000338868 

94 

. 000000000 

. 000338868 

95 

.000000000 

.000338868 

96 

. 000000000 

.000677736 

97 

. 000000000 

. 000677736 

98 

. 000000000 

.001016605 

99 

. 000000000 

.001355473 

=  100 

. 000000000 

.039308709 

P(X>  10)=  .189090909 

P(D>  10)=  .517451711 

mean  queue-size  (simulation)  is  :  6.394646464646464 

standard  deviation  of  queue-size  (simu.)  is:  7.603460993452924 

mean  delay  (simulation)  is:  22.430362588952900 

standard  deviation  of  delay  (simulation)  is:  28.780133446363980 


mean  queue-size  (analysis)  is  :  3.376208812762674 

mean  queue-size  (analysis,  lower  bound)  is:  3.350000000000001 

mean  queue-size  (analysis .upper  bound)  is:  4.350000000000001 

mean  delay  (analysis)  is:  1 6 . 1 C 4 4 C538  1  33 7 0 

mean  delay  (analysis,  lower  bound)  is:  16.750000000000000 

mean  queue-size  ( ana  1 ys i s , upper  bound)  is:  21.750000000000010 
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C>  DStdma2 

Enter  the  file_name  for  data: 

(no  more  than  12  characters) 

(in  PC,  QUOTE  the  ' f i1e_name' ) 

' pout2 ' 

What  would  you  like  to  do  ? 

1:  Simulation  only 
2:  Analysis  only 
3:  Simulation  and  analysis 
Please  enter  a  number: 

3 

Enter  the  start  time  (when  statistic  starts) 

10 

Enter  the  stop  time  (total  simulation  time) 

1000 

Enter  the  frame  duration  (m)  (<=500)  and  number 
of  slots  (n)  allocated  to  a  station  (n<=m) 
m  and  n  should  be  integer 
10  1 

Enter  the  batch  arrival  rate  [batch . /s 1 ot ]  pi 

0.01 

Enter  the  transmission  rate  [mess. /slot]  ql 

0.2 

Enter  the  batch  size  distribution  index 
1:  deterministic 
2:  geometric 
3 :  uni  form 


3 

Enter  the  range  of  the  batch  size: 

(e.g.,  “1.0,  3.0")  ( u 1 +u 2 ) *p 1 / 2 <n*q 1 /m 

1  .  0 . 2 . 6 

To  find  the  probabilities  P(X>x)  and  P(D>d),  please 
enter  x  and  d  (x<=100,  d<=100) 
x  and  d  should  be  integer 
20  10 
please  wait 
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**  Simulation  and  Analysis  ** 

PS-TDMA  System  Parameters: 
m ,  n  =  10  1 

the  mean  batch  arrival  rate  pi  1 . 000000000000000E-002 

the  mean  service  rate  ql  2 . OOOOOOOOOOOOOOOE-OO 1 

system  throughput  rho  9 . OOOOOOOOOOOOOOOE-OO 1 

the  batch  size  is  :  uniform 

the  range  of  the  batch  size  is:  [  1.000000000000000 

2 . 600000000000000] 


The  queue-size  probabilities  are: 


j 

P(J) 

D(j) 

0 

. 049898990 

.000000000 

1 

.033838384 

.202479339 

2 

.091414141 

.000000000 

3 

. 064848485 

. 000000000 

4 

.076969697 

.000000000 

5 

. 060909091 

.004132231 

6 

.096464646 

. 000000000 

7 

. 084242424 

.000000000 

8 

. 059191919 

. 000000000 

9 

.  1  08888889 

.000000000 

1  0 

.  068484848 

.004132231 

1  1 

.065959596 

. 000000000 

12 

.074949495 

. 000000000 

1  3 

.  028888889 

.000000000 

14 

.019494949 

.000000000 

1  5 

.  007474747 

. 000000000 

16 

.006060606 

.000000000 

17 

. 002020202 

. 000000000 

1  8 

.000000000 

.000000000 

1  9 

.  000000000 

. 000000000 

20 

.  000000000 

. 000000000 

2  1 

.  000000000 

.000000000 

22 

. 000000000 

.000000000 

23 

.000000000 

.000000000 

24 

.000000000 

.000000000 

25 

. 000000000 

.000000000 

26 

.000000000 

. 000000000 

27 

.  000000000 

. 000000000 

28 

.000000000 

.000000000 

29 

.  000000000 

. 000000000 

30 

.  000000000 

. 000000000 

31 

.  000000000 

.  000000000 

32 

.  000000000 

. 000000000 

33 

.  000000000 

. 000000000 

34 

.000000000 

. 000000000 
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35 

. 0000G0000 

.008264463 

36 

.000000000 

.000000000 

37 

.000000000 

.000000000 

38 

.000000000 

.000000000 

39 

. 000000000 

.000000000 

40 

. 000000000 

. 000000000 

41 

. 000000000 

. 000000000 

42 

. 000000000 

.000000000 

43 

.000000000 

. 000000000 

44 

. 000000000 

. 000000000 

45 

. 000000000 

. 000000000 

46 

. 000000000 

.000000000 

47 

.000000000 

.004 132231 

48 

. 000000000 

. 000000000 

49 

.000000000 

.000000000 

50 

. 000000000 

.008264463 

51 

.  000000000 

.000000000 

52 

. 000000000 

. 000000000 

53 

.000000000 

.000000000 

54 

.000000000 

. 000000000 

55 

. 000000000 

.000000000 

56 

.000000000 

.000000000 

57 

.000000000 

. 000000000 

58 

.  000000000 

.000000000 

59 

.  000000000 

. 000000000 

60 

.  000000000 

.004132231 

61 

. 000000000 

. 000000000 

62 

.  000000000 

.004132231 

63 

.000000000 

. ooooooooo 

64 

.000000000 

.000000000 

65 

.  000000000 

. ooooooooo 

66 

.000000000 

. ooooooooo 

67 

.  000000000 

.ooooooooo 

68 

. 000000000 

.ooooooooo 

69 

.  000000000 

.ooooooooo 

70 

.  000000000 

. ooooooooo 

71 

.  000000000 

.  ooooooooo 

72 

.  000000000 

.004 132231 

73 

. 000000000 

-OOOOOOOOO 

74 

.  000000000 

. ooooooooo 

75 

.  000000000 

.ooooooooo 

76 

.  000000000 

.  004  1  3223  1 

77 

.000000000 

.004 1 32231 

78 

.  000000000 

.  ooooooooo 

79 

.  000000000 

.ooooooooo 

80 

.000000000 

.  ooooooooo 

81 

.  000000000 

. ooooooooo 

82 

. 000000000 

.ooooooooo 

83 

.000000000 

.ooooooooo 

84 

. 000000000 

.  004  1  32231 

85 

.  000000000 

. ooooooooo 

86 

.  000000000 

. ooooooooo 

87 

.  000000000 

.  ooooooooo 
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88 

. 000000000 

.000000000 

89 

. 000000000 

. 000000000 

90 

.000000000 

.  000000000 

91 

. 000000000 

.  000000000 

92 

. 000000000 

. 000000000 

93 

. 000000000 

. 000000000 

94 

.000000000 

.000000000 

95 

■  .000000000 

.000000000 

96 

. 000000000 

.000000000 

97 

.000000000 

. 000000000 

98 

. 000000000 

. 000000000 

99 

. 000000000 

. 000000000 

=  100 

. coooooooo 

.743801653 

P(X>  20)=  .000000000 

P(D>  10)=  .789256198 

mean  queue-size  (simulation)  is  :  6.847373737373737 

standard  deviation  of  queue-size  (simu.)  is:  3.899230508336784 

mean  delay  (simulation)  is:  267.053719008264500 

standard  deviation  of  delay  (simulation)  is:  202.557489423080300 


mean  queue-size  (analysis)  is  :  2.565600000000002 

mean  queue-size  (analysis,  lower  bound)  is:  2.565600000000002 

mean  queue-size  (ana 1 ys i s , upper  bound)  is:  2.565600000000002 

mean  delay  (analysis)  is:  142.533333333333400 

mean  delay  (analysis,  lower  bound)  is:  142.533333333333400 

mean  queue-size  ( ana 1 ys i s , upper  bound)  is:  142.533333333333400 
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Appendix  B.  The  PS-TDMA  Program:  Source  Code 
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C: 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C; 

c : 

C: 

C: 

C: 

C: 


Packet-Switched  TDMA  System 

*  Version  1  Sept.  20,  1989  * 

*  Professor  Izhak  Rubin  * 

*  Department  of  Electrical  Engineering  * 

*  Univeristy  of  California  * 

*  Los  Angels,  CA  90024  * 


Descriptions  of  the  variables  (partial) : 
pi  :  mean  message  batch  arrival  rate 
b  :  mean  number  of  message  arrivals  per  slot 
ql  :  mean  message  service  rate 
types  :  the  type  of  batch  size  distribution 
u  :  deterministic  message  btach  size 

ul  :  lower  limit  of  the  uniformly  dist.  message  batch  size 

u2  :  uuper  limit  of  the  uniformly  dist.  message  batch  size 

seedt  :  the  seed  to  generate  the  number  of  batch  arrival 

seeds  :  the  seed  to  generate  the  number  of  messages  served 

ic(j)  :  total  #  of  messages  in  the  system  at  the  start  of  jth  slot 
(j-0,1,2, . . .,100) 

p(j)  :  the  prob.  that  there  are  j  messages  in  the  system  at  the 
start  of  the  jth  slot 

qmean  :  mean  queue-size  calculated  from  simulation 
wmean  :  mean  message  delay  calculated  from  little's  formula 
External  function:  rand (seed) 

return  a  value  that  is  uniformly  distributed  over  (0,  1) 


implicit  real*8 (a-h,o-z) 
complex*16  x(500) 
real* 8  mu 

integer  xn (500) , seedt, seeds, types, 

+  batch, bs, icd (0: 100,1: 500 ), id (0:100) , 

+  ict (0 : 100, 1 : 500) , ic (0 :100) , typesl, seedt 1, dn (500) 
common  mu, u, ul, u2, seeds 
character*12  filen 
mm-lOO 

print  *,  'Enter  the  file_name  for  data:  ' 
print  *,  '(no  more  than  12  characters)' 
print  *,  '  (in  PC,  QUOTE  the  "  f  ile_name"  ) ' 
read  *,  filen 
open (10, f ile-filen) 

print  *,  'What  would  you  like  to  do  ?  ' 
print  *,  '  1:  Simulation  only  ' 
print  *,  '  2:  Analysis  only  ' 
print  *,  '  3:  Simulation  and  analysis  ' 

print  *,  'Please  enter  a  number:' 
read  *,  typesl 
if (typesl .eq. 1)  then 
write(10,*)  ' 

write (10,*)  '**  Simulation  **' 

print  *,  'Enter  the  start  time  (when  statistic  starts) 
read  *,  istart 

print  *,  'Enter  the  stop  time  (total  simulation  time)' 
read  *,  istop 
elseif (typesl .eq. 3)  then 
write(10,*)  ' 

write  (10,*)  '**  Simulation  and  Analysis  **' 

print  *,  'Enter  the  start  time  (when  statistic  starts) 

read  *,  istart 

print  *,  'Enter  the  stop  time  (total  simulation  time)' 
read  *,  istop 
elseif  (typesl .eq. 2)  then 
write(10,*)  ' 


write (10, *) 


Analysis  ****' 
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9  +*** 


c 

88 


endif 


print  *,  'Enter  the  frame  duration  (m)  (<«500)  and  number' 
print  *,  '  of  slots  (n)  allocated  to  a  station  (n<-m) ' 
print  *,  '  m  and  n  should  be  integer  ' 
read  *,  m,n 

Pri°t  *,  'Enter  the  batch  arrival  rate  [batch. /slot]  pi  ' 
read  *,  pi 

print  *,  'Enter  the  transmission  rate  [mess. /slot]  ql  ' 
read  *,  ql 

1  print  *,  'Enter  the  batch  size  distribution  index' 
print  *,  '  1:  deterministic' 

print  *,  '  2:  geometric' 

print  *,  '  3:  uniform' 

print  *,  '  ' 

read  * ,  types 
if (types. eq.l)  then 

P^lnt  *t  'Enter  the  fixed  batch  size  u  (u*pl<n*ql/m) : ' 
read  *,  u 
rho«u*pl*m/ (n*ql) 
if  (rho.ge.l.dO)  then 
print  *,  '  ' 

print  *,  'system  not  stable,  please  reenter  the  parameters  ' 
go  to  88 
endif 
b~u*pl 

a2-pl*u* (u-1) 
elseif (types.eq.2)  then 

print  *,  'Enter  the  mean  batch  size,  b  (b*pl<n*ql/m)  : ' 
read  *,  b 
rho-b*pl*m/ <n*ql) 
if  (rho.ge.l.dO)  then 
print  '  ' 

print  *,  'system  not  stable,  please  reenter  the  parameters  ' 
go  to  88 
endif 

mu-1 .OdO/b 
b-b*pl 

a2-pl*2 . dO* (1-mu) /mu**2 
elseif (types. eq. 3)  then 

print  *,  'Enter  the  range  of  the  batch  size:' 
print  *,  '  (e.g.,  ”1.0,  3.0")  (ul+u2) *pl/2<n*ql/m  ' 

read  *,  ul,u2 
rho- (ul+u2) *pl*m/ (2*n*ql) 
if  (rho.ge.l.dO)  then 
print  '  ' 

print  *,  'system  not  stable,  please  reenter  the  parameters  ' 
go  to  88 
endif 

b-pl* (ul+u2) /2 .dO 

a2-pl* (u2* (u2+l) * (2*u2+l) /6-u2# (u2+l) /2- (ul-1) *ul* (2*ul-l) 
*  / 6+ul* (ul-1) /2) / (u2-ul+l) 

else 

go  to  1 

endif 

write  (10,*)  '  PS-TDMA  System  Parameters:  ' 
write(10,*)  'm,  n  -  ',m,n 

write  (10,*)  'the  mean  batch  arrival  rate  pi  ',pl 
write  (10,*)  'the  mean  service  rate  ql  ',ql 
write  (10,*)  'system  throughput  rho  ' , rho 
if (types. eq.l)  then 

write (10,*)  'the  batch  size  is  :  deterministic  ' 
write  (10,*)  'the  fixed  batch  size  is  :  ',u 
elseif (types .eq. 2)  then 

write  (10,*)  'the  batch  size  is  ;  geometric  ' 
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write (10,*)  'the  mean  batch  size  is  :  ',1/mu 
elseif (types. eq. 3)  then 

write  (10,*)  'the  batch  size  is  :  uniform  ' 

write  (10,*)  'the  range  of  the  batch  size  is:  [' ,ul, ' : ' ,u2, ' ] ' 
endif 
c 

if (typesl .eq. 2)  goto  125 

print  *,  'To  find  the  probabilities  P(X>x)  and  P (D>d) ,  please  ' 
print  *,  'enter  x  and  d  (x<-100,  d<-100) ' 
print  *,  'x  and  d  should  be  integer' 
read(5,*)  igx,igd 
print  *,  'please  wait  ' 
seedt-1099 
seedtl-1099 
seeds-99999 
do  14  i-0,irm 
do  15  j-l,m 
ict (i, j)-0 
icd(i, j)-0 

15  continue 
ic (i) -0 
id(i)-0 

14  continue 
c 

totalx-0 . OdO 

tolx2-0 . OdO 

ktd-0 

tdd-0 

tdd2-0 

do  16  i-l,m 

16  xn(i)-0 
c 

do  117  i-1,500 
117  dn(i)«0 

npl-n+1 

c 

c  **  simulation  starts  here 
c 

do  20  i-l,istop 
mn-nb ( seedt , pi ) 
if (xn (m) .eq. 0)  then 
if(mn.eq.O)  then 
xn(l)-0 
else 

bs-batch (types) 
xn (1) -mn*bs 
endif 
else 

ndn-nb ( seedt 1 , ql ) 
if(mn.eq.O)  then 
xn ( 1 ) -xn (m) -ndn 
else 

bs-batch (types) 
xn (1) -xn (m) -ndn+bs*mn 
endif 
endif 

11-xn (1) +ndn 
minx-min (500,11) 
do  149  kd-l,minx 
149  dn (kd)-dn (kd) +1 

c 

if(ndn.eq.l)  then 
kdd-dn ( 1 ) 
minxl-minx-1 
do  119  kd-l,minxl 
119  dn (kd)-dn (kd+1) 


68 


dn (minxl+1) -0 
endif 
c 

if  (i.gt.istart)  then 
if (xn (1) .gt .mm)  then 
ict (mm, 1) -ict (mm, 1) +1 
else 

ict (xn (1) , l)-ict (xn  (1) , 1) +1 
endif 

totalx-totalx+xn (1) 
tolx2-tolx2+xn (1) **2 
if(ndn.eq.l)  then 
ktd-ktd+1 
tdd-tdd+kdd 
tdd2«tdd2+kdd**2 
if (kdd . gt .mm)  then 
icd(mm, l)-icd(mm, 1)+1 
else 

icd(kdd, 1) -icd (kdd, 1) +1 
endif 
endif 
endif 
c 

do  17  j«2,n 
mn-nb(seedt,pl) 
if (xn ( j-l) . eq. 0)  then 
if(mn.eq.O)  then 
xn ( j) “0 
else 

bs-batch (types) 
xn ( j) -mn*bs 
endif 
else 

ndn-nb ( seedt 1 , ql ) 
if (mn.eq.O)  then 
xn ( j) -xn ( j-l) -ndn 
else 

bs-batch (types) 
xn(j)«  xn( j-l)-ndn+bs*mn 
endif 
enaif 

11-xn ( j) +ndn 
minx-min (500,11) 
do  139  kd-l,minx 
139  dn (kd) -dn (kd) +1 

c 

if(ndn.eq.l)  then 
kdd-dn ( 1 ) 
minxl-minx-1 
do  129  kd-l,minxl 
129  dn (kd) -dn (kd+1 ) 

dn (minxl+1) -0 
endif 
c 

if (i .gt . istart )  then 
if (xn ( j) .gt .mm)  then 
ict (mm, j ) -ict (mm, j ) +1 
else 

ict (xn ( j) , j ) -ict (xn ( j) , j) +1 
endif 

totalx-totalx+xn ( j ) 
tolx2-tolx2+xn ( j) **2 
if(ndn.eq.l)  then 
ktd-ktd+1 
tdd-tdd+kdd 
tdd2-tdd2+kdd**2 
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if (kdd.gt .mm)  then 
icd (mm, j ) -icd (mm, J ) +1 
else 

icd(kdd, j)-icd(kdd, j)+l 
endif 
endif 
endif 

17  continue 

do  18  j«npl,m 
mn-nb  (seedt,pl) 
bs-batch (types) 
xn ( j ) -xn ( j— 1 ) +mn*bs 
c 

minx-min ( 500 , xn ( j ) ) 
do  211  kd-l,Tninx 
211  dn(kd)-dn(kd)+l 

c 

if (i .gt . istart )  then 
if (xn ( j ) . gt .mm)  then 
ict (mm, j / -ict (mm, j ) +1 
else 

ict (xn ( j ) , j ) *ict ( xn ( j ) , j ) +1 
endif 

totalx-totalx+xn ( j ) 
tolx2«tolx2+xn ( j ) **2 
endif 

18  continue 
c 

20  continue 

do  22  j«0,mm 
do  21  k“l , m 
ic ( j)-ic ( j)+ict ( j,k) 
id( j)-id( j)+icd( j, k) 

21  continue 

22  continue 
c 

pd-0 

px-0 

last-istop-istart 

print  *,'The  queue-size  probabilities  are  as  follows:' 


write  (6,*)  'j  P(j)  D(j)  ' 

print  *,' - • 

write  (10,*)  'The  queue-size  probabilities  are:' 
write(10,*)  'j  P  ( j )  D(j>  ' 

write(10,*)  ' - • 


qmean-totalx/ (last*m) 

var-dsqrt ( (tolx2-m*last*qmean**2) / (last*m-l) ) 
dmean«tdd/ (ktd*l . OdO) 

vard-dsqrt(  (tdd2-ktd*dmean**2) / (ktd-1 . OdO)  ) 
do  110  i-0 , mm 
if(i.lt .mm)  then 

write (6, 200)  i, ic (i) *1 .d0/ (m*last) ,id(i) *1.0d0/ktd 
write (10, 200)  i, ic (i) *1 .d0/ <m*last ) , id (i) *1 . OdO/ktd 
else 

write (6, 210)  i,ic(i) *l.d0/ (m*last) , id (i) *1. OdO/ktd 
write  (10, 210)  i, ic (i) *1 . d0/ ( last *m) , id (i) *1 . OdO/ktd 
endif 

if (i.gt .igx)  px-px+ic (i) *1 .d0/ (m*last) 
if(i.gt.igd)  pd-pd+id (i) *1 . dO/ktd 
110  continue 

2CC  format (lx, i3, 2x, 2f 14 . 9) 

210  format (lx, '>-' , i3,2x, 2f 14 .9) 

write (6,*)  'mean  queue-size  (simulation)  is  :  ' , qmean 

write(6,*)  'standard  deviation  of  queue-size  (simu.)  is:',var 
write(6,*)  'mean  delay  (simulation)  is:  ',dmean 
write (6,*)  'standard  deviation  of  delay  (simulation)  is:  ',vard 
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write (6, 300)  igx,px 
write (6, 301)  igd, pd 
write (10, 300)  igx,px 
write  (10, 301)  igd,pd 

300  format (lx, 'P (X>' , i3, ' )-  ',fl4.9) 

301  format (lx, 'P(D>',i3,')-  ',fl4.9) 

write  (10,*)  'mean  queue-size  (simulation)  is  :  ',qmean 
write(10,*)  'standard  deviation  of  queue-size  (simu.)  is:  ',var 
write (10,*)  'mean  delay  (simulation)  is:  ',dmean 
write  (10,*)  'standard  deviation  of  delay  (simulation)  is:',vard 
if (typesl .eq. 1)  then 
goto  126 
else 

goto  127 
endif 

125  print  *,  'please  wait' 

127  call  tdma (x, n,m, pi, ql, suml, types) 

Q  ★★★★★★★ 

d-ql* (2*b* (1-b) +a2) / (2* (ql-b) ) +ql**2* (a2-b**2+b* (2-ql) ) * (m-n) 

*  / (2* (ql-b) * (n*ql-m*b) ) 
amean«d+ (m-n) *ql**2/ (2*m) *suml 

Q  ★  ★  ★  ★  ★ 

write  (6,*)  'mean  queue-size  (analysis)  is  :  ',amean 
write  (6,*)  'mean  queue-size  (analysis,  lower  bound)  is:  ',d 
dl-d+ (m-n) * (n-1) *1.0d0/2/m 

write  (6,*)  'mean  queue-size  (analysis, upper  bound)  is:  ',dl 
write  (6,*)  'mean  delay  (analysis)  is:  ',amean/b 
write  (6,*)  'mean  delay  (analysis,  lower  bound)  is:  ',d/b 
write  (6,*)  'mean  queue-size  (analysis, upper  bound)  is:  ',dl/b 
write  (10,*)  'mean  queue-size  (analysis)  is  :  ',amean 
write  (10,*)  'mean  queue-size  (analysis,  lower  bound)  is:  ',d 
write  (10,*)  'mean  queue-size  (analysis, upper  bound)  is:  ',dl 
write  (10,*)  'mean  delay  (analysis)  is:  ',amean/b 
write  (10,*)  'mean  delay  (analysis,  lower  bound)  is:  ',d/b 
write  (10,*)  'mean  queue-size  (analysis, upper  bound)  is:  ',dl/b 

126  close  (10) 
stop 

end 

c 

subroutine  tdma (x, n,m, pi, ql, suml, types) 

implicit  complex*16  (a-h,o-z) 

real*8  pi, q, ql, suml,mu,u,ul, u2 

dimension  x(l) 

integer  seeds, types 

common  mu, u, ul , u2 , seeds 

X0-.2 

1-m-n 

call  root (xO, 1, n, pi, ql, x, types) 
nl-n-1 
sum-0 . OdO 
do  15  i-l,nl 

q-ql+(l-ql)*x(i) 

if  (types.eq.l)  then 

a-1 . 0d0-pl+pl*x (i) **int (u*l .001) 

elseif (types .eq. 2)  then 

a«l-pl+pl*mu*x  (i)  /  (1-  (1-mu)  *x  (i) ) 

else 

a-l-pH-pl*  ( x  ( i )  **int  (ul*l .  01)  -x  (i)  **int  (u2+l  .01)  )  /  (l-x(i) ) 

*  /(u2-ul+l) 
endif 

sum-sum+ (a*q+x (i) ) / (a*q-x (i) ) 

15  continue 
c 

suml-real (sum) 

return 

end 
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c 

c 

subroutine  root (xO,l,n,p,q, x, types) 

implicit  complex*16  (a-h,o-z) 

real*8  p, q,mu, u, ul, u2 

integer  seeds, types 

dimension  x(l) 

common  mu, u,ul,u2, seeds 

nl-n-1 

do  15  k=l,nl 
xl-xO 

call  fix (k,l,n,p,q,xl,y, types) 
x(k)«y 
15  continue 
return 
end 

subroutine  f ix (k, 1, n,p, q, xO, y, types) 

implicit  complex*16  (a-h,o-z) 

real*8  p,delt, q,mu,u,ul,u2 

integer  types, seeds 

common  mu, u, ul, u2, seeds 

itmax-1000 

delt-1.0d-12 

i-1 

10  y-g(xO,k,l,n,p,q, types) 
if (i .eq. itmax)  goto  20 
if (cdabs (y-xO) . le . delt )  goto  15 
xO-y 
i-i+1 
go  to  10 

20  write (6,*)  'maximum  number  reached' 

stop 

write  (6, *)  '  x  -  ' , y 
write (6, *)  'f(x)  -  ',x0-y 
15  return 

end 
c 
c 

function  g (x, k, 1, n,p, q, types) 
implicit  complex*16  (a-h,o-z) 
common  mu, u,ul,u2, seeds 
real*8  p,pl,q,ql,t, tl,mu,u,ul,u2 
integer  seeds, types 
ql-q+(l-q) *x 
if  ( types. eq.l)  then 
a-1 . 0d0-p+p*x**int (u*l . 001) 
elseif (types .eq. 2)  then 
a-l-p+p*mu*x/ (1- (1-mu) *x) 
else 

a»l-p+p* (x**int (ul*l . 01) -x**int (u2+l . 01) ) / (1-x) 
*  / (u2-ul+l) 

endif 

pl-2 . 0d0*3 . 14159d0*k/n 
t-dcos (pi) 
tl— dsin  (pi) 

g-dcmplx (t,tl)*ql*a**(l. 0d0+l*l . OdO/n) 

return 

end 

c 


integer  function  batch (types) 
integer  seeds, types 
real*8  mu,u,ul,u2,rand 
common  mu, u,ul,u2, seeds 
if (types. eq.l)  then 


72 


ooooo  oooo 


batch-int (u*l .001) 
elseif (types. eq. 2)  then 

batch=igeom(l .  OdO/mu,  seeds) 
else 

batch-ul+int (  (u2-ul+l) *rand (seeds) ) 
endif 

c  write (6,*)  'batch  -  func  re  batch 
return 
end 

- Random  Generator  Function - 

This  function  will  take  the  argument  "seed"  as  an  "input"  seed, 

and  return  the  value  which  is  uniformly  distributed  over  [0,  1] , 
and  the  output  seed  to  be  used  as  the  next  input  seed, 
function  rand (seed) 
integer  seed 

double  precision  pp,i7, i2,rand 
i7-7**5 
i2-2**31-l 
pp-i7*seed 
seed-dmod (pp, i2 ) 
rand-seed/ (i2+l) 
return 
end 


***********************  function  igeom  **************************** 

FUNCTION  IGEOM  PRODUCES  A  GEOMETRICALLY  DISTRIBUTED  PSEUDO-RANDOM 
OBSERVATION  WITH  AVERAGE  AVG  ,  FROM  R.N.  STREAM  ISTRM 
FUNCTION  IGEOM (AVG, ISTRM) 

REAL* 8  THETA, GI, FI, avg,uu, rand 
Uu-  RAND (istrm) 

THETA  -  1 . OdO/AVG 
NN— INT ( 1 0 . 0  d0  * AVG ) 

GI-THETA 
FI-THETA 
1-1 

IF (Uu.LE.FI)  GOTO  20 
DO  10  1-2, NN 
GI  -  (1 .-THETA) *GI 
FI  -  FI  +  GI 
IF (Uu.LE.FI)  GOTO  20 
10  CONTINUE 
C  WRITE (*,600)  TNOW, AVG, Uu 

c600  FORMAT (IX, ' ! ! ! !  ERROR  IN  %  IGEOM  %,  TIME  -  ',F12.4, 
c  */2X, '  AVG  — ' , F12 . 4 , '  Uu  -  ',F8.5) 

20  IGEOM  -  I 
RETURN 
END 


c  This  function  will  take  the  value  1  with  proba.  p  and  0  with  prob.  1-p 
function  nb(seed,p) 
real*8  t,p,rand 
integer  seed 
t-rand(seed) 
if(t.le.p)  then 
nb-1 
else 
nb-0 
endif 
return 
end 
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Appendix  C. 


The  CS-TDMA  Program:  Input  and  Output 
Examples 
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c> 

C>cstdma2 

Enter  the  f i1e_name  for  data: 

(no  more  than  12  characters) 

(in  PC,  QUOTE  the  'file_name') 

' cout  1  ' 

What  would  you  like  to  do  ? 

1:  Simulation  only 
2:  Analysis  only 
3:  Simulation  and  analysis 
Please  enter  a  number: 

3 

Enter  the  start  time  (when  statistic  starts) 

10 

Enter  the  stop  time  (total  simulation  time) 

1000 

Enter  the  frame  duration  (m)  (<=500)  and  number 
of  slots  allocated  to  a  station  (n) 
n<=m,  m  and  n  should  be  integers 
10  5 

Enter  the  batch  arrival  rate  pi 

0.01 

Enter  the  transmission  rate  al 
0 . 1 

Enter  the  botch  size  distribution  index 
1 :  deterr inistic 
2:  geometric 
3:  unifo  m 


2 

Enter  the  mean  batch  size,  b  ( b*p1 <n*q 1 /m) : 

4 

To  find  the  probabilities  P ( X>x )  and  P(D>d),  please 
enter  x  anr  d  (  x,  d<=100.  both  should  be  integers) 
10  10 
please  wait 
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CS-TDMA  EXAMPLE  1 


**  Simulation  and  Analysis  ** 

CS-TDMA  System  Parameters: 
m ,  n=  10  5 

the  mean  batch  arrival  rate  pi  1 . 000000Q00QQ0000E-002 

the  mean  service  rate  ql  1 . OOOOOOOOOOOOOOOE-OOI 

system  throughput  rho  =  7 . 9 9 9 9 9 9 9 9 9 99 9 9 9 9 E- 00 1 

the  batch  size  is  :  geometric 

the  mean  batch  size  is  :  4.000000000000000 

The  queue-size  probabilities  are: 


j 

P(j) 

D(j) 

0 

. 062424242 

.000000000 

1 

. 067777778 

.357392317 

2 

.068787879 

. 1 12922002 

3 

. 085353535 

.011 641444 

4 

. 066262626 

. 002328289 

5 

. 047575758 

. 000000000 

6 

. 058686869 

. 000000000 

7 

. 030606061 

.0221 18743 

8 

. 029696970 

. 005820722 

9 

. 028484848 

.002328289 

10 

.034343434 

.001164144 

1  1 

. 030000000 

.001 164144 

12 

.014949495 

. 002328289 

1  3 

.028787879 

. 002328289 

14 

.017070707 

.  000000000 

1  5 

.027676768 

.000000000 

16 

.021414141 

.000000000 

17 

.  043333333 

.001164144 

18 

. 029898990 

.002328289 

1  9 

.  026464646 

.004656577 

20 

.  026767677 

.000000000 

21 

.01 1010101 

.001 164144 

22 

.012121212 

.002328289 

23 

.  026969697 

.  000000000 

24 

.010808081 

.  000000000 

25 

.  008282828 

.001 164144 

26 

.  003131313 

.  000000000 

27 

.005151515 

.001164144 

28 

.015454545 

.005820722 

29 

.  023232323 

.000000000 

30 

.010303030 

.  002328289 

31 

.  0  13838384 

.001164144 

32 

.  004040404 

.001 164144 

33 

.  005757576 

.  000000000 

J4 

.  000909091 

.001164144 

35 

.  000202020 

.000000000 

77 


36 

. 002424242 

. 000000000 

37 

.000000000 

. 000000000 

38 

.000000000 

. 005820722 

39 

. 000000000 

. 002328289 

40 

.000000000 

.001 164144 

4  1 

.000000000 

.001164144 

42 

.000000000 

.001164144 

43 

.000000000 

. 000000000 

44 

. 000000000 

.000000000 

45 

. 000000000 

.001 164144 

46 

. 000000000 

. 002328289 

47 

. 000000000 

.001 164144 

48 

.000000000 

.002328289 

49 

. 000000000 

. 002328289 

50 

.000000000 

.000000000 

51 

. 000000000 

. 000000000 

52 

. 000000000 

.001 164144 

53 

. 000000000 

.000000000 

54 

.000000000 

. 000000000 

55 

. 000000000 

. OOOOGuOOO 

56 

. 000000000 

.  002328289 

57 

.  000000000 

. 002328289 

58 

.  000000000 

.001164144 

59 

. 000000000 

. 002328289 

60 

.  000000000 

.002328289 

61 

. 000000000 

.  003492433 

62 

.  000000000 

.  000000000 

63 

.  000000000 

.  002328289 

64 

.  000000000 

.001164144 

65 

.000000000 

.0011 64144 

66 

.000000000 

.001164144 

67 

.  000000000 

.002328289 

68 

.000000000 

.001164144 

69 

.  000000000 

.  003492433 

70 

. 000000000 

.001 164144 

71 

.000000000 

.003492433 

72 

.  000000000 

.000000000 

73 

.  000000000 

.  000000000 

74 

. 000000000 

.002328289 

75 

.  000000000 

.  000000000 

76 

.  000000000 

.001164144 

77 

.  000000000 

.  000000000 

78 

.  000000000 

.001 164144 

79 

. 000000000 

.002328289 

80 

.  000000000 

. 000000000 

81 

.  000000000 

.0011 64144 

82 

.  000000000 

.001164144 

83 

.  000000000 

.001164144 

84 

.  000000000 

. 000000000 

85 

.  000000000 

.  002328289 

86 

.  000000000 

. 000000000 

87 

.  000000000 

.001 164144 

88 

.  000000000 

.000000000 
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89 

. 000000000 

.001 164144 

90 

.000000000 

. 000000000 

91 

. 000000000 

.001 164144 

92 

. 000000000 

.001164144 

93 

.000000000 

. 002328289 

94 

. 000000000 

.001 164144 

95 

.000000000 

.000000000 

96 

.000000000 

.001 164144 

97 

. 000000000 

. 003492433 

98 

. 000000000 

.008149010 

99 

.000000000 

.001 164144 

=  100 

.000000000 

.364377183 

mean  queue-size  (simulation)  is  :  10.655757575757580 

standart  deviation  (queue-size,  simu.)  is:  8.976742237234827 


mean  delay  (simulation)  is:  113.440046565774200 

standard  deviation  (delay,  simulation)  is:  195.965372555913300 


P(X> 

P(D> 

II  II 

/~\  / — s 

o  o 
« —  <— 

.42000000000 

.48428405122 

mean 

mean 

queue- 

queue- 

size  (analysis,  lower  bound) 
size  (analysis,  upper  bound) 

i  s  : 
i  s  : 

1  5 . 740000000000000 
1 9 .700000000000000 

mean 

mean 

del  ay 
delay 

(analysis,  lower  bound)  is: 
(analysis,  upper  bound)  is: 

393.499999999999900 

492.499999999999900 
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cstdma2 

Enter  the  file_name  for  data: 

(no  more  than  12  characters) 

(in  PC,  QUOTE  the  'file_name') 

' cout  2 ' 

What  would  you  like  to  do  ? 

1:  Simulation  only 
2:  Analysis  only 
3:  Simulation  and  analysis 
Please  enter  a  number: 

3 

Enter  the  start  time  (when  statistic  starts) 

10 

Enter  the  stop  time  (total  simulation  time) 

1000 

Enter  the  frame  duration  (m)  (<=500)  and  number 
of  slots  allocated  to  a  station  (n) 
n<=m,  m  and  n  should  be  integers 
1 0  1 

Enter  the  batch  arrival  rate  pi 

0.01 

Enter  the  transmission  rate  ql 
0 . 2 

Enter  the  batch  size  distribution  index 
1:  deterministic 
2:  geometric 
3:  uniform 


2 

Enter  the  mean  batch  size,  b  ( b*p 1 <n*q 1 /m ) : 

1  .  8 

To  find  the  probabilities  P(X>x)  and  P(D>d),  please 
enter  x  and  d  (  x,  d<=100,  both  should  be  integers) 
10  10 
please  wait 
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CS-TDMA  EXAMPLE  2 


**  Simulation  and  Analysis  ** 

CS-TDMA  System  Parameters: 
m ,  n=  10  1 

the  mean  batch  arrival  rate  pi  1 . OOOOOOOOOOOOOOOE-O02 

the  mean  service  rate  ql  2 . OOOOOOOOOOOOOOOE-OO 1 

system  throughput  rho  =  9 . OOOOOOOOOOOOOOOE-OO 1 

the  batch  size  is  :  geometric 

the  mean  batch  size  is  :  1.800000000000000 


The  queue-size  probabilities  are: 


J 

P(j) 

D(j) 

0 

. 140707071 

.  000000000 

1 

.091717172 

.419871795 

2 

.  133434343 

.000000000 

3 

. 1 09292929 

. 000000000 

4 

.111717172 

. 003205128 

5 

. 096060606 

. 000000000 

6 

. 060505051 

.  000000000 

7 

.057575758 

.000000000 

8 

.  036262626 

.  000000000 

9 

.032828283 

.000000000 

10 

.014040404 

. 000000000 

1  1 

.017373737 

. 000000000 

12 

.010707071 

.000000000 

1  3 

.013535354 

. 000000000 

14 

.015353535 

.006410256 

1  5 

.014141414 

.  000000000 

16 

.011414141 

.  000000000 

1  7 

. 007676768 

. 000000000 

18 

. 002020202 

. 000000000 

19 

. 003737374 

.000000000 

20 

.012121212 

.000000000 

21 

.  005353535 

.0064  10256 

22 

.  002424242 

.003205128 

23 

.000000000 

.006410256 

24 

.  000000000 

.000000000 

25 

.000000000 

.003205128 

26 

.000000000 

. 000000000 

27 

. 000000000 

. 000000000 

28 

.  000000000 

. 000000000 

29 

. 000000000 

. 000000000 

30 

.000000000 

.  000000000 

31 

.000000000 

.000000000 

32 

.000000000 

. 000000000 

33 

.000000000 

.003205128 

34 

.  000000000 

.  000000000 

35 

.  000000000 

.003205128 

81 


36 

.000000000 

. 003205128 

37 

. 000000000 

.000000000 

38 

. 000000000 

.000000000 

39 

. 000000000 

. 000000000 

40 

.000000000 

. 003205128 

4  1 

. 000000000 

. 000000000 

42 

.000000000 

.000000000 

43 

.  000000000 

.  000000000 

44 

. 000000000 

.000000000 

45 

.  000000000 

. 000000000 

46 

.  000000000 

.  0032051  28 

47 

. 000000000 

.003205128 

48 

.000000000 

. 003205128 

49 

.000000000 

.000000000 

50 

.000000000 

. 003205128 

51 

.000000000 

.  000000000 

52 

.  000000000 

.000000000 

53 

.  000000000 

.003205128 

54 

.000000000 

.000000000 

55 

.  000000000 

.000000000 

56 

.000000000 

. 003205128 

57 

.  ooooocooo 

.003205128 

58 

.  000000000 

.003205128 

59 

.  000000000 

.000000000 

60 

.  000000000 

. 003205128 

61 

.  000000000 

.000000000 

62 

.000000000 

.003205128 

63 

. 000000000 

.  003205128 

64 

.000000000 

.000000000 

65 

. 000000000 

. 003205128 

66 

.000000000 

.000000000 

67 

. 000000000 

.006410256 

6e 

.000000000 

.012820513 

69 

. 000000000 

. 000000000 

70 

.000000060 

.003205128 

71 

. 000000000 

.000000000 

72 

. 000000000 

. 000000000 

73 

.  oocoooooo 

.  0032051  28 

74 

.  000000000 

.003205128 

75 

.  000000000 

.000000000 

76 

.000000000 

.000000000 

77 

.000000000 

.  0032051  28 

78 

. 000000000 

.003205128 

79 

.  000000000 

. 003205128 

80 

.000000000 

. 000000000 

81 

.000000000 

.  000000000 

82 

.  000000000 

.  000000000 

83 

.000000000 

.000000000 

84 

.000000000 

. 000000000 

85 

.  000000000 

.0064 10256 

86 

.  000000000 

.000000000 

87 

. 000000000 

.003205128 

88 

.000000000 

. 000000000 
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89 

.000000000 

. 000000000 

90 

.000000000 

. 000000000 

91 

.000000000 

.003205128 

92 

.  000000000 

.000000000 

93 

.000000000 

. 003205128 

94 

.000000000 

. 000000000 

95 

. 000000000 

.006410256 

96 

.000000000 

. 000000000 

97 

.000000000 

.006410256 

98 

.000000000 

.  000000000 

99 

.000000000 

.003205128 

>=100 

.  000000000 

.429487179 

mean  queue-size  (simulation)  is  : 

standart 

deviation  (queue-size,  simu 

4 . 856969696969697 
is:  4.608791216124809 


mean  delay  (simulation)  is:  153.708333333333300 

standard  deviation  (delay,  simulation)  is:  247.684721276557900 


P(X> 

10)  = 

. 1 1585858586 

P(D> 

10)  = 

. 57692307692 

mean 

queue- 

size  (analysis. 

1  ower 

bound)  is 

mean 

queue- 

size  (analysis. 

upper 

bound)  is 

mean 

del  ay 

(analysis,  lower 

bound ) 

i  S  : 

mean 

delay 

(analysis,  upper 

bound ) 

i  s  : 

15.228000000000010 

15.228000000000010 

846.000000000000200 
846 . 000000000000200 
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Delay-Throughput  Evaluator,  IRI  Corp. 


Appendix  D.  The  CS-TDMA  Program:  Source  Code 
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Circuit-Switched  TDMA  System 


C: 

c 

c 

c 

c 

c 

c 

c 

c 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C: 

C; 

C: 

C: 

C: 

C: 

C: 
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*  Version  1  Sept.  20,  1989  * 

*  Professor  Izhak  Rubin  * 

*  Department  of  Electrical  Engineering  * 

*  Univeristy  of  California  * 

*  Los  Angels,  CA  90024  * 

*  ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★■Jr 


Descriptions  of  the  variables: 


itbs 

Pi 

mu 


total  arrivals  during  a  frame 
mean  session  arrival  rate 
mean  session  service  rate 


types 

u 

ul 

u2 

seedt 

seedtl 

ic(  j) 

P(j) 

qmean 

External 


:  the  type  of  batch  size  distribution 
:  deterministic  batch  size 

:  lower  limit  of  the  uniformly  dist .  batch  size 
:  upper  limit  of  the  uniformly  dist.  batch  size 
:  the  seed  to  generate  the  batch  arrival 
:  the  seed  to  generate  the  number  of  session  served 
per  slot 

:  total  number  of  sessions  in  that  system  at  the  start  of 
jth  slot  ( j-0, 1, 2, . . . , 100) 

:  the  prob.  that  herer  are  j  sessions  in  the  system  at  the  start  of 
a  slot  ( j-0, 1, 2, . . . , 100) 

:  mean  queue-size  calculated  from  simulation 
function:  rand (seed) 

return  a  value  that  is  uniformly  distributed  over  (0,  1) 


implicit  real*8 (a-h, o-z) 

integer  xn(500) , seedt, seeds, types, batch, bs 
integer  icd (0 : 100, 1 : 500 ) , id (0 : 100 ) , dn (500 ) 

integer  ict (0:100, 1:500), ic (0:100), typesl, seedtl, index (500) 
real*8  mu 

common  mu, u, ul ,  u2 , seeds 
character*12  filen 

print  *,  'Enter  the  file_name  for  data:  ' 
print  *,  '(no  more  than  12  characters)' 
print  *,  '(in  PC,  QUOTE  the  "  f  ile_name"  )  ' 
read  *,  filen 
open (10, f ile-filen) 

print  *,  'What  would  you  like  to  do  ?  ' 

print  *,  '  1:  Simulation  only  ' 

print  *,  '  2:  Analysis  only  ' 

print  * ,  '  3:  Simulation  and  analysis  ' 

print  *,  'Please  enter  a  number:' 

read  *,  typesl 

if (typesl .eq. 1)  then 

write(10,*)  '  ' 

write(10,*)  '**  Simulation  **' 

print  *,  'Enter  the  start  time  (when  statistic  starts)' 
read  *,  istart 

print  *,  'Enter  the  stop  time  (total  simulation  time)' 
read  *,  istop 
elseif (typesl .eq. 3)  then 
write(10,*)  '  ' 

write  (10,*)  '**  Simulation  and  Analysis  **' 

print  *,  'Enter  the  start  time  (when  statistic  starts)' 

read  *,  istart 

print  *,  'Enter  the  stop  time  (total  simulation  time)' 
read  *,  istop 
'dseif  (typesl .eq. 2)  then 
write(10,*)  '  ' 

write(10,*)  '****  Analysis  ****' 

endif 


print  *,'  ' 
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c 


print  *,  'Enter  the  frame  duration  (m)  (<-500 )  and  number' 
print  *,  'of  slots  allocated  to  a  station  (n)  ' 
print  *,  '  n<-m,  m  and  n  should  be  integers' 
read  *,  m,  n 

print  *,  'Enter  the  batch  arrival  rate  pi  ' 
read  *,  pi 

print  *,  'Enter  the  transmission  rate  ql  ' 
read  *,  ql 

1  print  *,  'Enter  the  batch  size  distribution  index' 
print  *,  '  1:  deterministic' 

print  *,  '  2:  geometric' 

print  *,  '  3:  uniform' 

print  *,  '  ' 

read  *,  types 
if (types. eq.l)  then 

print  *,  'Enter  the  fixed  batch  size  u  (u*pl<n*ql/m) : ' 
read  *,  u 
rhO“*u*pl*m/ (n*ql) 
if  (rho.ge.l.dO)  then 

print  *,  'system  not  stable,  please  reenter  the  parameters  ' 
go  to  88 
endif 
b-u*pl 

a2-pl*u* (u-1) 
elseif (types. eq. 2)  then 

print  *,  'Enter  the  mean  batch  size,  b  (b*pl<n*ql/m) : ' 
read  *,  b 
rho-b*pl*m/ (n*ql ) 
if  (rho.ge.l.dO)  then 

print  *,  'system  not  stable  ,  please  reenter  the  parameters  ' 
go  to  88 
endif 

mu»l . OdO/b 
b-b*pl 

a2«pl*2 .d0* (1-mu) /mu**2 
elseif (types . eq . 3 )  then 

print  *,  'Enter  the  range  of  the  batch  size:' 

print  *,  '  (e.g.,  "1.0,  3.0")  (ul+u2) *pl/2<n*ql/m  ' 

read  *,  ul,u2 

rho~  (ul+u2) *pl*m/ (2*n*ql) 

if  (rho.ge.l.dO)  then 

print  *,  'system  not  stable,  please  reenter  the  parameters  ' 
go  to  88 
endif 

b«pl*(ul+u2)/2.d0 

a2«pl* (u2* (u2+l ) * (2*u2+l) /6-u2* (u2+l) 12- (ul-1) *ul* (2*ul-l) 

*  /6+ul* (ul-1) /2) / (u2-ul+l) 

else 

go  to  1 

endif 


'  CS-TDMA  System  Parameters:  ' 
'm,  n-  ',m, n 

'the  mean  batch  arrival  rate  pi 
'the  mean  service  rate  ql  '^ql 
'system  throughput  rho  -  ',rho 


wriw  (10,  *) 
write  (10, *) 
write (10, *) 
write (10, *) 
write  (10, *) 
if (types. eq.l)  then 

write (10,*)  'the  batch 
write (10,*)  'the  fixed 
elseif (types. eq. 2)  then 

write (10,*)  'the  batch  size 
write (10,*)  'the  mean  batch 
elseif (types. eq. 3)  then 

write  (10,*)  'the  batch  size  is  :  uniform  ' 
write  (10,*)  'the  range  of  the  batch  size  is:  l',ul, 
endif 


size  is  : 
batch  size 


is  : 
size 


deterministic 
is  :  ',u 

geometric  ' 
is  :  ',1/mu 


u2 , '  ] ' 
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15 

14 

c 


16 

117 

c 

c  ** 
c 


149 

c 

119 


if (typesl .eq.2)  goto  125 

print  *,  'To  find  the  probabilities  P(X>x)  and  P(D>d),  please' 

print  *,  'enter  x  and  d  (  x,  d<-100,  both  should  be  integers)' 

read{5,*)  igx,igd 

print  *,  'please  wait' 

seedt-1099 

seedtl-1099 

seeds-99999 

nun-100 

do  14  i— 0 , mm 
do  15  j-l,m 
ict (i, j ) -0 
icd (i, j ) -0 
ic (i) -0 
id (i) -0 
continue 

totalx-0 . OdO 

tolx2-0 . OdO 

ktd-0 

tdd-0 

tdd2-0 

px-0 

pd-0 

do  16  i-l,m 
index (i) -0 
xn (i)-0 

do  117  i-1,500 
dn (i)-0 

simulation  starts  from  here 
npl*=n+l 

do  20  i-l,istop 
itbs-0 

mn«nb ( seedt , pi ) 
if  (index (1) .eq. 0)  then 
if(mn.eq.l)  then 
bs-batch (types) 
xn (1) =xn (m) +bs 
itbs-itbs+bs 
else 

xn (1) -xn (m) 
endif 
else 

ndn-nb (seedtl, ql) 
if(ndn.eq.l)  index(l)-0 
if(mn.eq.O)  then 
xn (l)-xn (m) -ndn 
else 

bs-batch (types ) 
xn (1) -xn (m) -ndn+bs*mn 
itbs-itbs+bs 
endif 
endif 

11-xn (1) +ndn 
do  149  kd-1,11 
dn (kd)-dn (kd) +1 

if(ndn.eq.l)  then 
kdd-dn ( 1 ) 
do  119  kd-1, xn (1 ) 
dn (kd) -dn (kd+1 ) 
dn (xn(l) +l)-0 
end?  f 
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if (i . gt . istart )  then 
if (xn (1) . gt .mm)  then 
ict (mm, 1 ) -ict (mm, 1 ) +1 
else 

ict  (xn  (1) ,  1)  -ict  (xn  (1) ,  1)  +1 
endif 

totalx-totalx+xn ( 1 ) 
tolx2-tolx2+xn (1) **2 
if(ndn.eq.l)  then 
ktd-ktd+1 
tdd-tdd+kdd 
tdd2-tdd2+kdd**2 
if (kdd. gt .mm)  then 
icd (mm, 1 ) -icd (mm, 1 ) +1 
else 

icd(kdd, 1) -icd (kdd, 1) +1 
endif 
endif 
endif 
c 

do  17  j-2,n 
mn-nb ( seedt , pi ) 
if  (index (j) .eq.O)  then 
if(mn.eq.l)  then 
bs=batch (types) 
xn  ( j) -xn ( j  —  1 ) +bs 
itbs-itbs+bs 
else 

xn ( j) -xn ( j— 1 ) 
endif 
else 

ndn-nb ( seedt 1 , ql ) 
if(ndn.eq.l)  index (j)«0 
if (mn. eq.O)  then 
xn ( j)-xn ( j-1) -ndn 
else 

bs-batch (types) 
xn ( j)-xn ( j-1) -ndn+bs*mn 
itbs-itbs+bs 
endif 
endif 

11-xn ( j) +ndn 
do  139  kd-1,11 
139  dn (kd) -dn (kd) +1 

c 

if(ndn.eq.l)  then 
kdd-dn (1) 
do  129  kd«l,xn(1) 

129  dn (kd) -dn (kd+i) 

dn(xn(  j)+l)-0 

endif 

c 

if (i . gt . istart)  then 
if (xn ( j) .gt .mm)  then 
ict (mm, j ) -ict (mm, j ) +1 
else 

ict (xn( j) , j) -ict (xn ( j) , j) +1 
endif 

totalx-totalx+xn ( j ) 
tolx2-tolx2+xn ( j)  **2 
if(ndn.eq.l)  then 
ktd-ktd+1 
tdd-tdd+kdd 
t  dd2- 1  dd2  +kdd  *  *  2 
if (kdd. gt .mm)  then 
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icd (mm, 1 ) -icd (mm, 1 ) +1 
else 

icd (kdd, 1 ) -icd (kdd, 1 ) +1 
endif 
endif 
endif 

17  continue 
c 

do  18  j-npl,m 
mn-nb ( seedt , pi ) 
if(mn.eq.l)  then 
bs-batch (types) 
xn( j)-xn( j-l)+bs 
itbs=itbs+bs 
else 

xn( j)-xn( j-1) 
endif 

do  211  kd-l,xn(j) 

211  dn (kd) -dn (kd) +1 

if (i.gt .istart)  then 
if (xn ( j) .gt .mm)  then 
ict (mm, j ) -ict (mm, j ) +1 
else 

ict (xn ( j) , j) -ict (xn ( j) , j)+l 
endif 

totalx-totalx+xn ( j ) 
tolx2«tolx2+xn ( j )  **2 
endif 

18  continue 

if (xn (m) . ge . n)  then 
do  78  k-l,n 
78  index (k)«l 
else 

do  23  k-l,n 
if (itbs.gt.O)  then 
if (index (k) . eq.O)  then 
index (k) -1 
itbs-itbs-1 
endif 
else 
goto  20 
endif 

23  continue 
endif 

20  continue 

do  22  j-0,mm 
do  21  k-1 , m 
ic ( j)-ic ( j)+ict ( j,k) 
id ( j ) -id ( j ) +icd ( j , k ) 

21  continue 

22  continue 
c 

last-istop-istart 

print  *,'The  queue-size  probabilities  are  as  follows:' 


write  (6, *)  'j  P(j)  D ( j) 

print  *,' - ' 

write  (10,*)  'The  queue-size  probabilities  are:' 

write  (10,  *)  'j  P(j)  D  ( j ) 

write(10,*)  ' - ' 


a.uean-totalx/  (last*m) 

var-dsqrt (  tolx2/ (last*m)  -qmean**2  ) 

wmean-tdd/ (ktd*1.0d0) 

vard-dsqrt(  tdd2/ (ktd*l .OdO) -dmean**2  ) 

do  110  i-0,mm 

if (i. It .mm)  then 
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write(6,200)  i, ic (i) *1 .dO/ (m*last) , id(i) *1 .dO/ktd 
write {10, 200)  i,ic (i)*l.dO/ (m*last) , id (i) *1 .dO/ktd 
else 

write  (6,  210)  i,ic(i)  *1  .dO/  (m*last)  ,id(i)  *1. dO/ktd 
write  (10, 210)  i,  ic (i) *1 . dO/ (last *m) , id (i) *1 .dO/ktd 
endif 

if (i.gt.igx)  px-px+ic (i) *1 . 0/ (m*last) 
if(i.gt.igd)  pd=pd+id (i) *1 .dO/ktd 
110  continue 

200  format (lx, 13, 2x, 2f 14 . 9) 

210  format (lx,' >-' , i3, 2x, 2f 14 . 9) 

write  (6,*)  'mean  queue-size  (simulation)  is  :  ',qmean 

write  (6,*)  'standart  deviation  (queue-size,  simu.)  is  : 

write  (6,*)  'mean  delay  (simulation)  is:  ',wmean 

write  (6,*)  'standart  deviation  (delay,  simulation)  is:  ', 

write (10,*)  'mean  queue-size  (simulation)  is  :  ',qmean 

write  (10,*)  'standart  deviation  (queue-size,  simu.)  is:  ' 

write  (10,*)  'mean  delay  (simulation)  is:  ' , wmean 

write  (1'.,*)  'standard  deviation  (delay,  simulation)  is:  ' 

write (6,300)  igx,px 

write (6, 301)  igd,pd 

write  (10, 300)  igx,px 

write  (10, 301)  igd,pd 

300  format (lx, 'P(X>',i3,' )-  ',fl4.11) 

301  format (lx, 'P (D>' , i3, ' )-  ',fl4.11) 
if (typesl .eq. 1)  goto  126 

c  *****  analysis 

125  amean-  m* (  m*b+ (m* (m-1) *b**2+m*a2- (n-1) *m*b*ql 

*  +2*m*b* (1-ql) ) / (2* (n*ql-m*b) )  ) 

write  (6,*)  'mean  queue-size  (analysis,  lower  bound)  is  : 

*  amean/m- (m-1 ) *b/2 

write  (6,*)  'mean  queue-size  (analysis,  upper  bound)  is  : 

*  amean/m- (m-1) *b/ 2  + (n-1) * (2-ql) /2+  (n-1  ) *b 
write(6,*)  'mean  delay  (analysis,  lower  bound)  is:  ', 

*  amean/ (m*b)- (m-1) /2.d0 

write (6,*)  'mean  delay  (analysis,  upper  bound)  is:  ', 

*  (amean+m* (n-1) * (2-ql) /2) / (m*b) -  (m-l)/2.d0  + 
write (10,*)  'mean  queue-size  (analysis,  lower  bound)  is  : 

*  amean/m- (m-1 ) *b/ 2 

write (10,*)  'mean  queue-size  (analysis,  upper  bound)  is  : 

*  amean/m  -(m-1) *b/2  + (n-1) * (2-ql) /2+  (n-1  ) *b 
write(10,*)  'mean  delay  (analysis,  lower  bound)  is:  ', 

*  amean/ (m*b)-(m-l)/2.d0 

write  (10,*)  'mean  delay  (analysis,  upper  bound)  is:  ', 

*  (amean+m* (n-1) * (2-ql) /2) / (m*b) -  (m-l)/2.d0  + 
126  close(10) 

stop 

end 


integer  function  batch (types) 
integer  seeds, types 
real*8  mu,u,ul,u2,rand 
common  mu, u,ul,u2, seeds 
if (types. eq.l)  then 
batch-int (u*l .001) 
elseif (types. eq. 2)  then 

batch-igeom (1 . OdO/mu, seeds) 

else 

batch-ul+int (  (u2-ul+l) *rand (seeds) ) 
endif 
return 
end 

C - Random  Generator  Function - 


' ,  var 
vard 
,  var 
,  vard 


9 

9 

9 

9 

n-1 

9 

9 

9 

9 

n-1 


91 


nonno 


This  function  will  take  the  argument  "seed"  as  an  "input"  seed, 

and  return  the  value  which  is  uniformly  distributed  over  [0,  1], 
and  the  output  seed  to  be  used  as  the  next  input  seed, 
function  rana(seed) 
integer  seed 

double  precision  pp, i7, i2, rand 

i7-7**5 

i2-2**31-l 

pp-i7*seed 

seed-dmod (pp,  i2 ) 

rand-seed/ (i2+l) 

return 

end 


FUNCTION  IGEOM  **************************** 


FUNCTION  IGEOM  PRODUCES  A  GEOMETRICALLY  DISTRIBUTED  PSEUDO-RANDOM 
OBSERVATION  WITH  AVERAGE  AVG  ,  FROM  R.N.  STREAM  ISTRM 
FUNCTION  IGEOM (AVG, ISTRM) 

REAL*8  THETA, GI, FI, avg, uu, rand 
Uu-  RAND(istrm) 

THETA  -  1 . OdO/AVG 
NN-INT (10 . 0d0*AVG) 

GI-THETA 
FI-THETA 
1-1 

IF (Uu.LE.FI)  GOTO  20 
DO  10  1-2, NN 
GI  -  ( 1. -THETA) *GI 
FI  -  FI  +  GI 
IF (Uu.LE.FI)  GOTO  20 
10  CONTINUE 
C  WRITE (*,6000)  TNOW, AVG, Uu 

C6000  FORMAT (IX, ' ! ! ! !  ERROR  IN  %  IGEOM  %,  TIME  -  ',F12.4, 
c  */2X, '  AVG  — ' , F12 .  4 ,  '  Uu  -  ',F8.5) 

20  IGEOM  -  I 
RETURN 
END 
C 


c  This  function  will  take  the  value  1  with  proba.  p  and  0  with  prob.  1-p 
function  nb(seed,p) 
real*8  t,p,rand 
integer  seed 
t-rand(seed) 
if(t.le.p)  then 
nb-1 
else 
nb-0 
endif 
return 
end 
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ITDMA  EXAMPLE  1 


**  Simulation  and  Analysis  ** 


Integrated  CS/PS  System  Parameters: 


m , 

n  = 

1  0 

5 

t  he 

mean 

batch  arrival  rate  pi  1 . 000000000000000E-00 1 

t  he 

mean 

service  rate  ql 

1 .000000000000000 

system  throughput  rho  = 

6. 522882542734397 E- 001 

t  he 

batch  size  is  ;  geometric 

t  he 

mean 

batch  size  is  : 

5 . 000000000000000 

The 

queue-size  probabilities  are: 

j 

P(J) 

D(j) 

0 

.300333333 

. 000000000 

1 

. 04 1666667 

.000000000 

2 

. 040666667 

.  049743475 

3 

. 037444444 

.046174437 

4 

. 034888889 

. 04 1 043944 

5 

.032555556 

. 040374749 

6 

. 032555556 

.  040374749 

7 

.  027555556 

. 040374749 

8 

. 028777778 

. 036805710 

9 

.030555556 

.038367165 

10 

. 028  1  1  1  1  1  1 

.033236672 

1  1 

.026222222 

. 030336828 

12 

.024777778 

.030113763 

13 

.019555556 

.028106179 

14 

.  024000000 

.027660049 

15 

.017333333 

.027213919 

16 

.019555556 

.024760205 

17 

.014333333 

.023644881 

1  8 

.  0  1  6333333 

.  024314075 

1  9 

.013222222 

. 022306491 

>=  20 

.  189555556 

.395047959 

mean 

sys 

tem-size  (simulation)  is  :  10.965111111111110 

standard 

deviation  (system-size,  simu.)  is:  14.562551077216630 

mean 

delay  (simulation)  is:  22.974124470220830 

standard 

deviation  (delay 

,  simulation)  is:  23.174030509119400 

P(X> 

10) 

=  .36488888889 

P(D> 

10) 

=  .63350434977 

Blocking 

probability  of  a 

CS  session:  7 . 2 9 3 2 94 7 38 84 7 554 E-002 

mean 

queue-size  (up  and  lower  bounds) 

72 

.280568408166620 

4 . 500000000000000 

mean  delay  (up  and  lower 

bounds  ) 

144 

.561 136816333200 

9 . 000000000000000 
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C> i tdma  2 

Enter  the  file_name  for  data 
(no  more  than  12  characters) 

(in  PC,  QUOTE  the  ’file_name’) 

’ i out  1  ' 

Wnat  would  you  like  to  do  ? 

1 :  Simulation  only 
2:  Analysis  only 
3:  Simulation  and  analysis 
Please  enter  a  number: 

3 

Enter  the  start  time  (when  statistic  starts) 

100 

Enter  the  stoo  time  (total  simulation  time) 

1000 

Enter  the  frame  duration  (m)  (m<=50)  and  the  maximum 
number  of  slots  for  CS  support  (n) 

(n<=m)  m  and  n  should  be  integers 
10  5 

Enter  the  session  (circuit)  arrival  rate 

0.01 

Enter  the  session  transmission  rate 
0.04 

Enter  the  packet  batch  arrival  rate  [mess. /slot]  pi 
0  .  1 

please  wait 

average  number  of  CS  slots  used  per  frame  (V):  2.334677242395972 

Enter  the  batch  size  distribution  index 
1 :  deterministi c 
2:  geometric 
3:  uniform 

2 

Enter  the  mean  batch  size,  b  (m*b*p1 <m-V) : 

5 

To  find  the  probabilities  P(X>x)  and  P(D>d),  please 
enter  x  and  d  (x,  d<=  20.  both  should  be  integers) 

10  10 
please  wait 


Delay-Throughput  Evaluator,  IRI  Cotp. 
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oooonooooooooooooooooo 


Integrated  CS/PS  TDMA  System 


C: 
c 

c  *  Version  1  Sept.  20,  1989  * 

c  *  Professor  Izhak  Rubin  * 

c  *  Department  of  Electrical  Engineering  * 

c  *  Univeristy  of  California  * 

c  *  Los  Angels,  CA  90024  * 

Q  ★★★★★★**★★★*★★★★★★★★★★*★★★★★★★★★★★★★★★★★*★★★ 

C 

c  This  program  finds  the  mean  system  size  and  delay  for 
c  a  movable  boundary  integrated  data/vioce  system.  For 
c  the  CS  subsystem,  an  M/Geom/N/N  queueing  model  is  assumed, 
c  While  for  the  PS  part,  the  data  transmission  rate 
c  is  assumed  one. 


c 


Descriptions  of  the  variables: 
m  :  frame  duration 

n  :  maximum  number  of  slot  for  CS  support 
n<-m 

q(i, j)  :  transition  matrix  of  the  CS  connection  system 
aim  :  mean  session  arrival  rate 
uc  :  mean  session  transimission  rate 
pi  :  mean  packet  arrival  rate 

ql  :  mean  packet  transimssion  rate,  here  we  set  ql-1 
types  :  the  type  of  service  time  distribution 
u  :  deterministic  batch  size 

ul  :  lower  limit  of  the  uniformly  dist.  batch  size 

u2  :  upper  limit  of  the  uniformly  dist .  batch  size 

seedt  :  the  seed  to  generate  the  number  of  batch  arrivals  per  slot 

seeds  :  the  seed  . . .  service 

ic(j)  :  total  number  of  packets  in  the  system  at  the  start  of  jth  slot 
< j-0,1,2, . . .,20) 

p(j)  :  the  prob.  that  there  are  j  packets  in  the  system 
at  the  start  of  a  slot  ( j-0, 1, 2, . . . , 20) 
qmean  :  mean  queue-size  obtained  from  simulation 
External  function:  rand {seed) 

return  a  value  that  is  uniformly  distributed  over  (0,  1) 

imolicit  real*8 (a-h,o-z) 

integer  xn (50) , types, ict (0 : 20, 1 : 50) , ic (0 : 20) , typesl 

integer  seed, seeds, icd (0:20, 1:50) ,id(0:20) ,dn(500) 

dimension  a(0: 50,0: 50 ),q(0:50, 0:50), g(0:50) 

dimension  x(0:50),v(0:50)p(0:20),dw(0:20) 

common  amu, u, ul,u2, seeds, n,m, nl,ml, seed, istart, istop 

character*12  filen 

nl-50 

ml-50 

mm- 20 

print  * ,  'Enter  the  file_name  for  data  ' 

print  * ,  ' (no  more  than  12  characters)  ' 

print  *,  '(in  PC,  QUOTE  the  "  f  ile_name"  ) ' 

read  *,  filen 

open (10, f ile-filen) 

ql-1 

print  *,  'What  would  you  like  to  do  ?  ' 

print  '  1:  Simulation  only  ' 

print  *,  '  2:  Analysis  only  ' 

print  *,  '  3:  Simulation  and  analysis  ' 

print  *,  'Please  enter  a  number:' 

read  *,  typesl 

if (typesl .eq. 1)  then 

write(10,*)  '  ' 

write (10,*)  '**  Simulation  **' 

print  *,  'Enter  the  start  time  (when  statistic  starts)' 
read  *,  istart 
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print  *,  'Enter  the  stop  time  (total  simulation  time)' 
read  *,  istop 
elseif (typesl .eq.3)  then 
write (10,*)  '  ' 

write (10,*)  '**  Simulation  and  Analysis  **' 

print  *,  'Enter  the  start  time  (when  statistic  starts)' 

read  *,  istart 

print  *,  'Enter  the  stop  time  (total  simulation  time)' 
read  *,  istop 
elseif  (typesl .eq. 2)  then 
write(10,*)  '  ' 

write(10,*)  '****  Analysis  ****' 

endif 

c 

88  print  *,  '  ' 

print  *,  'Enter  the  frame  duration  (m)  (m<-50)  and  the  maximum' 

print  *,  '  number  of  slots  for  CS  support  (n)  ' 
print  *,  '  (n<-m)  m  and  n  should  be  integers' 
read  *,  m, n 

print  *,  'Enter  the  session  (circuit)  arrival  rate  ' 
read  *,  aim 

print  *,  'Enter  the  session  transmission  rate  ' 
read  *,  uc 

print  *,  'Enter  the  packet  batch  arrival  rate  (mess . /slot ]  pi  ' 
read  *,  pi 

c 

print  *,  'please  wait  ' 
c  **begin  of  computing  v(i) 

call  qmat (q, v,m, n,alm,uc,nl,ml) 
delt  -  1.0d-14 
call  vi (a,q,n,nl,ml,delt, v) 
sum-0 .d0 
do  215  i-1, n 

215  sum=sum+i*v(i) 

vbar-sum 

write(6,*)  'average  number  of  CS  slots  used  per  frame  (V):',vbar 
c  **  end  of  computing  sum  of  v(i) 
c 

1  print  *,  'Enter  the  batch  size  distribution  index' 
print  *,  '  1:  deterministic' 

print  *,  '  2:  geometric' 

print  *,  '  3:  uniform' 

print  *,  '  ' 

read  *,  types 
if (types .eq. 1)  then 

print  *,  'Enter  the  fixed  batch  size  u  (m*u*pl<m-V) : ' 
read  *,  u 
rho-u*pl*m/ (m-sum) 
if  (rho . ge . 1 .d0)  then 
print  *,  '  ' 

print  *,  'system  not  stable,  please  reenter  the  parameters' 
go  to  88 
endif 
b-u*pi 

a2-pl*u* (u-1 ) 
elseif (types. eq. 2)  then 

print  *,  'Enter  the  mean  batch  size,  b  (m*b*pl<m-V) : ' 
read  *,  b 
rho-b*pl*m/ (m-sum) 
if  (rho.ge.l.dO)  then 
print  *,  '  ' 

print  *,  'system  not  stable,  please  reenter  the  parameters' 
go  to  88 
endif 

amu-1 . OdO/b 
b-b*pl 
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a2«pl*2 .dO* (1-amu) /amu**2 
elseif (types. eq. 3)  then 

print  *,  'Enter  the  range  of  the  batch  size:' 
print  *,  '  (e.g.,  "1.0,  3.0")  m* (ul+u2) *pl/2 (m-V)  ' 

read  *,  ul,u2 

rho- (ul+u2) *pl*m/ (2* (m-sum) ) 
if  (rho.ge.l .dO)  then 
print  *,  '  ' 

print  *,  'system  not  stable,  please  reenter  the  parameters' 
go  to  88 
endif 

b-pl*(ul+u2)/2.d0 

a2-pl* (u2* (u2+l) * (2*u2+l) /6-u2* (u2+l) /2- (ul-1) *ul* (2*ul-l) 

/ 6+ul* (ul-1) /2) / (u2-ul+l) 

else 

go  to  1 

endif 


c  **  write  the  parameters  to  the  output  file 

write (10,*)  '  Integrated  CS/PS  System  Parameters:  ' 
write (10,*)  '  m,  n  -  ',m,n 

write (10,*)  'the  mean  batch  arrival  rate  pi  ',pl 
write (10,*)  'the  mean  service  rate  ql  ',q 1 

write  (10,*)  'system  throughput  rho  -  ' , rho 
if (types. eq.l)  then 

write (10,*)  'the  batch  size  is  :  deterministic  ' 
write (10,*)  'the  fixed  batch  size  is  :  ',u 
elseif (types. eq. 2)  then 

write (10,*)  'the  batch  size  is  ;  geometric  ' 
write (10,*)  'the  mean  batch  size  is  :  ',1/amu 
elseif (types. eq. 3)  then 

write (10,*)  'the  batch  size  is  :  uniform  ' 
write(10,*)  'the  range  of  the  batch  size  is:  [' ,ul, ' : ' ,u2, 
endif 


if (typesl .eq.2)  then 

write (6,*)  'please  wait' 
goto  127 
endif 

print  *,  'To  find  the  probabilities  P(X>x)  and  P(D>d),  please' 
print  *,  'enter  x  and  d  (x,  do  20,  both  should  be  integers)' 
read (5,*)  igx,igd 
print  *,  'please  wait' 
c  **  start  simulation 

call  simu (pl,p, qmean, q, xn, types, ic, ict, var, 

*  dn,wmean,dw,vard,icd, id, mm) 


c 


pd-0 

px-0 

print  *,'The  queue-size  probabilities  are  as  follows:' 


write  (6, *)  'j  P(j)  D ( j) 

print  *,' - ' 

write  (10,*)  'The  queue-size  probabilities  are:' 
writedO,*)  'j  P(j)  D(j) 

write  CO,*)  ' - : - ' 

do  110  i-0, mm 
if (i. It .mm)  then 


write(6,200)  i,p(i),dw(i) 
write (10, 200)  i,p(i),dw(i) 
else 

write(6,210)  i,p(i),dw(i) 
writedO, 210)  i,p(i),dw(i) 
endif 

if(i.gt.igx)  px-px+p(i) 
if(i.gt.igd)  pd-pd+dw(i) 

110  continue 


100 


200  format (lx,  i3, 2x,  f  14 . 9, 2x,  f  14 . 9) 

210  format (lx, '  >-' ,  i3, 2x,  fl4 . 9,  2x,  f  14 . 9) 

write(6,*)  'mean  system-size  (simulation)  is  :  ',qmean 
write (6,*)  'standard  deviation  (system-size,  simu.)  is  :  ',var 
write(6,*)  'mean  delay  (simulation)  is:  ',wmean 
write  (6,*)  'standard  deviation  (delay,  simulation)  is:  ',vard 
write (10, *)  'mean  system-size  (simulation)  is  :  ' , qmean 

write  (10,*)  'standard  deviation  (system-size,  simu.)  is:  ',var 
write (10,*)  'mean  delay  (simulation)  is:  ' , wmean 
write (10,*)  'standard  deviation  (delay,  simulation)  is:  ',vard 
write (6, 50)  igx,px 
write (10, 50)  igx,px 
write (6, 301)  igd,pd 
write (10, 301)  igd,pd 
50  format (lx, 'P(X>',i3,')-  ',fl4.11) 

301  format (lx, ' P (D>' , i3, ' ) -  ',fl4.11) 
if (typesl .eq. 1)  goto  126 
c 

c  **  analysis  here 
c 

127  p2-a2 

pd-b 

rho-1- (m-m*pd-vbar ) 
npl-n+1 
do  20  i-0,n 
do  21  j-0,n 

21  a(i,  j)—  q(i,  j)+v(  j) 

20  continue 

c 

do  22  i-0, n 
a(i,i)-l+a(i,i) 
a(i,npl)=m*pd- (m-l)+i 
22  continue 

call  gauss (a,x,n,npl,nl,ml,delt) 
c 

do  25  i->0, n 

25  g(i)-0 
c 

do  26  i-0,n 
do  27  j-0,n 

27  g(i)-g(i)+q(i,  j)*x(j) 

26  continue 
c 

c  write (6, *)  'g(i)  1  ' , (g(i) ,i-0,n) 

vh-0 . 0 
do  30  i-0,n 

30  vh-vh+v (i) * (m* (1-pd) -i) *g (i) 
c 

call  smax (g,n, sm, smx) 
c 

v2-0.0 
do  29  i-2,n 

29  v2-v2+i*(i-l)*v(i) 

c 

vh2-m* (m-1) * (pd-1) **2  +m*p2  +2* (m*pd-m+l)  *vbar  +v2 
sum-0 

do  31  i-0,n 

31  sum-s\im->  (m-i-1)  *  (m-i)  *v  (i) 
c 

up-m*pd- (m-1)  +(m-l)*pd/2  -sum/(2*m)  +  vh2/ (2* (1-rho) ) 
bll-up-vh/ (1-rho)  +  sm 

upl-up-vh/ (1-rho) + (1-pd) * (m-1) + (m-1) * (1-rho) /m  +  smx 
c 

c  **  to  find  the  maximum  of  two  lower  bounds 
bl2-pd+p2/ (2* (1-pd) ) 
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if (bl2 . gt . bll )  bll-bl2 
c  **  to  check  if  it  is  locally  stable 
if (m*pd.lt .m-n)  then 
up2- (2** (1-pd) +p2) 121 (1-pd) 

*  +n* (p2-pd**2+pd) 121 (1-pd) / (m-n-m*pd) 
up2-up2+ (m-n-1) *n/2/m 

if (up2 .It .upl)  upl-up2 
endif 

if (m*pd.gt .m-n)  write(6,*)  'not  locally  stable' 
write (6,*)  'Blocking  probability  of  a  CS  session:  ',v(n) 
write (10,*)  'Blocking  probability  of  a  CS  session:  ',v(n) 
write (6,*)  '  mean  queue-size  (up  and  lower  bounds)' 
write (6,  *)  upl,'  ',bll 

write(6,*)  '  mean  delay  (up  and  lower  bounds)' 
write  (6,*)  upl/b, '  ',bll/b 

write (10,*)  '  mean  queue-size  (up  and  lower  bounds)  ' 
write  (10,*)  upl,'  ',bll 

write (10,*)  '  mean  delay  (up  and  lower  bounds)' 
write (10,*)  upl/b,'  ',bll/b 
c 

126  close (10) 
stop 
end 
c 
c 

SUBROUTINE  smax (x, n, sm, smx) 

IMPLICIT  REAL *8  (A-H,0-Z) 

DIMENSION  x  (0 : *) 
sm-x ( 0 ) 
smx-x (0) 

DO  10  I-l,n 
if  (x(i) .gt.smx)  then 
smx-x (i) 
else 

if  (x(i) .lt.sm)  sm-x(i) 
endif 

10  CONTINUE 
RETURN 
END 
c 
c 

subroutine  simu (pl,p, qmean, q, xn, types, ic, ict, var, 

*  dn,wmean,dw,vard,icd, id, mm) 
implicit  real*8 (a-h,o-z) 

integer  seed, seeds, seedt, seedtl, seeds, xn(l) , ic (0:mm) 
integer  ict (0 :mm,  1 :ml) , types, batch, bs,dn (1) , icd (0 :mm, 1 :ml) 
dimension  q (0 :nl, 0 :ml) , p(0 :mm) , dw (0 :mm) , id (0 :mm) 
common  amu, u, ul,u2, seeds, n,m, nl,ml, seed, istart, istop 
seed-1099 
seedt-1099 
seedtl-1099 
seeds-99999 
do  14  i-0,mm 
do  15  j-l,m 
ict (i, j)-0 
15  icd(i,j)-0 
ic(i)-0 
id(i)-0 
14  continue 
c 

totalx-0 . OdO 

tolx2-0 .0d0 

ktd-0 

tdd-0 

tdd2-0 

do  16  i-l,m 
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16 


xn(i)-0 

do  117  i-1,500 
117  dn (i)-0 
ni-0 
c 

c  *  siumlation  starts  from  here 
c 

do  20  i-l,istop 
c 

nj«nij (seed,n,ni,q,nl,ml) 
ni-n  j 
mmni-m-ni 
mmnil-mmni+l 
c 

mn-nb ( seedt , pi ) 
if (xn (m) .eq. 0)  then 
ndn=0 

if(mn.eq.O)  then 
xn (1) -0 
else 

bs-batch (types ) 
xn (1)  -mn*bs 
endif 
else 
ndn-1 

if(mn.eq.O)  then 
xn (1) -xn (m) -1 
else 

bs-batch (types) 
xn (1) -xn (m) -l+bs*mn 
endif 
endif 

11-xn (1) +ndn 
do  149  kd-1,11 
149  dn (kd) -dn (kd) +1 
c 

if(ndn.eq.l)  then 
kdd-dn ( 1 ) 
do  119  kd-1, xn (1) 

119  dn (kd) -dn (kd+1) 

dn (xn (1) +l) -0 
endif 
c 

if (i .gt . istart )  then 
if (xn (1) . gt .mm)  then 
ict (mm, 1) -ict (mm, 1) +1 
else 

ict (xn (1) , l)-ict (xn (1) , 1) +1 
endif 

totalx-totalx+xn ( 1 ) 
tolx2-tolx2+xn (1) **2 
if(ndn.eq.l)  then 
ktd-ktd+1 
tdd-tdd+kdd 
tdd2-tdd2+kdd**2 
if (kdd.gt .mm)  then 
icd (mm, 1 ) -icd (mm, 1 ) +1 
else 

icd(kdd, 1) -icd(kdd, 1) +1 
endif 
endif 
endif 
c 

do  17  j-2,mmni 
mn-nb ( seedt ,  pi ) 
if (xn ( j-1) .eq. 0)  then 
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ndn-0 

if(mn.eq.O)  then 
xn(j)-0 
else 

bs-batch (types) 
xn( j)«mn*bs 
endif 
else 
ndn-1 

if(nm.eq.O)  then 
xn(j)-xn( j-l)-l 
else 

bs-batch (types) 
xn(j)-  xn ( j-1) -l+bs*mn 
endif 
endif 
c 

11-xn ( j) +ndn 
do  139  kd-1,11 
139  dn(kd)-dn(kd)+l 
c 

if(ndn.eq  1)  then 
kdd-dn ( 1 ) 
do  129  kd-l,xn(j) 

129  dn (kd) -dn (kd+1) 

dn(xn( j)+l)-0 
endif 
c 

if (i.gt .istart)  then 
if  (xn  ( j)  .gt  .non)  then 
ict (mm, j ) -ict (mm, j ) +1 
else 

ict (xn ( j) , j)-ict (xn ( j ) , j) +1 
endif 

totalx-totalx+xn ( j ) 
tolx2-tolx2+xn ( j ) **2 
if(ndn.eq.l)  then 
ktd-ktd+1 
tdd-tdd+kdd 
tdd2-tdd2+kdd**2 
if (kdd. gt .mm)  then 
icd(mm, l)-icd(mm, 1) +1 
else 

icd (kdd, 1 ) -icd (kdd, 1 ) +1 
endif 
endif 
endif 
c 

17  continue 
c 

if (mmnil .gt .m)  goto  20 
do  18  j-mmnil,m 
mn-nb (seedt,pl) 
bs-batch (types) 
xn( j)-xn( j-1) +mn*bs 
do  211  kd-l,xn( j) 

211  dn(kd)-dn(kd)+3 

c 

if (i.gt .istart)  then 
if (xn( j) .gt.mm)  then 
ict (mm, j ) -ict (mm, j ) +1 
else 

ict (xn ( j) , j ) —ict (xn ( j) , j) +1 
endif 

totalx-totalx+xn ( j ) 
tolx2-tolx2+xn  ( j) **2 
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endif 


18  continue 

20  continue 

last-istop-istart 
qmean-totalx/ (last*m) 

var-dsqrt  (  tolx2/ (last*m)  -qmean**2  ) 
wmean-t dd/ ( kt d* 1 . OdO ) 

vard-dsqrt (  tdd2/ (ktd*1.0d0)-wmean**2  ) 

do  22  j«0,mm 
do  21  k-1 ,  m 
ic(j)-ic(j)+ict( j,k) 
id ( j)-id ( j) +icd(  j, k) 

21  continue 

p ( j ) -ic ( j ) * 1 . dO / ( last  *m) 
dw ( j ) —id ( j ) *1 .dO/ktd 

22  continue 

return 

end 


function  batch (types) 
integer  seeds, types, batch, seed 
real*8  amu, u, ul, u2, rand 

common  amu, u, ul,u2, seeds, n,m, nl,ml, seed, i start, istop 
if (types. eq.l)  then 
batch-int (u*l . 001 ) 
elseif (types.eq.2)  then 

batch-igeom(l .OdO/amu, seeds) 

else 

batch-int (ul) +int (  (u2-ul+l) *rand (seeds) ) 
endif 
return 
end 

- Random  Generator  Function - 

This  function  will  take  the  argument  "seed"  as  an  "input"  seed, 

and  return  the  value  which  is  uniformly  distributed  over  [0,  1], 
and  the  output  seed  to  be  used  as  the  next  input  seed, 
function  rand (seed) 
integer  seed 

double  precision  pp, i7, i2, rand 

i7— 7**5 

i2-2**31-l 

pp-i7*seed 

seed-dmod (pp, i2 ) 

rand-seed/ (i2+l) 

return 

end 


***********************  FUNCTION  IGEOM  **************************** 

FUNCTION  IGEOM  PRODUCES  A  GEOMETRICALLY  DISTRIBUTED  PSEUDO-RANDOM 
OBSERVATION  WITH  AVERAGE  AVG  ,  FROM  R.N.  STREAM  ISTRM 
FUNCTION  IGEOM (AVG, ISTRM) 

REAL* 8  THETA, GI, FI, avg,uu, rand 
Uu-  RAND (istrm) 

THETA  -  1. OdO /AVG 
NN-INT (10 . 0d0*AVG) 

GI-THETA 

FI-THETA 

1-1 
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IF (Uu.LE.FI)  GOTO  20 
DO  10  1-2, NN 
GI  -  (1. -THETA) *GI 
FI  -  FI  +  GI 
IF (Uu.LE.FI)  GOTO  20 
10  CONTINUE 

c  WRITE (*,21)  TNOW, AVG, Uu 

c21  FORMAT (IX, ' ! ! ! !  ERROR  IN  %  IGEOM  %,  TIME  -  ',F12.4, 
c  */2X, '  AVG  — ' , F12 .  4,  '  Uu  -  ',F8.5) 

20  IGEOM  -  I 
RETURN 
END 
C 

c  This  function  will  take  the  value  1  with  proba.  p  and  0  with  prob.  1-p 
function  nb (seed, pi) 
real*8  t, pi, rand 
integer  seed 
t-rand(seed) 
if (t. le.pl)  then 
nb-1 
else 
nb-0 
endif 
return 
end 

c 

function  ni j (seed,nc, i, q, nl,ml) 
implicit  real*8 (a-h, o-z) 
dimension  q(0:nl,0:ml) 
integer  seed 
t-rand(seed) 
sum-0 . 0 
do  20  j-0,nc 
sum-sum+q(i, j) 
if (  t.le.sum  )  then 
nij-j 
goto  15 
endif 

20  continue 

15  return 

end 

c 

c 

subroutine  vi (a,q,n,nl,ml,delt,v) 
c  this  subroutine  finds  the  steady  probability  v(i) 
implicit  real*8 (a-h, o-z) 

dimension  a(0:nl,0:ml) , q(0 : nl, 0 :ml) , v (0 : * ) 
c 

npl-n+1 
Do  10  i-0,npl 
10  a(0,i)-l.d0 

do  15  i-1, n 
do  15  j-0,n 

15  a(i,  j)—  q(j,i) 
do  16  i-l,n 

a(i,i)-l.d0+a(i,i) 

16  a  (i,npl)-0.d0 
c 

call  gauss(a,v,n,npl,nl,ml,delt) 

return 

end 

subroutine  qmat (q, x,m, n, aim, uc, nl,ml ) 
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implicit  real*8 (a-h,o-z) 
dimension  q(0 :nl, 0 :ml) , x (0 : *) 
c 

do  10  k-0,n 

10  x  ()c)-  (m*alm)  **k*dexp(-m*alm)  /ifac(k) 
c 

nml-n-1 
do  12  i-0,n 
do  11  j-0, nml 

11  q(i,j)-0.d0 
q(i,n)-l.d0 

12  continue 
c 

do  15  i-0,  n 
do  18  j-C'.nml 
km-max  <0, i- j) 
do  17  k-km, i 

17  q(i, j)-q(i, j) +cnk (i,k) *uc**k* (1-uc) ** (i-k) *x( j- (i-k) ) 

16  continue 

15  continue 

c 

do  20  i-0,n 
do  25  k-0, nml 

25  q(i,  n)  -q(i,n)  -q(i,  k) 

20  continue 

return 
end 
c 

function  ifac(k) 
if(k.eq.O)  then 
ifac-l 
else 
i-1 

do  10  j-l,k 
10  i-i*j 

ifac-i 
endif 
return 
end 
c 

function  cnk(n,k) 
implicit  real*8 (a-h, o-z) 
cnk-(l.d0*ifac(n) ) / (ifac(k) *ifac(n-k) ) 
return 
end 
C 

SUBROUTINE  GAUSS (A, X, N, NP 1 , N1 , Ml , DELT) 

THIS  ALGORITHM  SOLVES  THE  LINEAR  SYSTEM  EQUATION 
USING  THE  GAUSS -ELIMINATION  METHOD  WITH  PARTIAL  PIVOTING 
IT  ALSO  FINDS  THE  DETERMINENT  OF  A 
DOUBLE  PRECISION  IS  USED  IN  HTIS  SUBROUTINE 
SUBROUTINE  PVT  IS  USED  TO  INTERCHANGE  ROWS  IF  NECESSARY 
SUBROUTINE  UX  IS  CALLED  TO  FIND  THE  SOLUTION 
IF  THE  ERROR  IS  NEEDED,  THEN  CALL  SUBROUTINE  ERROR 

IMPLICIT  REAL*8  (A-H, O-Z) 

DIMENSION  A (0 :nl, 0 :ml) , X (0 : *) 

C 

NM1-N-1 
DO  15  I-0,NM1 
IP1-I+1 

CALL  PVTGUS (A,N,NP1, IP, I,N1,M1) 

C 

IF  (  DABS (A (I, I))  .LT.  DELT  )  GOTO  20 
C 

DO  61  K-IP1,N 
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C 


R1-A(K,I)/A(I,I) 

DO  62  J-IP1,NP1 
62  A (K, J) -A (K, J) -R1*A(I, J) 

61  CONTINUE 

15  CONTINUE 

IF  (  DABS (A (N, N) )  .LT.  DELT  )  GOTO  20 
CALL  UX ( A, N , NP 1 , X , N1 , Ml ) 

GO  TO  103 
C 

20  WRITE  (6,  *)  '  THE  MATRIX  IS  SINGULAR  ' 
103  RETURN 
END 


SUBROUTINE  PVTGUS ( A, N , NP1 , IP , J, N1 , Ml ) 
IMPLICIT  REAL *8  (A-H,0-Z) 

DIMENSION  A(0:nl,0:ml) 

IP-J 

BIG-DABS (A (J,  J)  ) 

DO  52  I— J, N 

IF  (  BIG  .LT.  DABS (A (I, J) )  )  THEN 
BIG-DABS (A (I , J) ) 

IP-I 
END  IF 
52  CONTINUE 

IF  (IP.EQ.J)  RETURN 
DO  61  K— 0 , NP1 
SAV— A ( J, K) 

A  ( J,  K)  —A  ( IP,  K) 

A  (IP,  K) -SAV 
61  CONTINUE 
RETURN 
END 


SUBROUTINE  UX(A,N,NP1,X,N1,M1) 

IMPLICIT  REAL* 8  (A-H,0-Z) 

THIS  SUBROUTINE  USES  BACKWARD  SUBSTITUTION  TO  FIND 
THE  SOLUTION  OF  LINEAR  SYSTEM  EQUATION 
A*X— B,  USING  THE  MATRIX  COMPUTED  BY  GUASS  OR  GAUSSP 
OR  GAUSSM 

DIMENSION  A(0 :nl, 0 :ml) , X(0 : *) 

X (N)-A(N,NP1) /A (N,N) 

DO  70  1-1, N 
IN-N-I 
SUM-0 . 0 
IN1-IN+1 
DO  75  J-IN1,N 

75  SUM-SUM+A (IN, J) *X ( J) 

X  (IN)  -  (A  ( IN,  NP1 )  -SUM)  /A  ( IN,  IN) 

70  CONTINUE 
C 

RETURN 

END 
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Appendix  G. 


The  FDDI  Program:  Symmetric  Case  Example 


fddi.l.inp 

**  Simulation  of  a  Timed  Token  Rotation  Protocol  (FDDI-I  Type) 
program:  FDD I 

Professor  Izhak  Rubin,  UCLA 

Enter  Select  Input  Mode 

1.  from  KEYBOARD 

2.  from  DATA  FILE 

1 

Enter  the  output  data  file  name 

ttl  .2 

Enter  Feature  Selection 

1.  Symmetric  System 

2.  2  Classes  of  Stations 

3.  Different  Loading  from  Station  to  Station 

1 

Enter  the  statistics  collection  start  time  (msec) 

4000 

Enter  the  stop  time  (msec) 

10000 

Enter  w,d,x,  for  computing  P (W>w) ,  P (D>d) ,  P (X>x) 

111 

Enter  the  number  of  stations  (N) 

6 

Enter  the  walk  time  from  a  station  to  its  neighboring  station  (in  msec;  r) 

0.1 

Enter  Target  Token  Rotation  Time  (msec;  TTRT) 

200 


Synchronous  traffic  arrival  rate  (packets/msec/station;  as) 
0.001 

Mean  synchronous  packet  transmission  time  (msec;  plena) 
0.05 

Enter  the  bandwidth  time  (msec;  <  (TTRT-N*r) /N) 

(the  max  synchronous  message  transmission  per  visit) 

0. 

Enter  number  of  message  priority  classes  per  station  (Np) 

3 

Enter  priority-1  threshold  (msec;  T_pri(l)) 

100 

Enter  priority-1  arrival  rate  (aa;  packets/msec/station) 

0.1 

Enter  priority-1  mean  transmission  time  (msec;  plena) 

0.3 

Enter  priority-2  threshold  (msec;  T_pri(2)) 

76.5 

Enter  priority-2  arrival  rate  (aa;  packets/msec/station) 

0.1 

Enter  priority-2  mean  transmission  time  (msec;  plena) 

0.3 

Enter  priority-3  threshold  (msec;  T__pri(3)) 

56.2 

Enter  priority-3  arrival  rate  (aa;  packets/msec/station) 

0.1 

Enter  priority-3  mean  transmission  time  (msec;  plena) 

0.3 
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This  is  an  example  of  the  input  file  for  the  symmetric  system. 


ff7 .2 

output  file 

1 

ifeat 

4000 

cl 

10000 

c2 

111 

w  d  x 

6 

N 

0.1 

r 

200 

TTRT 

0.001 

as 

0.05 

plens 

0. 

BWT 

3 

Np 

100. 

T_pri<l) 

0.1 

aa  (1) 

0.3 

plena (1) 

76.5 

T_pri<2) 

0.1 

aa  (2) 

0.3 

plena (2) 

56.2 

T_pri (3) 

0.1 

aa  (3) 

0.3 

plena (3) 

•  P*rfor«anca  of  a  TlMd  Tok#n  Rotation  Protocol  (FDDI-typa)  Ring  Networka  *•  Station  E<V)>  *lgraa(V))  Pr(W>  1.000)  Pr  (D>  1.000)  Pr  (X>  1) 
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Delay-Throughput  Evaluator,  IRI  Coip. 


Appendix  H. 

The  FDDI  Program:  An  Example  for  the  2  Classes  of 

Stations  Case 
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fddi.2.inp 

**  Simulation  of  a  Timed  Token  Rotation  Protocol  (FDDI-1  Type) 
program:  FDD I 

Professor  Izhak  Rubin,  UCLA 

Enter  Select  Input  Mode 

1.  from  KEYBOARD 

2.  from  DATA  FILE 

1 

Enter  the  output  data  file  name 
f  f  8 . 1 

Enter  Feature  Selection 

1.  Symmetric  System 

2.  2  Classes  of  Stations 

3.  Different  Loading  from  Station  to  Station 

2 

Enter  the  statistics  collection  start  time  (msec) 

400 

Enter  the  stop  time  (msec) 

2000 

Enter  w,d,x,  for  computing  P (W>w) ,  P (D>d) ,  P (X>x) 

10  10  2 

Enter  the  number  of  class-1  stations  (Nl) 

3 

Enter  the  number  of  class-1  stations  (N2) 

3 

Enter  the  walk  time  from  a  station  to  its  neighboring  station  (in  msec;  r) 

1 

Enter  number  of  message  priority  classes  per  station  (Np) 

3 

Enter  Target  Token  Rotation  Time  (msec;  TTRT) 

200 

Synchronous  traffic  arrival  rate  for  class-1  and  class-2  stations 
(packets/msec/station)  (asl,  as2) 

0.1  1 

Mean  synchronous  packet  transmission  time  for  class-1  and  class-2  stations 
(msec)  (plensl,  plens2) 

0.1  0.1 

Enter  the  Bandwidth  Time  for  class-1  and  class-2  stations  (msec) 

(BWT1*N1+BWT2*N2  <  (TTRT-walk  time) 

(the  max  synchronous  message  transmission  per  visit) 

20  40 

Asynchronous  Traffic 

Enter  priority-1  threshold  for  both  classes  (msec;  T_pri (classl, 1)  T_pri (class2, 1) ) 
40  75 

Enter  priority-1  arrival  rate  for  both  classes  (packets/msec/station)  (aal  aa2) 

0.01  0.1 

Enter  priority-1  mean  transmission  time  (msec)  (plenal  plena2) 

0.001  0.001 

Enter  priority-2  threshold  for  both  classes  (msec;  T_pri (classl, 2)  T_pri (class2, 2) ) 
20  30 

Enter  priority-2  arrival  rate  for  both  classes  (packets/msec/station)  (aal  aa2) 

0.01  0.1 

Enter  priority-2  mean  transmission  time  (msec)  (plenal  plena2) 

0.001  0.001 

Enter  priority-3  threshold  for  both  classes  (msec;  T_pri (classl, 3)  T_pri (class2, 3) ) 
10  16 

Enter  priority-3  arrival  rate  for  both  classes  (packets/msec/station)  (aal  aa2) 

n  n  i  n  1 

Enter  priority-3  mean  transmission  time  (msec)  (plenal  plena2) 

0.001  0.001 
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I 


fddi.inp2 


This  is  an  example  of  the  input  file  for  2  classes  of  stations. 


f  f  8 . 1 
2 

400 

2000 

111 

3 

3 

1 

3 

200 
0.1  1 
0.1  0.1 
20  40 
40  75 
0.01  0.1 
0.001  0.001 
20  30 
0.01  0.1 
0.001  0.001 
10  16 
0.01  0.1 
0.001  0.001 


output  file 

ifeat 

cstart 

cend 

w  d  x 

N1 

N2 

r 

Np 

TTRT 

asl,  as2 
plensl,  plens2 
BWT1,  BWT2 

T_pri (1,1)  T_pri (2,1) 
aa(l,l)  aa (2, 1) 
plena ( 1 , 1 )  plena (2,1) 
T_pri(l,2)  T_pri (2,2) 
aa (1/ 2)  aa (2,2) 
plena  (1,2)  plena (2, 2) 
T_pri (1,3)  T_pri (2,3) 
aa(l,3)  aa(2,3) 
plena (1,3)  plena (2, 3) 
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Delay-Throughput  Evaluator,  IRI  Corp. 


Appendix  I. 

The  FDDI  Program:  Asymmetric  Case  Example 
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fddL3.inp 


**  Simulation  of  a  Timed  Token  Rotation  Protocol 
program:  FDD I 

Professor  Izhak  Rubin,  UCLA 


(FDDI-I  Type) 


Enter  Select  Input  Mode 

1.  from  KEYBOARD 

2.  from  DATA  FILE 

1 

Enter  the  output  data  file  name 
f ddi . 3 . out 

Enter  Feature  Selection 

1.  Symmetric  System 

2.  2  Classes  of  Stations 

3.  Different  Loading  from  Station  to  Station 
3 

Enter  the  statistics  collection  start  time  (msec) 

400 

Enter  the  stop  time  (msec) 

2000 

Enter  w,d,x,  for  computing  P (W>w) ,  P (D>d) ,  P (X>x) 

111 

Enter  the  number  of  stations  (N) 

€ 

Enter  walk  times  from  station  to  station  (msec)  (r  (i) ,  i-1 .  .N) 

0.1  0.1  0.1  0.1  0.1  0.1 

Enter  Target  Token  Rotation  Time  (msec;  TTRT) 

20 

Synchronous  traffic  arrival  rates  (packets/msec/station;  as (i) , i-1 . .N) 
0.02  0.2  0.02  0.2  0.02  0.2 

Mean  synchronous  packet  transmission  time  (msec;  plens (i) , i-1 . .N) 

0.5  0.5  0.5  0.5  0.5  0.5 

Enter  the  bandwidth  times  (msec;  <TTRT-walktime)  (BWT (i) , i-1 . -N) 

3  3  3  3  3  3 


Asynchronous  Traffic 

Enter  number  of  message  priority  classes  per  station  (Np) 

3 

Enter  priority-1  threshold  for  stations  1..N  (msec) 

(T _pri  (i,  1) ,  i-1.  .N) 

10  10  10  10  10  10 

Enter  priority-1  arrival  rate  for  stations  1,..N  (packets/msec/station) 
(aid, 1), i-1.  .N) 

0.1  0.1  0.1  0.1  0.1  0.1 

Enter  priority-1  mean  transmission  time  for  stations  1..N  (msec) 

(plena  (i,  1) ,  i-1 .  .N) 

0.3  0.3  0.3  0.3  0.3  0.3 

Enter  priority-2  threshold  for  stations  1..N  (msec) 

(T__pri  (i,  2)  ,  i-1 .  .N) 

7.65  7.65  7.65  7.65  7.65  7.65 

Enter  priority-2  arrival  rate  for  stations  1, . .N  (packets/msec/station) 
(aa (i,2) , i-1 . .N) 

0.1  0.1  0.1  0.1  0.1  0.1 

Enter  priority-2  mean  transmission  time  for  stations  1..N  (msec) 

(plena (i, 2) , i-1. .N) 

0.3  0.3  0.3  0.3  0.3  0.3 

Enter  priority-3  threshold  for  stations  1..N  (msec) 

(T _pri(i,3),i-l..N) 

5.62  5.62  5.62  5.62  5.62  5.62 

Enter  priority-3  arrival  rate  for  stations  1,..N  (packets/msec/station) 
(aa (i, 3) , i-1 . .N) 

0.1  0.1  0.1  0.1  0.1  0.1 

Enter  priority-3  mean  transmission  time  for  stations  1..N  (msec) 

(plena  (i,  3) ,  i-1.  .N) 

0.3  0.3  0.3  0.3  0.3  0.3 
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fddi.inp3 


This  is  an  example  of  the  input  file 
to  station. 


f f 9 . 31 
3 

400 

2000 

111 

6 

3 

200 

0.1  30 
0.1 
0.5 

100.  76.5  56.2 

0.1  0.10.1 

0.3  0.30.3 

0.1  30 

0.1 

0.5 

100.  76.5  56.2 
0.1  0.10.1 
0.3  0.30.3 

0.1  30 
0.1 
0.5 

100.  76.5  56.2 

0.1  0.10.1 

0.3  0.30.3 

0.1  30 

0.1 

0.5 

100.  76.5  56.2 
0.1  0.10.1 
0.3  0.30.3 

0.1  30 
0.1 
0.5 

100.  76.5  56.2 
0.1  0.10.1 
0.3  0.3  0.3 

0.1  30 
0.1 
0.5 

100.  76.5  56.2 
0.1  0.10.1 
0.3  0.30.3 


for  different  loading  from  station 


output  file 

ifeat 

cl 

c2 

W0, DO, X0 

N 

Np 

TTRT 

r,  BWT 

as 

plens 
T_pri (1) 
aa  (1) 
plena (1) 
r,  BWT 
as 

plens 
T_pri (2) 
aa  (1) 
plena (2) 
r,  BWT 
as 

plens 
T_pri (3) 
aa  (3) 
plena (3) 
r,  BWT 
as 

plens 
Tjpri (4) 
aa  (4) 
plena (4) 
r,  BWT 
as 

plens 
T_pri(5) 
aa  (5) 
plena (5) 
r,  BWT 
as 

plens 
Tjpri ( 6 ) 
aa  (6) 
plena (6) 
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if  (i2.lt. m2)  goto  2 
if  (kl.lt.Nl)  goto  4 
if  (k2.1t.N2)  goto  5 
return 
end 

subroutine  blocks (nq, atail, t2) 
implicit  real*8  (A-H,0-Z) 
logical  ofs(20),zf 

real*8  tofs (20) ,cs (20) ,ys (20) ,wofs (20) 
integer  os (20) , ss (20) , seeds (20) 

common  /bs/ofs, ss, tofs, cs, os, seeds,  ys,  wof s  /bc/sp,B,zf 
temp  *=  atail 

1  si  =  dmod (sp*ss (nq) ,B) 

pi  -  -dlog (sl/B) *cs (nq) 
temp  -  temp+max (0.000001, pi) 
ss (nq)  «  si 
if  (temp. It. t2)  then 

if  (zf)  os (nq)  -  os (nq) +1 
goto  1 
end  if 

tofs (nq)  -  temp 

s2  -  dmod(sp*seeds (nq) ,B) 

wof s (nq)  -  max(0 . 000001, s2/B*ys (nq) ) 

seeds (nq)  -  s2 

ofs(nq)  =  .true. 

return 

end 

subroutine  blocka (nq, j, atail, t2) 
implicit  real*8  (A-H,0-Z) 
logical  ofa(20,3),zf 

real*8  tofa (20, 3) , ca (20, 3) ,ya (20, 3) , wof a (20, 3) 
integer  oa (20,3) ,sa (20,3) ,seeda (20,3) 

common  /ba/ofa, sa, tofa, ca, oa, seeda,ya, wof a  /bc/sp,B,zf 
temp  -  atail 

1  si  -  dmod(sp*sa (nq, j) ,B) 

pi  -  -dlog(sl/B) *ca (nq,  j) 
temp  -  temp+max (0.000001, pi) 
sa (nq, j)  -  si 
if  (temp. It. t2)  then 

if  (zf)  oa(nq, j)  -  oa (nq, j)+l 
goto  1 
end  if 

tofa(nq, j)  «  temp 

s2  -  dmod (sp*seeda (nq, j) , B) 

wofa(nq,j)  -  max (0 . 000001, s2/B*ya (nq, j) ) 

seeda(nq, j)  -  s2 

ofa(nq, j)  -  .true. 

return 

end 
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C> i tdma2 

Enter  the  file_name  for  data 
(no  more  than  12  characters) 

(in  PC,  QUOTE  the  'filename') 

' i out  1  ’ 

What  would  you  like  to  do  ? 

1:  Simulation  only 
2 :  Ana  lysis  only 
3:  Simulation  and  analysis 
Please  enter  a  number: 

3 

Enter  the  start  time  (when  statistic  starts) 

100 

Enter  the  stoo  time  (total  simulation  time) 

1000 

Enter  the  frame  duration  (m)  (m<  =  50)  and  the  max-imum 
number  of  slots  for  CS  support  (n) 

(n<=m)  m  and  n  should  be  integers 
10  5 

Enter  the  session  (circuit)  arrival  rate 

0.01 

Enter  the  session  transmission  rate 
0 . 04 

Enter  the  packet  batch  arrival  rate  [mess. /slot]  pi 
0  .  1 

please  wait 

average  number  of  CS  slots  used  per  frame  (V): 

Enter  the  batch  size  distribution  index 
1:  deterministic 
2 :  geometric 
3:  uniform 

2 

Enter  the  mean  batch  size,  b  ( m*b*p 1 <m-V ) : 

5 

To  find  the  probabilities  P(X>x)  and  P(D>d),  please 
enter  x  and  d  (x,  d<=  20,  both  should  be  integers) 

10  10 
please  wait 


.  334677242395972 
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******************************************************************************* 

* 

*  Simulation  of  Timed  Token  Protocol  MAC  (FDDI  Type) 


Version  1  Sep.  20,  1989 

Professor  Izhak  Rubin 

Department  of  Electrical  Engineering 

University  of  California 

Los  Angeles,  CA  90024 


Description  of  Input  Parameters 


TTRT 

N 

Np 

as  (i) 
aa  (i, j) 
plens 
plena 
Tpri (i, 
BWT 
r  (i) 
w,d,x 


Target  Token  Rotation  Time  (msec) 
number  of  stations 

number  of  priorities  in  each  station 

arrival  rate  for  synchronous  traffic  at  station  i 

arrival  rate  for  asynchronous  traffic  at  station  i,  priority  j 

mean  packet  length  for  synchronous  traffic 

mean  packet  length  for  asynchronous  traffic 

T-Pri,  priority  threshold  for  asynchronous  traffic 

synchronous  bandwidth  time  assignment 

walk  time  from  station  i  to  station  i+1 

specified  values  for  P (W>w) ,  P (D>d) ,  P(X>x) 


program  FDDI 

implicit  real*8  (A-H,0-Z) 

character  f i*12, fo*12, ch*l 

logical  gf , zf , sf , af , ofs (20) ,ofa(20,3) 

real  azs (20) ,dzs (20) ,aza (20,3) ,dza (20, 3) , tzs (20)  ,tza (20,3) 

*  , tds (20) , tws (20) , tda (20,3), twa (20,3) 

real* 8  Aas (20,0:200) , works (20, 0 : 200) , as (20) , plens (20) , 

*  Tpri (20, 3) ,aa (20, 3), plena (20, 3) ,Us (20) ,Xav(20,3) , 

*  Aaa (20,3, 0 : 100) , works (20, 3, 0 : 100) ,Ua (20, 3) , r (20) , res (20) , 

*  TRT,  tO (20) , BWT (20) ,v(20) , Wa (20, 3) , Ws (20) ,Ds(20) ,Da(20,3) 

*  , cs (20) , ca (20, 3) , ys (20) ,ya (20, 3) ,Xsm(20) ,Xsv(20) ,Xam(20, 3) 

*  ,sigws (20) ,sigds (20) ,sigwa (20,3) ,sigda (20,3) ,tya (20,3) 

*,tss (20,  0 : 50) , tsa (20,3,0:50) ,t jumps (20) ,t jumpa (20,3) ,tys (20) 

*, dwell (20) ,dwellv(20) ,ctime (20) ,ctimev(20) ,dpar (20) ,dparv(20) 

*, tqs (20) ,tqa (20,3), tqsv (20) , tqav (20,3) ,  tua (3) 

*, asc (2) ,plensc (2) ,BWTc (2) ,Tpric (2,8) ,aac (2,8) ,plenac (2,8) 

*, tofs (20) , tofa (20,3) , wofs (20) , wofa (20,3) 
integer  tails (20) , taila (20,3),ss(20),sa(20,3), seeds (20) , 

*seeda (20,3) , heads (20) ,heada (20,3) , txs (200, 0 : 50) ,txa (20,3,0:50), 
*zcs (20) , zca (20, 3), os (20), oa (20, 3), dxa (20,3,0:50) ,NN(2) 
*,Xs(20),Xa(20,3), axs (20, 0 : 50) ,dxs (20,0:50) ,axa (20,3,0:50) 

*,pds (20,0 : 50) ,pws (20, 0 : 50) ,pda (20, 3, 0 :50) ,pwa (20, 3, 0 : 50) , a (20) 
*,PrXs (20)  ,PrXa (20,3) ,PrWs (20) , PrWa (20,  3) ,PrDs (20) ,PrDa (20,3)  ,X0 
common  Aas, Aaa, works, worka, Us, Ua, heads, heads, tails, taila, N,Np, 

*  Xs , Xa , axs , dxs , axa , mind, tss,tsa, t jumps , t j umpa 

common  /syn/ v, res , Ws , Ds , s igws , s igds , scale , pws , pds , zcs , P rWs , PrDs , 

*  W0, DO  /bs/ofs, ss, tofs, cs, os, seeds, ys, wofs  /bc/sp,B,zf 
*/ba/ofa, sa, tofa, ca, oa, seeda, ya,  wofa 

m<nd  ■  30 

base  for  random  numbers 
sp  ■  dble(7**5) 
yz  -  1 . 

do  101  i  -  1,31 
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101 

*** 


211 


212 


yz  -  yz*2 
B  -  yz-1 

Input  Section  *** 

write <*,*)  ' **  Simulation  of  a  Timed  Token  Rotation  Protocol' 

r  , '  (FDDI-I  Type)' 

write (*,*)  'program:  FDD I' 

write (*,*)  'Professor  Izhak  Rubin,  UCLA' 

write (*, *) 

write (*,*)  'Enter  Select  Input  Mode' 
write (*,*)  '1.  from  KEYBOARD' 
write  (*,*)  '2.  from  DATA  FILE' 
read(*,*)  mode 
write (*,*)  mode 
if  (mode.eq.2)  then 

write (*,*)  'state  name  of  data  file  (<«12  characters)' 

read  (*,*)  fi 

open(7,file*fi) 

read  (7,*)  fo 

open(3,file-fo) 

read (7,*)  ifeat 

read (7,*)  cl 

read (7,*)  c2 

read (7, *)  W0,D0,X0 

if  (ifeat. eq.l)  then 
read (7,*)  N 
read(7, *)  r(l) 
read (7,*)  TTRT 
read(7, *)  as(l) 
read(7,*)  plens(l) 
rhos  *  N*as (1) *plens (1) 
if  (rhos-ge.l.)  then 

write (*,*)  'The  synchronous  traffic  normalized  throughput' 
write (*,*)  'Rhos  -  N*as*plens  >»  1' 
end  if 

read (7,  *)  BWT(l) 

if  (BWT(l) .gt. (TTRT-N*r(l) )/N)  then 
write (*,*)  ' BWT  >  (TTRT-N*r) /N' 
end  if 
read (7,*)  Np 
rho  rhos 
do  211  j  -  1, Np 

read (7,*)  Tpri(l,j) 
read (7,*)  aa(l,j) 
read(7,*)  plena (l,j) 
rho  -  rho+aa (1, j) ‘plena (1, j) *N 
if  (rho.ge.l)  write (*,508)  j,rho,j 
continue 
do  212  i  -  2,N 
r(i)  -  r (1) 
as (i)  -  as  (1) 
plena (i)  -  plena (1) 

BWT  (i)  -  BWT (1) 
do  212  j  -  1, Np 

Tpri (i, j)  -  Tpri (1, j) 
aa(i,j)  -  aa (1, j) 
plena (i,j)  -  plena (l,j) 

else  if  (ifeat. eq. 2)  then 
read (7,*)  N1 
read (7,*)  N2 
N  -  N1+N2 

call  pmt (N1,N2,N, a) 
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* 


531 


read (7,*)  rc 

rr  -  rc*N 

read (7,*)  Np 

read (7,*)  TTRT 

read(7,*)  asc (1) , asc (2) 

read(7,‘)  plensc (1) , plensc (2) 

rho  -  asc  (1)  ’‘plensc  (1)  *Nl+asc  (2)  ‘plensc  (2)  ‘N2 

If  (rho.ge.l.)  then 

write (*,*)  'The  synchronous  traffic  normalized  throughput' 
write  (*/*)  'Rhos  -  asl*plensl+as2*plens2  >*»  1' 
end  if 

read(7, *)  BWTc (1) ,BWTc (2) 

if  ( (BWTc (1) *Nl+BWTc (2) *N2+rr) .gt .TTRT) 

write (*,*)  'Synchronous  BWT  >  (TTRT-walk  time)' 
do  531  j  *  l,Np 

read(7,‘)  Tpric (1, j) ,Tpric (2, j) 
read(7,‘)  aac (1, j) ,aac (2, j) 
read(7, *)  plenac (1, j) ,plenac (2, j) 

rho  ■  rho+aac (1, j) ‘plenac (1, j) *Nl+aac (2, j) ‘plenac (2, j) *N2 
if  (rho.ge.l)  write(*,508)  j,rho,j 
continue 


NN  (1)  -  N1 
NN  (2)  -  N2 
do  533  i  -  1,N 
r (i)  «  rc 
as  (i)  -  asc  (a  (i) ) 
plens(i)  -  plensc (a (i) ) 

BWT (i)  -  BWTc (a (i) ) 
do  533  j  *  1, Np 

Tpri (i, j)  -  Tpric (a (i) , j) 
aa(i,  j)  -  aac(a(i)  ,  j) 

533  plena  (i,j)  «•  plenac  (a  <i) ,  j) 


else 

read ( 7 , * )  N 
read (7,*)  Np 
read (7,*)  TTRT 
rho  -  0. 
temp  -  0 . 
do  221  i  -  1,N 

read(7, *)  r(i),BWT(i) 
read (7,*)  as(i) 
read (7,*)  plens(i) 
rho  -  rho+as (i) ‘plens (i) 
temp  -  temp+BWT (i) +r (i) 
read (7, ‘)  (Tpri (i, j) , j-l,Np) 
read (7, ‘)  (aa (i, j) , j-l,Np) 
read (7, *)  (plena (i, j) , j-l,Np) 
do  221  j  -  l,Np 

221  rho  ••  rho+aa (i, j) ‘plena  (i, j) 

if  (rho.ge.l.)  write (*,*)  'The  normalized  throughput  >«  1' 
if (temp. gt .TTRT)  write (*,*)' Synchronous  BWT  >  TTRT-walktime' 
end  if 
goto  99 
end  if 


write (*,*)  'Enter  the  output  data  file  name' 
read  (*,*)  fo 
~~?r.  (3,  f ile»fo) 


write (•,*)  'Enter  Feature  Selection' 
write (*,*)  '1.  Symmetric  System' 
write (*,*)  '2.  2  Classes  of  Stations 
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write (*,*)  '3.  Different  Loading  from  Station  to  Station' 
read(*,*)  ifeat 

write (*,*)  'Enter  the  statistics  collection  start  time  (msec)' 
read(*,*)  cl 

write (*,*)  'Enter  the  stop  time  (msec)' 
read(*,*)  c2 

write(*,*)  'Enter  w,d,x,  for  computing  P (W>w) ,  P (D>d) ,  P (X>x) ' 
read<*,*)  W0,D0,X0 


502 


601 

503 


if  (ifeat. eq.l)  then 

write (*,*)  'Enter  the  number  of  stations  (N) ' 
read(*,*)  N 
write (*, 502) 

format ('  Enter  the  walk  time  from  a  station  to  its  neighboring' 

*  , '  station  (in  msec;  r) ' ) 
read(*,*)  r(l) 

write (*,*)  'Enter  Target  Token  Rotation  Time  (msec;  TTRT)' 
read (*, *)  TTRT 
write (*, 503) 

format (/,'  Synchronous  traffic  arrival  rate', 

*  '  (packets/msec/station;  as)') 
read(*,*)  as(l) 

write (*,*)  'Mean  synchronous  packet  transmission  time', 

*  '  (msec;  plens)' 
read(*,*)  plens (1) 
rhos  -  N*as (1) *plens (1) 
if  (rhos.ge.l.)  then 

write (*,*)  'The  synchronous  traffic  normalized  throughput' 
write (*,*)  'Rhos  ”  N*as*plens  >-  1' 

write (*,*)  'Re-select  synchronous  loading  parameters' 
goto  601 
end  if 


602  write <*,*)  'Enter  the  bandwidth  time  (msec;  <  (TTRT-N*r) /N) ' 
write (*,*)  '(the  max  synchronous  message  transmission', 

*  '  per  visit) ' 
read(*,*)  BWT(l) 

if  (BWT(l) .gt. (TTRT-N*r(l))/N)  then 

write <*,*)  ' Re-select :BWT  >  (TTRT-N*r) /N' 
goto  602 
end  if 
write (*, 504) 

504  format (/,'  Enter  number  of  message  priority  classes', 

*  '  per  station  (Np)') 
read(*,*)  Np 

rho  ■  rhos 
do  401  j  -  1, Np 
write (*, 505)  j,j 


505 

format ('  Enter  priority-' , il, ' 

threshold' , 

* 

'  (msec;  T_pri (' , i, ' ) ) ' ) 
read ( * , * )  Tpri(l,j) 

603 

write (*,506)  j 

506 

format ('  Enter  priority-' ,  il, ' 

arrival  rate', 

* 

'  (aa;  packets/msec/station)') 
read<*,*)  aa(l,j) 
write (*,507)  j 

507 

format ('  Enter  priority-' , il. 

*  '  mean  transmission  time  (msec;  plena)') 
read(*,*)  plena (l,j) 

ro  -  rho+aa  (1,  j)  *plena  (1,  j)  *N 
if  (ro.ge.l)  then 
write (*,508)  j,ro,j 

508  format ('  Total  normalized  throughput  for  synchronous  and  ' 

*  ,/,'  asynchronous  traffic  of  priority  no  lower  than  ',il. 
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'  is' ,  f5 .2, / , '  Re-select  priority-' , il, '  loading  parameters') 
goto  €03 
end  if 
rho  -  ro 

401  continue 

do  402  i  -  2,N 

r  (i)  -  r  (1) 
as(i)  -  as(l) 
plena  (i)  -  plens(l) 

BWT  <i)  -  BWT(l) 
do  402  j  -  l,Np 

Tpri (i, j)  -  Tpri(l,j) 
aa  (i,  j)  -  aa (1, j) 

402  plena (i,j)  -  plena (l,j) 

else  if  (ifeat.eq.2)  then 
***  ifeat  -  2 


write (*,*)  'Enter  the  number  of  class-1  stations  (Nl)' 
read(*,*)  Nl 

write (*,*)  'Enter  the  number  of  class-1  stations  <N2) ' 
read(*,*)  N2 
N  -  N1+N2 

call  pmt (N1,N2,N, a) 
write (*,502) 
read(*,*)  rc 
rr  -  rc*N 
write (*,504) 
read(*,*)  Np 

write (*,*)  'Enter  Target  Token  Rotation  Time  (msec;  TTRT) ' 
read ( * , * )  TTRT 

811  write (*, 813) 

813  format  (/,'  Synchronous  traffic  arrival  rate  for  class-1  and  ', 

*  'class-2  stations',/,'  (packets/msec/station)  (asl,  as2)') 
read(*, *)  asc  (1) ,asc  (2) 

write (*, 815) 

815  format ('  Mean  synchronous  packet  transmission  time  for  ', 

*  '  class-1  and  class-2  stations',/,'  (msec)  (plensl,  plens2)') 
read(*,*)  plensc (1) , plensc (2) 

rhos  -  asc (1) *plensc (1) *Nl+asc (2) *plensc (2) *N2 
if  (rhos.ge.l.)  then 

write (*,*)  'The  synchronous  traffic  normalized  throughput' 
write  (*,*)  'Rhos  ”  asl*plensl+as2*plens2  >-  1' 
write (*,*)  'Re-select  synchronous  loading  parameters' 
goto  811 
end  if 

812  write(*,816) 

816  format ('  Enter  the  Bandwidth  Time  for  class-1  and  class-2  ', 

*  'stations  (msec)',/,'  (BWT1*N1+BWT2*N2  <  (TTRT-walk  time) ' , /, 

*  '  (the  max  synchronous  message  transmission  per  visit)') 
read(*, *)  BWTc (1) ,BWTc (2) 

if  ( (BWTc (1) *Nl+BWTc (2) *N2+rr) .gt .TTRT)  then 

write(*,*)  'Re-select;  synchronous  BWT  >  (TTRT-walk  time)' 
goto  812 
end  if 

write (*,*)  'Asynchronous  Traffic' 
rho  *  rhos 
do  431  j  -  1 , Np 
write (*, 805)  j,j,j 

805  format ('  Enter  priority-' ,  il, '  threshold  for  both  classes', 

*  '  (msec;  T_pri (classl, ' ,  il, ' )  T_pri (class2, ' ,  il, ' ) ) ' ) 
read(*,*)  Tpric (1, j) ,Tpric (2, j) 

803  write(*,806)  j 

806  format ('  Enter  priority-' , il, '  arrival  rate  for  both  classes' 
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, '  (packeta/maec/atation)  <aal  aa2)') 
read(*,*)  aac (1, j) ,aac (2, j) 
write (*,807)  j 

format ('  Enter  priority-' , il, 

'  mean  tranamission  time  (maec)  (plenal  plena2)') 
read(*, *)  plenac (1, j) ,plenac (2, j) 

ro  -  rho+aac <1, j) *plenac (1, j) *Nl+aac (2, j) *plenac (2, j) *N2 
if  (ro.ge.l)  then 
write (*,508)  j,ro,j 
goto  803 
end  if 
rho  -  ro 

NN (1)  -  N1 
NN (2)  -  N2 
do  433  i  -  1,N 
r(i)  -  rc 
aa (i)  -  aac (a  (i) ) 
plena  (i)  -  plenac  (a  (i) ) 

BWT(i)  -  BWTc(a(i) ) 
do  433  j  -  l,Np 

Tpri  (i,  j)  -  Tpric(ad) ,  j) 
aa (i, j)  -  aac  (a  (i) , j) 

433  plena (i,j)  -  plenac (a (i) , j) 

elae 

***  ifeat  -  3 


807 


431 


write(*,*)  'Enter  the  number  of  atationa  (JJ) ' 
read ( * ,  * )  N 

write (*,*)  'Enter  walk  timea  from  atation  to  atation', 

*  '  (maec)  (r (i) , i-1 . .N) ' 

read(*, *)  (r(i),i-l,N) 

rr  -  0. 

do  921  i  -  1 , N 
921  rr  -  rr+r(i) 

write (*,*)  'Enter  Target  Token  Rotation  Time  (maec;  TTRT) ' 
read(*,*)  TTRT 

933  write (*,*)  ' Synchronoua  traffic  arrival  ratea', 

*  '  (packeta/maec/station;  aa (i) , i-1 . .N) ' 
read(*, *)  (aa (i) , i«l,N) 

write (*, *)  'Mean  aynchronoua  packet  tranamiaaion  time', 

*  '  (maec;  plena (i) , i-1 . .N) ' 
read(*, *)  (plena (i) , i-1, N) 
rhoa  -  0. 

do  923  i  -  1, N 

923  rhoa  -  rhoa+aa (i) *plena (i) 

if  (rhoa.ge.l.)  then 

write (*,*)  'The  aynchronoua  traffic  normalized  throughput' 
write(*,*)  'Rhoa  -  aum  of  a (i) *plena (i)  >-  1' 
goto  933 
end  if 

932  temp  -  rr 

write (*,*)  'Enter  the  bandwidth  timea  (maec;  ', 

*  ' <TTRT-walktime)  (BWT (i) , i-1 . .N) ' 
read ( * , * )  (BWT (i) , i-l,N) 

do  924  i  -  1,N 

924  temp  -  temp+BWT(i) 

if  ( temp. gt .TTRT)  then 

write (*,*)' aynchronoua  BWT  >  (TTRT-walk  time)' 
goto  932 
end  if 

write (*,*)  ' Aaynchronoua  Traffic' 
write (*, 504) 
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read(*,*)  Np 
rho  -  rhos 
do  931  j  -  1 » Np 
write (*,915)  j,j 

915  format ('  Enter  priority-' , il, '  threshold  for  stations  1' 

*  ,'..N  (msec)',/,'  (T_pri (i, ' , il, ' ) , i-1 . .N) ' ) 
read(*, *)  (Tpri (i, j) , i-l,N) 

903  write (*, 906)  j,j 

906  format ('  Enter  priority-' , il, '  arrival  rate  for  stations', 

*  '  1,..N  (packets/msec/station) ' , /, '  (aa (i, ' , il, ' ) , i-1 . .N) ' ) 
read (*, *)  (aa (i, j) , i-l,N) 

942  write (*, 907)  j,j 

907  format ('  Enter  priority-' , il, '  mean  transmission  time  for', 

*  '  stations  1..N  (msec)',/,'  (plena (i, ', il, '), i-1 . .N) ' ) 
read (*, *)  (plena <i, j) , i-1, N) 

ro  -  rho 
do  941  i  -  1,N 

941  ro  -  ro+aa (i, j) *plena (i, j) 

if  (ro.ge.l)  then 
write(*,508)  j,ro,j 
goto  942 
end  if 

931  rho  -  ro 


end  if 

************************  End  of  Input  Section  ************************** 

99  do  403  i  -  1,N 

k  -  i*32 
ss(i)  -  k 
seeds (i)  -  k+1 
do  403  j  -  1, Np 
sa(i,j)  -  k+2+j 

403  seeda(i,j)  -  k+10+j 

ccc  -  c2/10. 
do  21  nq  -  1,N 
os(nq)  -  0 
cs (nq)  -  l./as(nq) 
ys (nq)  -  plena (nq) *2. 
sl  -  dmod(  sp*ss (nq) ,  B) 
s2  -  dmod(  sp*seeds (nq) ,  B) 
pl  -  -dlog (sl/B) *cs (nq) 

Aas(nq,0)  -  max (0 . 000001, pl) 
works(nq,0)  -  max (0 .000001, s2/B*ys (nq) ) 
ss(nq)  -  sl 
seeds (nq)  -  s2 
do  21  kp  -  1, Np 
oa(nq, kp)  -  0 
ca(nq, kp)  -  l./aa(nq, kp) 
ya(nq, kp)  -  plena (nq, kp) *2 . 
sl  -  dmod(  sp*sa (nq, kp) ,  B) 
s2  -  dmod(  sp*seeda (nq, kp) ,  B) 
pl  -  -dlog (sl/B) *ca (nq, kp) 

Aaa(nq,kp,0)  -  max (0 . 000001, pl) 

worka (nq, kp, 0)  -  max (0. 000001, s2/B*ya (nq, kp) ) 

aa(nq, kp)  -  sl 

21  seeda(nq,kp)  -  s2 

' ‘ 1  cycle  *** 

***  no  station  is  allowed  to  transmit  during  1st  cycle 
Rs  -  0 . 
zf  -  .false. 

TRTmax  -  0 . 
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do  5  i  ■  1,H 
tails (i)  -  0 
heads  (i)  -  0 
res(i)  -  works (i,0) 
zcs(i)  -  0 
Ds(i)  -  0. 
do  6  j  ■  l,Np 
zca(i,j)  -  0 
Da(i,j)  -  0. 
heads (i/j)  -  0 
tails (i,j)  -  0 
Ua(i,j)  -  0. 

6  continue 

Us (i)  -  0. 

Rs  “  Rs+r(i) 

5  continue 

if  (rho.lt. 0.999)  then 

scale  ■  Rs*rho/ (1-rho) *3 . / (mind+1 . ) 
else 

scale  -  TTRT*2.0/ (mind+1. 0) 
end  if 

***  2nd  cycle  *** 

***  only  synchronous  traffic  is  allowed  during  2nd  cycle 
t  -  Rs 

do  10  nq  -  1,N 
call  unfw(nq,  t) 
tO  (nq)  -  t 
TRT  «  t-tO (nq) 
if  (TRT.gt . TRTmax)  then 
TRTmax  “  TRT 
nmax  -  it 

tmax  -  t 
end  if 

v (nq)  -  0. 
y  -  0. 

gf  ■=  .false. 

8  if  (Us (nq) .eq.O.O.or.gf)  goto  10 

diff  -  BWT (nq) -y 
if  (Us (nq) .ge.diff )  then 
call  qwus (nq,t,diff ) 
gf  ■  .true, 
else 

temp  -  Us tnq) 
call  qwus (nq,t, temp) 
y  -  y+temp 
goto  8 
end  if 

10  t  -  t+r(nq) 


***  afterwards-  to  the  end  *** 

***  both  synchronous  and  asynchronous  messages  are  allowed 
***  assume:  exhausted  service, 

***  synchronous  message  is  served  first 

***  priority  1  is  the  highest  priority 

ck  -  ccc 
it  -  2 

30  if  (t.le.c2)  then 


if  (t.gt.c2)  goto  77 
it  -  it+1 
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if  (.not.zf)  then 
if  (t.ge.cl)  then 
nstart  -  it 
tatart  “  t 
zf  -  .true, 
end  if 
end  if 

***  diaplay  in  progreaa 

if  (t.gt.ck)  then 
ck  “  ck+ccc 

write (*,71)  t,t/c2*100. 

71  format('  [Simulation  Time] ' , f 12 . 5, '  msec' , 3x, f 6 .2, ' %' ) 

end  if 

do  31  nq  -  1,N 
TRT  «=  t-tO  (nq) 
if  (TRT .gt .TRTmax)  then 
TRTmax  =  TRT 
nmax  *  it 
tmax  -  t 
end  if 

vtemp  “  TRT-v(nq) 

call  unfw(nq,t) 
tO (nq)  -  t 

if  (zf)  then 

if  (Xs (nq) .gt -X0)  PrXs (nq)  -  PrXs(nq)+l 
if  (Xs (nq) .gt .mind)  then 

txs(nq,mind)  -  txs (nq,mind) +1 
else 

txs (nq,Xs (nq) )  -  txs (nq, Xs (nq) ) +1 
end  if 

do  38  i  -  l»Np 

if  (Xa (nq, i) .gt .XO)  PrXa(nq,i)  -  PrXa(nq,i)+l 
if  (Xa (nq, i) .gt .mind)  then 

txa (nq, i,mind)  ■  txa (nq, i,mind) +1 
else 

txa (nq, i, Xa (nq, i) )  -  txa (nq, i, Xa (nq, i) ) +1 
end  if 

38  continue 

end  if 

do  496  i  -  l,Np 

496  tua(i)  «  0. 

v(nq)  -  0. 
y  -  0. 

gf  -  .false. 

32  if  (Us (nq) . It . 0 . OOCOOOl .or . gf )  goto  34 

33  diff  -  BWT (nq) -y 

if  (diff .It. 0.0000001)  then 
gf  -  .true, 
goto  34 
end  if 

if  (Us (nq) .ge.diff )  then 
call  qwus (nq,t,diff ) 
gf  -  .true, 
else 

temp  -  Us (nq) 
call  qwus (nq,t, temp) 
y  «  y+temp 
goto  32 
end  if 
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do  35  j  -  l,Np 

if  ( (TRT+v (nq) ) . ge . Tpri (nq, j ) )  then 
if  (gf .or. Us (nq) .eq.O.)  goto  35 
goto  33 

else  if  (Ua (nq, j) .gt. 0.0000001)  then 
ww  -  worka (nq, j,heada (nq,  j) ) 
if  (zf)  tua(j)  ■*  tua(j)+ww 

t2  “  t+ww 
if  (ofs(nq))  then 

if  (tofs (nq) .It .t2)  call  blocks (nq, tofs (nq) , t2) 
goto  54 
end  if 

ttl  -  Aas (nq, tails (nq) ) 
if  (ttl.lt .t2)  then 

Us (nq)  -  Us (nq) +works (nq, tails (nq) ) 
k  -  mod(tails (nq) +1, 201) 
if  ( k.eq. heads (nq) )  then 
call  blocks (nq, ttl, t2) 
goto  54 
end  if 

if  (zf)  then 

if  (Xs (nq) .gt .mind)  then 

axs (nq,mind)  -  axs (nq,mind) +1 
tss (nq,mind) “tss (nq,mind) +ttl-t jumps (nq) 
else 

axs (nq,Xs (nq) )  “  axs (nq, Xs (nq) ) +1 
tss (nq, Xs (nq) )«tss (nq,Xs (nq) ) +ttl-t jumps (nq) 
end  if 
end  if 

Xs (nq)  -  Xs (nq) +1 
t jumps (nq)  -  ttl 
si  -  dmod(sp*ss (nq) ,B) 
pi  -  -dlog (sl/B) *cs (nq) 
ttl  -  ttl+max(0. 000001, pi) 
tails (nq)  -  k 
ss (nq)  -  si 
Aas (nq, k)  -  ttl 
s2  -  dmod(sp*seeds (nq) ,  B) 
works (nq, k)  -  max (0 . 000001, s2/B*ys (nq) ) 
seeds (nq)  -  s2 
goto  52 
end  if 

do  51  i  “  l,Np 

if  (ofa (nq, i) )  then 

if  (tofa (nq, i) . It .t2)  call  blocka (nq, i, tofa (nq, i) , t2) 
goto  51 
end  if 

ttl  ■  Aaa (nq, i,taila (nq, i) ) 
if  (ttl . It . t2)  then 

Ua(nq, i)  -  Ua (nq, i) +worka (nq, i, tails (nq, i) ) 
k  -  mod(taila(nq,i)+l,101) 
if  (k.eq.heada (nq, i) )  then 
call  blocka(nq,i,ttl,t2) 
goto  51 
end  if 
if  (zf)  then 

if  (Xa (nq, i) .gt .mind)  then 

axa (nq, i,mind)  -  axa (nq, i,mind) +1 
tsa (nq, i,mind) “tsa (nq, i,mind) +ttl-t jumpa (nq, i) 
else 

axa (nq, i,Xa (nq, i) )  -  axa (nq, i,Xa (nq, i) ) +1 
tsa (nq, i,Xa (nq, i) )  -  tsa (nq, i, Xa (nq, i) ) +ttl- 

t jumpa (nq, i) 
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end  if 
end  if 

Xa(nq, i)  -  Xa(nq,i)+1 
t jumpa (nq, i)  -  ttl 
si  -  dmod(sp*sa (nq, i)  ,B) 
s2  -  dmod(sp*seeda (nq, i)  ,B) 
pi  =  -dlog (sl/B) *ca (nq, i) 

Aaa(nq,i,k)  -  ttl+max (0 . 000001, pi) 
works  (nq,  i,  k)  m  max  (0 . 000001,  s2/B*ya  (nq,  i) ) 
sa(nq,i)  -  si 
seeda(nq,i)  =  s2 
taila (nq, i)  -  k 
goto  53 
end  if 
continue 


Xa (nq, j)  -Xa(nq, j ) — 1 
if  (ofa(nq, j))  then 
ofa(nq, j)  -  .false, 
k  <■  mod  (taila  (nq,  j)  +1, 101) 

Aaa (nq, j,k)  «  tof a (nq, j) 
works (nq, j,k)  -  wofa(nq, j) 
taila (nq, j)  -  k 
end  if 

if  (zf)  then 

temp  -  t-Aaa (nq, j, heada (nq,  j) ) 

if  (temp. gt. W0)  PrWa(nq, j)  -  PrWa (nq, j)+l 

if  ( (temp+ww) .gt .DO)  PrDa (nq, j)  ■  PrDa (nq, j)+l 

sigwa (nq, j)  -  sigwa (nq, j) +temp*temp 

sigda (nq, j)  -  sigda (nq, j )+ (temp+ww) * (temp+ww) 

kd  -  min ( (temp+ww) /scale, mind) 

kw  -  min (temp/scale, mind) 

pda (nq, j, kd)  -  pda (nq, j, kd) +1 

pwa (nq, j, kw)  -  pwa (nq, j, kw) +1 

Wa (nq, j)  -  Wa (nq, j ) +temp 

Da (nq, j)  ■  Da (nq, j ) +temp+ww 

if  (Xa (nq, j) .ge.mind)  then 

dxa  (nq,  j,mind)  ■>  dxa  (nq,  j,mind)  +1 
tsa (nq, j,mind)  -  tsa (nq, j, mind) +t2-t jumpa (nq, j) 
else 

dxa (nq, j,Xa (nq, j) )  -  dxa (nq, j,Xa (nq, j) ) +1 
tsa (nq, j,Xa (nq, j ) +1 )  -  tsa (nq, j,Xa (nq, j) +1) +t2- 
*  t jumpa (nq, j) 

end  if 

zca(nq, j)  -  zca (nq, j)+l 
end  if 

t jumpa (nq, j)  -  t2 
Ua(nq, j)  -  Ua(nq,j)-ww 
v(nq)  «  v(nq)+ww 
t  -  t2 

heada (nq,j)  -  mod (heada (nq, j) +1, 101) 
goto  32 
end  if 
35  continue 


The  token  leaves  station  nq,  takes  walk  time  r(nq) 
and  goes  to  station  (nq+1) 


if  (zf)  then 

tqs(nq)  -  tqs(nq)+y 
tqsv(nq)  -  tqsv(nq)+y*y 
do  499  i  -  l,Np 
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tqa(nq,i)  -  tqa (nq, i) +tua (i) 
tqav(nq,  i)  =  tqav (nq, i) +tua (i) *tua (i) 
dwell (nq)  —  dwell (nq) +v (nq) 
dwellv(nq)  -  dwellv (nq) +v (nq) *v (nq) 
ctime (nq)  -  ctime (nq) +vtemp+v (nq) 

ctimev(nq)  ■  ctimev (nq) + (vtemp+v (nq) ) * (vtemp+v (nq) ) 
dpar (nq)  -  dpar (nq) +vtemp 
dparv(nq)  -  dparv (nq) +vtemp*vtemp 
end  if 
t  -  t+r(nq) 
goto  30 
end  if 


***  Output  Section  *** 

77  write (3, 701) 

701  format (/,'  **  Performance  of  a  Timed  Token  Rotation' 
*,'  Protocol  (FDDI-type)  Ring  Networka  **') 

cyclen  -  it-nstart+1 
dtr  -  0. 
do  555  i  -  1,N 
555  dtr  -  dtr+dwell(i) 

thn  -  dtr/ (t-tstart) 
dtr  -  dtr/cyclen 
sf  -  .false, 
af  -  .false, 
do  577  i  -  1,N 

if  (BWT(i) .gt.O.)  sf  -  .true, 
do  577  j  -  l,Np 

if  (Tpri (i, j) .gt . 0 . )  af  -  .true. 

577  continue 


if  (ifeat.eq.l)  then 

write (3, 702)  'Symmetric  Systems  ' 

else  if  (ifeat.eq.2)  then 

write (3, 702)  '2  Classes  of  Stations  ' 

else 

write (3, 702)  'Different  Loading  from  Station  to  Station' 
end  if 

format (/,'  Feature  Selected:  ',a41) 

write  (3, 703)  cl, c2,TTRT,N,Np, 1 .O-Rs/2 . 0/TTRT, rho, thn, 

*  (t-tstart) /cyclen, dtr, TRTmax,nmax,tmax 
format ('  Statistics  Start  (msec) : ‘  , f 9 . 1, 

*  /,'  Statistics  Stop  (msec) . 1  £ -. 1, 


/,'  TTRT 


(msec)  : ' ,  £ 9 .3, 


/,'  Number  of  Stations  (N):',i9, 

/,'  Number  of  Priorities  (Np):',i9, 

/,'  Max  Throughput  (l-walktime/2TTRT) :' ,f9.4, 

/,'  Normalized  Throughput  (specif ied) : ' , f 9 . 4, 

/,'  Normalized  Throughput  (realized) :', f 9 . 4, 

/,'  Realized  Mean  Cycle  Time  (msec) : ' , f 9 . 4, 

/,'  Realized  Mean  Dwell  Time  (msec) : ' ,  f 9. 4, 

/,'  Max  Cycle  Time  (msec) :' ,f 9.4, 

*  '  at',i6,'th  cycle,  t-',f9.1) 

if  (ifeat.eq.l)  then 
write (3,704)  r(l) 

if  (sf)  write (3,758)  BWT(l) ,as (1) ,plens (1) 
if  (af)  then 

write (3,759)  'i','i' 

write (3, 705)  (i,Tpri (1, i) , aa (1, i) , plena (1, i) ,  i-l,Np) 
end  if 

else  if  (ifeat.eq.2)  then 
do  778  j  -  1,2 
write (3, 776)  j 
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ch  -  char(48+j) 
write(3,704)  r(j) 

if  (sf)  write(3,758)  BWT ( j ) , as ( j ) , plena ( j ) 
if  (af)  then 

write(3,759)  ch,ch 

write  (3,705)  (i,Tpri ( j,i) ,aa ( j,i) , plena ( j, i) , i-l,Np) 
end  if 
778  continue 

write (3, 747) 
k  -  mod (N, 20) 
if  (mod(N,20) .eq.O)  then 
k  -  N/20 
else 

k  -  N/20+1 
end  if 

do  742  i  -  l,k 
j  «  min(i*20,N) 

write  (3,749)  (ii, ii- (20* (i-1) +1) , j) 

742  write (3,746)  (a (ii) , ii- <20* (i-1) +1) , j) 

747  format (/,'  System  Configuration') 

749  format  (/,'  station' , 20i4) 

746  f ormat (  '  class:  ',20i4) 

else 

do  777  j  -  1,N 
write (3, 776)  j 
ch  ■  char(48+j) 
write(3,704)  r(j) 

if  (sf)  write(3,758)  BWT( j) , as (j) , plena (j) 
if  (af)  then 

write(3,759)  ch,ch 

write (3,705)  (i,Tpri ( j,i) ,aa ( j, i) , plena ( j,i) ,i»l,Np) 
end  if 
777  continue 
776  format (/,'  Station', i3) 
end  if 

704  format ('  Walk  Time  (r;msec) : ' , f 10 . 4) 

758  format('  Bandwidth  Time  (BWT;msec) : ' , f 10 . 4, 

*  /,'  Arrival  Rate  for  Synchronous  Traffic  (packets/msec', 

*  ' /station) : ' , f 9 .4, 

*  /,'  Mean  Packet  Length  for  Synchronous  Traffic  ', 

*  '  (msec) :',f9.4) 

759  format('  Asynchronous  Traffic:  Prioriry- j' , 7x, 'T_pri  aa(' 

*  ,Al,',j)  plena (',Al,',j)') 

705  format (33x, il, 3f 12 . 4) 
write (3, 65)  it 

65  format (/,'  Token  goes', i7,'  cycles  in  simulation') 

kk  -  1 
999  dm  -  0 . 
wm  «  0 . 
xm  -  0 . 
dv  -  0. 
wv  -  0 . 
xv  -  0. 
am  -  0. 

do  721  i  -  1,N 

if  (ifeat.eq.2 .and. a (i) .ne.kk)  goto  721 
if  (sf)  then 

if  (zcs (i) .eq.O)  then 
Ws(i)  -  -1 
Ds  (i)  -  -1 
else 

Ws(i)  ■  Ws(i)/zcs(i) 

Ds(i)  *  Ds(i)/zcs(i) 
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tqs(i)  -  tqs (i) /cyclen 

tqsv(i)  -  sqrt (tqsv (i) /cyclen-tqa (i) *tqs (i) ) 

XI  -  0. 

X2  -  0. 

do  67  j  -  0,mind 
XI  -  X1+ j*txs (i, j) 

X2  -  X2+j* j*txs (i, j) 
tds(i)  -  tds (i) +pds (i, j) 
tws(i)  -  tws (i) +pws <i, j) 
tzs(i)  *  tzs (i) +txs (i, j) 
tya(i)  -  tys (i) +tss (i, j) 
azs(i)  -  azs (i) +axs (i, j) 

67  dzs(i)  -  dzs (i) +dxs (i, j) 

Xsm(i)  =  Xl/tzs(i) 

Xav(i)  -  aqrt (X2/tzs (i) -Xam(i) *Xsm(i) ) 

am  -  am+aa (i) 

dm  »  dm+Ds (i) *as (i) 

wm  *  wm+Ws (i) *as  (i) 

xm  -  xm+Xsm(i) 

xv  -  xv+X2/tza(i) 

if  (zcs  (i)  .gt .  0)  then 

dv  «  dv+sigds (i) /zcs (i) *as  (i) 
vrv  -  wv+sigws (i) /zcs (i) *as (i) 
sigws(i)  -  sqrt (sigws (i) /zcs (i) -Ws (i) *Ws (i) ) 
sigds(i)  -  sqrt (sigds (i) /zcs (i) -Da (i) *Ds (i) ) 
end  if 
end  if 
721  continue 


* 


if  (ifeat.eq.2)  then 
m  -  NN(kk) 

write (3, 725)  '  **  Class-' //char (48+kk) 
else 
m  -  N 

write (3, 725)  '  ' 

end  if 
dm  -  dm/am 
wm  -  wm/am 
xm  -  xra/m 

725  format (/, all, 9x, 'E (X)  sigma (X) 6x, 'E (W)  sigma (W) 6x, 

*  'E(D)  sigma (D)') 

if  (af)  write (3, 722)  xm, aqrt (xv/m-xm*xm) , wm, sqrt (wv/am~wm*wra) , 

*  dm, sqrt (dv/am-dm*dm) 

722  formatt'  Sync  Traf f ic ! ' , 6f 10 . 4) 


do  723  j  -  l,Np 
am  -  0 . 
dm  -  0. 
wm  -  0. 
xm  -  0. 
dv  -  0 . 
wv  -  0 . 
xv  »  0. 

do  727  i  -  1,N 

if  (ifeat .eq.2 .and. a (i) .ne . kk)  goto  727 
tqa(i,j)  -  tqa (i, j) /cyclen 

tqav(i,j)  -  sqrt (tqav(i, j) /cyclen-tqa (i, j) *tqa (i, j) ) 
if  (zca (i, j) .eq.0)  then 
Wa(i,j)  -  -1 
Da  (i, j )  -  -1 
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Wa  (i,  j)  -  Wa  (i,  j)  /zca  (i,  j) 

Da  (i, j)  -  Da (i, j) /zca  (i, j) 
end  if 
XI  -  0. 

X2  -  0. 

do  764  k  -  0,mind 
XI  -  Xl+k*txa (i, j,k) 

X2  *  X2+k*k*txa <i, j, k) 
tda(i,j)  -  tda (i, j) +pda (i, j, k) 
twa (1/ j)  -  twa(i, j)+pwa(i, j,k) 
tza (i,  j)  -  tza(i, j)+txa(i, j,k) 
tya(i,j)  -  tya(i, j)+tsa(i, j,k) 
aza  (i, j)  -  aza (i, j) +axa (i, j, k) 
dza(i,j)  •»  dza  (i,  j) +dxa  (i,  j,  k) 

Xam(i,j)  -  Xl/tza(i,j) 

Xav(i,j)  =  sqrt  (X2/tza  (i,  j)  -Xam(i,  j)  *Xam(i,  j) ) 

am  -  am+aa (i, j) 

dm  -  dm+Da  (i,  j)  *aa  (i,  j) 

wm  -  wm+Wa  (i,  j)  *aa  (i,  j) 

xm  -  xm+Xam(i, j) 

if  (zca (i, j) .gt.O)  then 

dv  «  dv+sigda  (i,  j) /zca  (i,  j)  *aa  (i,  j) 
wv  -  wv+sigwa  (i,  j) /zca  (i,  j)  *aa  (i,  j) 

sigwa(i,j)  -  sqrt (aigwa (i, j) /zca (i, j) -Wa (i, j) *Wa (i, j) ) 
sigda(i,j)  *  sqrt (sigda (i, j) /zca (i, j) -Da (i, j) *Da (i, j) ) 
end  if 

xv  «  xv+X2/tza (i, j) 
continue 


dm  *  dm/am 
wm  -  wm/ am 
xm  *■  xm/m 

write  (3,724)  j,xm,sqrt  (xv/m-xm*xm)  ,wm,sqrt  (wv/am-wm*wm) , 
*  dm, sqrt (dv/am-dm*dm) 

724  format  ( '  Async:  pri-' , il, 6f 10 . 4) 

723  continue 

if  (ifeat .eq.2 .and.kk.eq. 1)  then 
kk  -  kk+1 
goto  999 
end  if 


write(3,70) 

70  format  (/,'  Station  #stat  E(X)  sigma (X)  E (W)  ', 

*' sigma (W)  E (D)  sigma <D) ',/, 75 ('*') ) 

do  68  i  -  1,N 

dwell (i)  -  dwell (i) /cyclen 

dwellv(i)  ■  sqrt (dwellv (i) /cyclen-dwell (i) *dwell (i) ) 
ctime(i)  -  ctime (i) /cyclen 

ctimev(i)  -  sqrt (ctiroev(i) /cyclen-ctime (i) *ctime (i) ) 
dpar(i)  -  dpar (i) /cyclen 

djparv(i)  -  sqrt  (dparv(i) /cyclen-dpar  (i)  *dpar  (i) ) 
if  (sf)  then 

write (3, 66)  i, 0, zcs (i) ,Xsm(i) ,Xsv(i) ,  Ws (i) ,  sigws (i)  ,Ds (i) , 
*  sigds(i) 

66  format  (/,  14,  i2, 18,  6f  10 . 4) 

*»r'd  if 


if  (.not.af)  goto  68 
do  968  j  -  l,Np 

if  (sf .or . j .gt . 1)  then 
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write (3, €9)  j, zca  (i, j) ,Xam(i, j) ,Xav(i, j) ,Wa (i, j) , 
sigwa (i, j) ,Da  <i, j) , sigda (i, j) 

else 

write (3, 66)  i, j, zca (i, j) , Xam(i, j) , Xav (i, j) , Wa (i, j) , 
sigwa (i, j) ,Da (i, j) , sigda (i, j) 

end  if 

format (5x, il,i8,6f!0.4) 


968  continue 


68  continue 


write (3,78)  W0,D0,X0 

78  format (/, '  Station  E(Vj)  sigma(Vj)  Pr(W>',f7.3, 

*')  Pr(D>',f7.3,')  Pr(X>',i2, ')',/, 73('*')) 

do  290  i  -  1,N 

if  (sf)  write<3/79)  i,  0,  tqs (i) , tqsv (i) , 1 . 0*PrWs (i) /zcs (i) , 

*  1.0‘PrDs (i) /zcs (i) , PrXs (i) /tzs (i) 
if  (.not.af)  goto  290 

do  291  j  ”  l,Np 

if  (sf .or. j .gt . 1)  then 

write (3, 279)  j,tqa(i, j) ,tqav(i, j) , 1 . 0*PrWa (i,  j) /zca(i,  j) , 

*  1 .0*PrDa  <i,  j)  /  zca  (i,  j)  ,PrXa  (i,  j)  /tza  (i,  j) 
else 

write (3,79)  i, j,tqa (i, j) , tqav(i, j) , 1 .0*PrWa (i, j) /zca (i, j) , 

*  1 .0*PrDa (i, j) /zca (i, j) ,PrXa (i, j) /tza (i, j) 
end  if 

291  continue 

79  format (/,i4,i2,2fl2.5,3fl4.5) 

279  format (i6, 2f 12 . 5, 3f 14 . 5) 

290  continue 


write (3, 765) 

765  format (/, '  Station' , 8x, ' E (V)  sigma (V) 8x, ' E (C)  sigma(C)', 

*  6x, 'E(C-V)  sigma (C-V)',/,  80 ('*')) 
do  768  i  -  1,N 

768  write(3,596)  i, dwell (i) ,dwellv (i) , ctime (i) , ctimev (i)  , 

*  dpar (i) ,dparv (i) 

596  format (i4, 4x, 6fl2 .7) 


do  91  i  -  1,N 
write (3, 571)  i 

571  format  (/,'  **  Station', i4) 
if  (sf)  then 

write(3,706)  as (i) *plens (i) 

706  format  (/,'  Synchronous  Traffic:', 

*  '  Normalized  throughput  (Rho_s) , f 9 . 6) 
write (3, 708) 

708  format (/,'  j  P(X-j)  P(Xa-j)  P(Xd-j)  P(Xt-j)',6x, 

*  'w  P(W<-w)',6x, 'd  P(D<-d)') 
g  -  0 . 

gl  -  0. 
g2  -  0. 

do  92  j  ■  0,mind 
g  ”  g+scale 

gl  -  gl+pws (i, j) /tws (i) 
g2  -  g2+pds (i, j) /tds (i) 

92  write (3, 93)  j,txs (i, j) /tzs (i) ,axs (i, j) /azs (i) ,dxs (i, j) /dzs (i) , 

*  tss (i, j) /tys (i),g,gl,g,g2 

93  format (i3,4f9.5,2(f9.4,f9.5) ) 

end  if 


if  (.not.af)  goto  91 
do  791  j  -  1, Np 

write  (3,726)  j,aa (i, j) ‘plena (i, j) 
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format  (/,'  Async:  priority-' , il, 5x, 

'  Normalized  throughput  (Rho_a) , f 9 . 6) 
write (3, 708) 
g  -  0. 
gl  -  0. 
g2  -  0. 

do  791  k  “  0,mind 
g  -  g+scale 

gl  -  gl+pwaU,  j,k)/twa(i,  j) 
g2  =  g2+pda (i, j,k) /tda (i, j) 

791  write (3,93)  k, txa (i, j, k) /tza (i, j) , axa (i, j, k) /aza (i, j) , 

*  dxa (i, j, k) /dza (i, j) ,  tsa (i, j, k) /tya (i, j) , g,gl, g,g2 

91  continue 


end 


subroutine  qwus (nq, t,  dif ) 

implicit  real*8  (A-H,0-Z) 

logical  zf,ofs (20) ,ofa (20, 3) 

real*8  Aas (20, 0 : 200) , works (20, 0 : 200) ,Us (20) , 

*  Aaa (20,3,0:100) , worka (20, 3, 0:100), Ua (20, 3) ,Ws(20) ,Ds(20) , 

*  cs (20) , ca (20,3) ,ys (20) ,ya (20, 3) , res (20) ,  v(20) 

*, sigws (20) , sigds (20) ,tofs(20) , t of a (20, 3) , wof s (20) , wofa (20,3) 
*,tss(20,0:50) ,tsa (20, 3, 0 : 50) , t jumps (20) , t jumpa (20, 3) 
integer  tails (20) , taila (20, 3) , ss (20) , sa (20, 3) , zcs (20) 

*, heads (20) ,heada (20,3) , seeds (20) , seeda (20, 3) , os (20)  ,oa (20,3) 

*,Xs (20) ,Xa(20,3) , axs (20,0:50) ,  dxs (20,0:50) , axa (20,3,0:50) 

*,pds (20,0:50) ,pws (20,0:50) ,PrWs (20) ,PrDs (20) 
common  Aas, Aaa, works, worka, Us, Ua, heads, heada, tails, taila, N,Np, 

*  Xs,Xa, axs, dxs, axa, mind, tas, tsa, t jumps, t jumpa 

common  /syn/v, res, Ns, Ds, sigws, sigds, scale, pws,pds, zcs,PrWs,PrDs, 

*  W0, DO  /bs/ofs,ss,tofs,cs,os,seeds,ys,wofs  /bc/sp,B,zf 
*/ba/ofa, sa, tofa, ca, oa, seeda, ya, wofa 

v (nq)  -  v (nq) +dif 
t2  •  t+dif 

ttl  Aas  (nq,  tails  (nq)  ) 
tt2  •  t+res (nq) 

2  if  ( (ttl) ,ge.t2 .and. 

*  tt2 .ge . (0 . 0000001 +t2) )  goto  4 
if  (ttl . It .tt2)  then 

if  (zf)  then 

if  (Xs (nq) .gt .mind)  then 

axs (nq, mind)  -  axs (nq, mind) +1 
tss(nq,mind)  -  tss (nq, mind) +ttl-t jumps (nq) 
else 

axs (nq,Xs (nq) )  -  axs (nq, Xs (nq) ) +1 
tss (nq, Xs (nq) )  ■  tss (nq,Xs (nq) ) +ttl-t jumps (nq) 
end  if 
end  if 

Xs (nq)  -  Xs(nq)+1 
t jumps (nq)  -  ttl 
j  -  mod (tails (nq) +1,201) 
if  ( j - eq . heads (nq) )  then 
temp  -  Aas (nq, tails (nq) ) 

6  if  (zf)  os(nq)  -  os(nq)+l 

si  “  dmod (sp*ss (nq) , B) 
pi  ■  -dlog (sl/B) *cs (nq) 
temp  ■  temp+max (0.000001, pi) 
ss (nq)  -  si 

if  (temp. It. t2)  goto  6 
goto  4 
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